Cod sursa(job #4229)

Utilizator gigi_becaliGigi Becali gigi_becali Data 1 ianuarie 2007 16:50:58
Problema Distante Scor 0
Compilator cpp Status done
Runda Arhiva de probleme Marime 1.97 kb
#include <cstdio>
#include <map>
#include <algorithm>
#include <vector>
#include <string>
#include <deque>
#define oo 0x3f3f3f3f
#define left(i) (((i)<<1)+1)
#define right(i) (((i)<<1)+2)
#define tata(i) ((i-1)>>1)
#define maxn 1<<16
using namespace std;

int H[maxn], d[maxn], t[maxn], T, n, m, poz[maxn], d1[maxn];

int nh;
struct nod{ int nd, w; nod(int _nd, int _w){ nd=_nd; w=_w;};};

vector< vector<nod> >a;

void swap(int i, int j)
{
	int t;
	t=H[i]; H[i]=H[j]; H[j]=t;
	poz[H[i]]=i; poz[H[j]]=j;
}

void HeapDw(int i, int nh)
{
	int min=i;
	
	if(left(i)<=nh && d[H[left(i)]]<=d[H[i]]) min=left(i);
	if(right(i)<=nh && d[H[right(i)]]<=d[min]) min=right(i);
	if(min!=i) swap(min, i), HeapDw(min, nh);

}

void HeapUp(int k)
{
	int t;
	if(k<=0) return ;
	t=(k-1)>>1;
	if(d[H[k]]<=d[H[t]])
	{
		swap(k, t);
		HeapUp(t);
	}
}

void BuildH(int k)
{
	for(int i=1;i<k;i++) HeapUp(i);
}

int scoate_heap()
{
	swap(0, nh-1);
	poz[H[nh-1]]=0; nh--;
	HeapDw(0, nh-1);
	return H[nh];
}

void dijkstra(int s)
{
	int i, u;
	memset(d, oo, sizeof(d));
	nh=n;
	d[s]=0;
	for(i=0;i<n;i++) 
	{
		H[i]=i+1; poz[i+1]=i;
	}
	
	 BuildH(n); nh=n;
	
	while(nh>0)
	{
		u=scoate_heap();
		for(i=0;i<a[u].size();i++)
			if(d[a[u][i].nd]>d[u]+a[u][i].w)
			{
				d[a[u][i].nd]=d[u]+a[u][i].w;
				HeapUp(poz[a[u][i].nd]);
			}
	}
}



int main()
{
	freopen("distante.in", "r", stdin);
	freopen("distante.out", "w", stdout);
	int i, j,p, q, w, s;
	scanf("%d", &T);
	for(i=1;i<=T;i++)
	{	
		a.clear();
		scanf("%d %d %d\n", &n, &m, &s);
		a.resize(n+1);
		for(j=1;j<=n;j++) scanf("%d ", d1+j);
		for(j=1;j<=m;j++)
		{
			scanf("%d %d %d\n", &p, &q, &w);
			a[p].push_back(nod(q, w));
			a[q].push_back(nod(p, w));
		}
		dijkstra(s);
		int ok=1;
	//	for(j=1;j<=n;j++) printf("%d ", d[j]);
		//printf("\n");
		for(j=1;j<=n;j++) if(d1[j]!=d[j]){ ok=0; break;}
		if(ok) printf("DA\n");
			else printf("NU\n");
	}
	return 0;
}