hdu5735

很美妙的一题

官方题解 http://www.cnblogs.com/duoxiao/p/5777632.html

感觉有meet in middle的思想

 #include<bits/stdc++.h>

 using namespace std;
 typedef long long ll;
 const int mo=1e9+;
 struct way{int po,next;} e[];
 unsigned int mx[][],f[],tmp[][];
 int n,w[],fa[],p[],v[],ans,len;
 char s[];
 void add(int x,int y)
 {
     e[++len].po=y;
     e[len].next=p[x];
     p[x]=len;
 }

 void inc(int &a,int b)
 {
     a+=b;
     if (a>mo) a-=mo;
 }

 int calc(int a,int b)
 {
     if (s[]=='X') return a^b;
     else return a&b;
 }

 void get(int x)
 {
     if (x==) f[x]=;
     else f[x]=f[fa[x]]+(w[x]|w[fa[x]]);
     inc(ans,(ll)x*(f[x]+w[x])%mo);
     for (int i=p[x]; i; i=e[i].next) get(e[i].po);
 }

 void dfs(int x)
 {
     f[x]=;
     int a=w[x]>>,b=w[x]&;
     for (int i=; i<; i++)
         if (v[i]) f[x]=max(f[x],mx[i][b]+(calc(a,i)<<));
     inc(ans,(ll)x*(f[x]+w[x])%mo);
     v[a]++;
     for (int i=; i<; i++)
     {
         tmp[x][i]=mx[a][i];
         mx[a][i]=max(mx[a][i],f[x]+calc(b,i));
     }
     for (int i=p[x]; i; i=e[i].next) dfs(e[i].po);
     memcpy(mx[a],tmp[x],sizeof(tmp[x]));
     v[a]--;
 }

 int main()
 {
     int cas;
     scanf("%d",&cas);
     memset(mx,,sizeof(mx));
     while (cas--)
     {
         len=; memset(p,,sizeof(p));
         scanf("%d%s",&n,s);
         for (int i=; i<=n; i++) scanf("%d",&w[i]);
         for (int i=; i<=n; i++)
         {
             scanf("%d",&fa[i]);
             add(fa[i],i);
         }
         ans=;
         memset(v,,sizeof(v));
         if (s[]=='O') get();
         else dfs();
         for (int i=; i<=n; i++) p[i]=;
         printf("%d\n",ans);
     }
 }