[bzoj1316] 树上的询问

　　裸的点分治。。

　　及时把已经确定的询问清掉就能快不少。时间复杂度O(nlogn*p)

 1 #include<cstdio>
 2 #include<iostream>
 3 #include<cstring>
 4 #include<algorithm>
 5 using namespace std;
 6 const int maxn=10023;
 7 struct zs{int len,id;bool done;}q[103];
 8 struct zs1{int too,pre,dis;}e[maxn<<1];int tot,last[maxn];
 9 int st[maxn],top,pre,sz[maxn],mx[maxn],POI,dis[maxn];
10 bool gg[103],can[1000023],del[maxn];
11 int i,j,k,n,m,rt,mm;
12
13 int ra;char rx;
14 inline int read(){
15     rx=getchar(),ra=0;
16     while(rx<‘0‘||rx>‘9‘)rx=getchar();
17     while(rx>=‘0‘&&rx<=‘9‘)ra*=10,ra+=rx-48,rx=getchar();return ra;
18 }
19 inline int max(int a,int b){return a>b?a:b;}
20 void getrt(int x,int fa){
21     sz[x]=1,mx[x]=0;
22     for(int i=last[x];i;i=e[i].pre)if(!del[e[i].too]&&e[i].too!=fa)
23         getrt(e[i].too,x),sz[x]+=sz[e[i].too],mx[x]=max(mx[x],sz[e[i].too]);
24     mx[x]=max(mx[x],POI-sz[x]);
25     if(mx[x]<mx[rt])rt=x;
26 }
27 void getpoi(int x,int fa){
28     st[++top]=dis[x];
29     for(int i=last[x];i;i=e[i].pre)if(!del[e[i].too]&&e[i].too!=fa)
30         dis[e[i].too]=dis[x]+e[i].dis,getpoi(e[i].too,x);
31 }
32 inline void updQ(){
33     mm=m;m=0;
34     for(int i=1;i<=mm;i++)if(!q[i].done)q[++m]=q[i];
35 }
36 void work(int x){
37     int i,j,k;bool flag=0;
38     rt=0,getrt(x,0);
39     can[0]=1;pre=1;
40     for(i=last[rt];i;i=e[i].pre)if(!del[e[i].too]){
41         dis[e[i].too]=e[i].dis,getpoi(e[i].too,rt);
42         for(j=pre;j<=top;j++)
43             for(k=1;k<=m&&q[k].len>=st[j];k++)if(can[q[k].len-st[j]])q[k].done=1,flag=1;
44         while(pre<=top)
45             if(st[pre]<=q[1].len)
46                 can[st[pre++]]=1;
47             else pre++;
48     }
49     if(flag)updQ();
50     if(!m)return;
51     while(top)
52         if(st[top]<=1000000)can[st[top--]]=0;else top--;
53     del[rt]=1;
54     for(i=last[rt];i;i=e[i].pre)if(!del[e[i].too])POI=sz[e[i].too],work(e[i].too);
55 }
56
57 inline void insert(int a,int b,int c){
58     e[++tot].too=b,e[tot].dis=c,e[tot].pre=last[a],last[a]=tot;
59     e[++tot].too=a,e[tot].dis=c,e[tot].pre=last[b],last[b]=tot;
60 }
61 bool cmp(zs a,zs b){return a.len>b.len;}
62
63 int main(){mx[0]=1e9;
64     n=read(),m=read();int m1=m;
65     for(i=1;i<n;i++)j=read(),k=read(),insert(j,k,read());
66     for(i=1;i<=m;i++)q[i].len=read(),q[i].id=i;
67     sort(q+1,q+1+m,cmp);
68     while(m&&q[m].len==0)q[m].done=1,m--;
69     POI=n,work(1);
70     for(i=1;i<=m;i++)gg[q[i].id]=1;
71     for(i=1;i<=m1;i++)puts(gg[i]?"No":"Yes");
72     return 0;
73 }