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;
}
 


  • 0
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值