Cod sursa(job #256681)

Utilizator andrei-alphaAndrei-Bogdan Antonescu andrei-alpha Data 11 februarie 2009 23:47:41
Problema Kdrum Scor 80
Compilator cpp Status done
Runda Arhiva de probleme Marime 2.42 kb
using namespace std;

#include <set>
#include <map>
#include <list>
#include <deque>
#include <stack>
#include <queue>
#include <cstdio>

#define pb push_back
#define sz size
#define f first
#define s second
#define II inline
#define ll long long
#define db double
#define FOR(i,a,b) for(int i=a;i<=b;++i)
#define all(v) v.begin() , v.end()
#define CC(v) memset((v),0,sizeof((v)))
#define CP(v,w) memcpy((v),(w),sizeof((w)))
#define mp make_pair


typedef vector<int> VI;
typedef pair<int,int> pi;
typedef vector<string> VS;

#define IN "kdrum.in"
#define OUT "kdrum.out"
#define N_MAX 1<<6

int N,M,K,X1,Y1,X2,Y2;
int V[N_MAX][N_MAX];
int A[N_MAX][N_MAX][N_MAX];
int nr,stv[1<<8],poz[1<<14];

II void scan()
{
	freopen(IN,"r",stdin);
	freopen(OUT,"w",stdout);
	scanf("%d %d %d\n",&N,&M,&K);
	scanf("%d %d %d %d\n",&X1,&Y1,&X2,&Y2);
	FOR(i,1,N)
	FOR(j,1,M)
		scanf("%d",&V[i][j]);
}

II int gcd(int x, int y)  
{  
    for(int aux;y;)  
    {  
        aux = x%y;  
        x = y;  
        y = aux;  
    }  
    return x;  
}  

II void compute()
{
	int aux = K;
	aux = K;
	stv[0] = 1;
	FOR(i,2,K)
		if(!(aux % i))
		{	
			stv[++nr] = i;
			poz[i] = nr;
		}	
}

vector< pair<int,pi> > Q;

II void solve()
{
	int aux;
	memset(A,100,sizeof(A));
	A[X1][Y1][ poz[ gcd(V[X1][Y1],K) ] ] = 1;
	Q.pb( mp(poz[ gcd(V[X1][Y1],K) ],mp(X1,Y1)) );
	
	for(int ii=0;ii<(int)Q.sz();++ii)
	{
		int i = Q[ii].s.f;
		int j = Q[ii].s.s;
		int k = Q[ii].f;
		if(A[i][j][k] > 1<<20 || !V[i][j])
				continue;
		if(i == X2 && j == Y2 && k == poz[K])
			break;
		
		if(i!=1 && V[i-1][j])
		{
			aux = poz[gcd(stv[k] * V[i-1][j],K)];
			if(A[i-1][j][aux] > A[i][j][k]+1)
				A[i-1][j][aux] = A[i][j][k]+1,
				Q.pb( mp(aux,mp(i-1,j)) );
		}	
		if(i!=N && V[i+1][j])
		{
			aux = poz[gcd(stv[k] * V[i+1][j],K)];
			if(A[i+1][j][aux] > A[i][j][k]+1)
				A[i+1][j][aux] = A[i][j][k]+1,
				Q.pb( mp(aux,mp(i+1,j)) );
		}
		if(j!=1 && V[i][j-1])
		{	
			aux = poz[gcd(stv[k] * V[i][j-1],K)]; 
			if(A[i][j-1][aux] > A[i][j][k]+1)
				A[i][j-1][aux] = A[i][j][k]+1,
				Q.pb( mp(aux,mp(i,j-1)) );
		}
		if(j!=M && V[i][j+1])
		{
			aux = poz[gcd(stv[k] * V[i][j+1],K)];
			if(A[i][j+1][aux] > A[i][j][k]+1)
				A[i][j+1][aux] = A[i][j][k]+1,
				Q.pb( mp(aux,mp(i,j+1)) );
		}	
	}
	printf("%d\n",A[X2][Y2][poz[K]]);
}

int main()
{
	scan();
	compute();
	solve();
	return 0;
}