Cerc2014 parades

題目大意：一個樹，每個數的度數，給定m條path，找出最多不共邊的path，不需要方案。  n<=1000,m<=n*(n-1)/2

 

這題就是一個treedp，先定一個根，然後從下到上dp就行了

對於一個點x，要經過x且在x子樹中的路徑要在不破壞兒子子樹內的路徑爲基礎，因爲如果你破壞的話，也只能多一條路徑，而可能會破壞多條

於是先求出每個兒子的子樹中自己可以完成的最多的路徑

然後再考慮某個可以上來的點到x的路徑

再考慮可以上來的點的兩兩配對的路徑

最後再求可以從x上去的點就OK了

 

 具體看代碼

#include<cmath>
#include<cstdio>
#include<cstring>
#include<cassert>
#include<iostream>
#include<algorithm>
using namespace std;
const intmaxn=2011,maxm=maxn<<1,maxd=10,maxchar=20000000;
int pre[maxm],son[maxm],now[maxn],tot; intlowbit[1<<maxd];
void add(int a,int b){pre[++tot]=now[a]; now[a]=tot;son[tot]=b;}
char buf[maxchar],*wbuf=buf;
inline void read(int &x){
               x=0;while(!isdigit(*wbuf)) ++wbuf;
               while(isdigit(*wbuf)) x=x*10+*wbuf++-'0';
}
void init(){
               int n;read(n); tot=0; memset(now,0,sizeof(now)); int a,b;
               for (inti=1;i<n;++i) read(a),read(b),add(a,b),add(b,a);
}
int up[maxn][maxn],match[1<<maxd]; boolpr[maxn][maxn],link[maxd][maxd];
int dfs(int x,int fa=0){
               intans=0,sum=0,s[maxd];
               for (intp=now[x];p;p=pre[p]) if (son[p]!=fa) ans+=dfs(son[p],x),s[sum++]=son[p];
               for (inti=0;i<sum;++i) for (int j=1;j<=up[s[i]][0];++j)
                                              if(pr[up[s[i]][j]][x]) ++ans,up[s[i]][0]=0;
               memset(link,0,sizeof(link));
               for (inti=0;i<sum;++i) for (int j=i+1;j<sum;++j){
                               intx=s[i],y=s[j];
                               for(int a=1;a<=up[x][0];++a) for (int b=1;b<=up[y][0];++b)
                                              if(pr[up[x][a]][up[y][b]]) {link[i][j]=link[j][i]=1; goto end;}
                               end:;
               }
               match[0]=0;int m=(1<<sum)-1;
               for (inti=1;i<=m;++i){
                               inta=lowbit[i]; match[i]=match[i-(1<<a)];
                               for(int b=a+1;b<sum;++b) if (link[a][b] && (i>>b)&1)
                                              match[i]=max(match[i],match[i-(1<<a)-(1<<b)]+1);
               }
               ans+=match[m];up[x][0]=1; up[x][1]=x;
               for (inti=0;i<sum;++i) if (match[m]==match[m-(1<<i)])
                               for(int t=1;t<=up[s[i]][0];++t) up[x][++up[x][0]]=up[s[i]][t];
               return ans;
}
void work(){
               int m,a,b;read(m); memset(pr,0,sizeof(pr));
               for (inti=1;i<=m;++i) read(a),read(b),pr[a][b]=pr[b][a]=1;
               printf("%d\n",dfs(1));
}
int main(){
               intm=1<<maxd; for (int i=1;i<m;++i)lowbit[i]=(i&1)?0:lowbit[i>>1]+1;
               fread(buf,1,maxchar,stdin);int T; read(T); while (T--) init(),work();
               return 0;
}