Cod sursa(job #639139)

Utilizator indestructiblecont de teste indestructible Data 22 noiembrie 2011 16:01:42
Problema Guvern Scor 40
Compilator cpp Status done
Runda Arhiva de probleme Marime 1.91 kb
#include <stdio.h>
#include <string.h>
#include <algorithm>
#include <vector>
#include <set>
#define NMAX 200005
#define pii pair <int,int>
#define pb push_back
#define mp make_pair
#define f first
#define s second
using namespace std; 
int n,B[NMAX],dad[NMAX],st[NMAX],dr[NMAX],r,best[NMAX],E[NMAX],poz,rez;
vector <int> A[NMAX],C[NMAX];
vector < pair <pii,int> > D[NMAX];
set <pii> H,G;
set <pii> :: iterator it;
char viz[NMAX];
void read()
{
	scanf("%d",&n);
	int i,x,y;
	for (i=1; i<n; i++)
	{
		scanf("%d%d",&x,&y);
		A[x].pb(y);
		A[y].pb(x);
	}
	for (i=1; i<=n; i++)
		scanf("%d",&B[i]);
}
void dfs(int nod)
{
	viz[nod]=1;
	st[nod]=++r;
	it=H.lower_bound(mp(B[nod],0));
	if (it!=H.end())
		C[(*it).s].pb(nod);
	H.insert(mp(B[nod],nod));
	
	int i,vec;
	for (i=0; i<A[nod].size(); i++)
	{
		vec=A[nod][i];
		if (!viz[vec])
			dfs(vec);
	}
	H.erase(mp(B[nod],nod));
	dr[nod]=r;
}
void prepare()
{
	int i,j,vec;
	for (i=1; i<=n; i++)
	{
		for (j=0; j<C[i].size(); j++)
		{
			vec=C[i][j];
			D[i].pb(mp(mp(st[vec],dr[vec]),vec));
		}
		sort(D[i].begin(),D[i].end());
	}
}
inline int max(int x,int y)
{
	return x>y ? x : y;
}
void dfs2(int nod)
{
	viz[nod]=1;
	int i,j,vec,vec2;
	for (i=0; i<A[nod].size(); i++)
	{
		vec=A[nod][i];
		if (!viz[vec])
			dfs2(vec);
	}
	for (i=0; i<D[nod].size(); i++)
		E[i]=0;
	for (i=0; i<D[nod].size(); i++)
	{
		vec=D[nod][i].s;
		for (j=0; j<i; j++)
		{
			vec2=D[nod][j].s;
			if (dr[vec2]<st[vec])
				E[i]=max(E[i],E[j]);
		}
		E[i]+=best[vec];
	}
	for (i=0; i<D[nod].size(); i++)
		best[nod]=max(best[nod],E[i]);
	best[nod]++;
}
void find_sol()
{
	int i;
	for (i=1; i<=n; i++)
		rez=max(rez,best[i]);
}
int main()
{
	freopen("guvern.in","r",stdin);
	freopen("guvern.out","w",stdout);
	read();
	dfs(1);
	prepare();
	memset(viz,0,sizeof(viz));
	dfs2(1);
	find_sol();
	printf("%d\n",rez);
	return 0;
}