题解:

树上Hash

首先重心在边上就把边分裂

以重心为根建树,这样两个根一定对应

然后f[i][j]表示i匹配另一棵的j节点的最小代价

把他们的儿子摘出来做最小权匹配即可

#include<iostream>
#include<cstdio>
#include<cstring>
#include<vector>
#include<queue>
#include<algorithm>
using namespace std;
const int maxn=1500;
typedef unsigned long long uLint;
const int oo=100000000;// int n;
int rx[maxn],ry[maxn];
int mi[maxn],ming[maxn]; int cntedge=0;
int head[maxn]={0};
int to[maxn<<1],nex[maxn<<1];
void Addedge(int x,int y){
nex[++cntedge]=head[x];
to[cntedge]=y;
head[x]=cntedge;
} int root1=0,root2=0;
int siz[maxn]={0},g[maxn]={0};
void Dfs(int x,int fa){
siz[x]=1;g[x]=0;
for(int i=head[x];i;i=nex[i]){
if(to[i]==fa)continue;
Dfs(to[i],x);
siz[x]+=siz[to[i]];
g[x]=max(g[x],siz[to[i]]);
}
g[x]=max(g[x],n-siz[x]); if(g[x]<=g[root1]){
root2=root1;root1=x;
}else if(g[x]<=g[root2]){
root2=x;
}
} uLint hh[maxn];
int dep[maxn];
int father[maxn];
void Gethash(int x,int fa){
hh[x]=666;
father[x]=fa;
dep[x]=dep[fa]+1;
for(int i=head[x];i;i=nex[i]){
if(to[i]==fa)continue;
Gethash(to[i],x);
hh[x]+=hh[to[i]];
}
hh[x]=hh[x]*hh[x];
} int a[maxn];
bool cmp(const int &rhs1,const int &rhs2){
if(dep[rhs1]==dep[rhs2])return hh[rhs1]<hh[rhs2];
else return dep[rhs1]>dep[rhs2];
} struct NetworkFlow{
int totn,s,t; struct Edge{
int from,to,cap,flow,cost;
};
vector<int>G[maxn];
vector<Edge>edges;
void Addedge(int x,int y,int z,int w){
Edge e;
e.from=x;e.to=y;e.cap=z;e.flow=0;e.cost=w;
edges.push_back(e);
e.from=y;e.to=x;e.cap=0;e.flow=0;e.cost=-w;
edges.push_back(e);
int c=edges.size();
G[x].push_back(c-2);
G[y].push_back(c-1);
} int inq[maxn];
int d[maxn];
int p[maxn];
queue<int>q;
int Spfa(int &nowflow,int &nowcost){
for(int i=1;i<=totn;++i){
d[i]=oo;inq[i]=0;
}
d[s]=0;inq[s]=1;q.push(s);
while(!q.empty()){
int x=q.front();q.pop();inq[x]=0;
for(int i=0;i<G[x].size();++i){
Edge e=edges[G[x][i]];
if((e.cap>e.flow)&&(d[x]+e.cost<d[e.to])){
d[e.to]=d[x]+e.cost;
p[e.to]=G[x][i];
if(!inq[e.to]){
q.push(e.to);
inq[e.to]=1;
}
}
}
} if(d[t]==oo)return 0; int f=oo,x=t;
while(x!=s){
Edge e=edges[p[x]];
f=min(f,e.cap-e.flow);
x=e.from;
}
nowflow+=f;nowcost+=f*d[t];
x=t;
while(x!=s){
edges[p[x]].flow+=f;
edges[p[x]^1].flow-=f;
x=edges[p[x]].from;
}
return 1;
} int Mincost(){
int flow=0,cost=0;
while(Spfa(flow,cost)){
}
return cost;
}
void MCMFinit(){
totn=n+n+2;s=n+n+1;t=s+1;
edges.clear();
G[s].clear();G[t].clear();
}
}W; int f[maxn][maxn];
int main(){
scanf("%d",&n);
for(int i=1;i<=n-1;++i){
int x,y;
scanf("%d%d",&x,&y);
Addedge(x,y);
Addedge(y,x);
rx[i]=x;ry[i]=y;
}
for(int i=1;i<=n;++i)scanf("%d",&ming[i]);
for(int i=1;i<=n;++i)scanf("%d",&mi[i]);
mi[n+1]=ming[n+1]=0; g[0]=0x7fffffff;
Dfs(1,0);
if(g[root1]!=g[root2])root2=0;
if(root1&&root2){
cntedge=0;
memset(head,0,sizeof(head));
for(int i=0;i<n;++i){
if((rx[i]==root1&&ry[i]==root2)||(rx[i]==root2&&ry[i]==root1))continue;
Addedge(rx[i],ry[i]);
Addedge(ry[i],rx[i]);
}
Addedge(root1,n+1);
Addedge(n+1,root1);
Addedge(root2,n+1);
Addedge(n+1,root2);
root1=++n;
}
Gethash(root1,0); for(int i=1;i<=n;++i)a[i]=i;
sort(a+1,a+1+n,cmp); for(int i=1;i<=n;++i){
for(int j=1;j<=n;++j){
int x=a[i];
int y=a[j];
if(hh[x]!=hh[y]){
f[x][y]=oo;
continue;
}
W.MCMFinit();
for(int ii=head[x];ii;ii=nex[ii]){
if(to[ii]==father[x])continue;
W.G[to[ii]].clear();
}
for(int ii=head[y];ii;ii=nex[ii]){
if(to[ii]==father[y])continue;
W.G[to[ii]+n].clear();
}
for(int ii=head[x];ii;ii=nex[ii]){
if(to[ii]==father[x])continue;
for(int jj=head[y];jj;jj=nex[jj]){
if(to[jj]==father[y])continue;
if(hh[to[ii]]!=hh[to[jj]])continue;
W.Addedge(to[ii],to[jj]+n,1,f[to[ii]][to[jj]]);
}
}
for(int ii=head[x];ii;ii=nex[ii]){
if(to[ii]==father[x])continue;
W.Addedge(W.s,to[ii],1,0);
}
for(int ii=head[y];ii;ii=nex[ii]){
if(to[ii]==father[y])continue;
W.Addedge(to[ii]+n,W.t,1,0);
}
f[x][y]=(ming[x]^mi[y])+W.Mincost();
}
} // for(int i=1;i<=n;++i){
// for(int j=1;j<=n;++j){
// cout<<f[i][j]<<' ';
// }
// cout<<endl;
// }
printf("%d\n",f[root1][root1]);
return 0;
}

  

BZOJ 3197 [Sdoi2013]assassin的更多相关文章

  1. bzoj 3197 [Sdoi2013]assassin(Hash+DP+KM)

    Description Input Output Sample Input 4 1 2 2 3 3 4 0 0 1 1 1 0 0 0 Sample Output 1 HINT [思路] Hash,D ...

  2. BZOJ 3197: [Sdoi2013]assassin 树形DP + 最小费用流 + 树的同构

    Description Input Output 其实就是给出两颗树,求一种两种树同构的方式,使得不同颜色个数最少$.$树的重新构建,其实就是指定不同的点为根节点$.$ 好在树的重心有一个重要的性质: ...

  3. 【BZOJ3197】[Sdoi2013]assassin 树同构+动态规划+KM

    [BZOJ3197][Sdoi2013]assassin Description Input Output Sample Input 4 1 2 2 3 3 4 0 0 1 1 1 0 0 0 Sam ...

  4. [BZOJ 3198] [Sdoi2013] spring 【容斥 + Hash】

    题目链接:BZOJ - 3198 题目分析 题目要求求出有多少对泉有恰好 k 个值相等. 我们用容斥来做. 枚举 2^6 种状态,某一位是 1 表示这一位相同,那么假设 1 的个数为 x . 答案就是 ...

  5. [BZOJ 3129] [Sdoi2013] 方程 【容斥+组合数取模+中国剩余定理】

    题目链接:BZOJ - 3129 题目分析 使用隔板法的思想,如果没有任何限制条件,那么方案数就是 C(m - 1, n - 1). 如果有一个限制条件是 xi >= Ai ,那么我们就可以将 ...

  6. BZOJ 3123: [Sdoi2013]森林 [主席树启发式合并]

    3123: [Sdoi2013]森林 题意:一个森林,加边,询问路径上k小值.保证任意时刻是森林 LCT没法搞,树上kth肯定要用树上主席树 加边?启发式合并就好了,小的树dfs重建一下 注意 测试点 ...

  7. BZOJ 3198: [Sdoi2013]spring [容斥原理 哈希表]

    3198: [Sdoi2013]spring 题意:n个物品6个属性,求有多少不同的年份i,j满足有k个属性对应相等 一开始读错题了,注意是对应相等 第i个属性只能和第i个属性对应 容斥一下 \[ 恰 ...

  8. 洛谷 P3307: bzoj 3202: [SDOI2013] 项链

    题目传送门:洛谷P3307.这题在bzoj上是权限题. 题意简述: 这题分为两个部分: ① 有一些珠子,每个珠子可以看成一个无序三元组.三元组要满足三个数都在$1$到$m$之间,并且三个数互质,两个珠 ...

  9. bzoj 3197 DP

    这道题我们可以看成给定两个黑白树,可以修改其中一棵树的颜色,问最少修改多少颜色可以使两棵树同构. 首先我们知道在树的同构中树上最长链中点(如果是偶数的话就是中间两个点)是不变的,我们把这个点叫做树的重 ...

随机推荐

  1. ubnutu18.4 修改交换分区大小

    ubuntu18.04默认的swap文件在根目录/下,名字是swapfile 1.查看交换分区大小  也可以使用系统监视器查看 free -m 2.创建一个swap文件 :大小为8g count= 3 ...

  2. python isdigit()函数

    isdigit() 函数是作用于字符串的,用来判断字符串是否全部由数字组成. x = '123456' y = 'iloveyou123' print(x.isdigit(),y.isdigit()) ...

  3. 15.Pythonic与python杂记

    switcher ={ :'sunday', :'monday', :'thuesday' } day = day_name=switcher.get(day,'Unknow') print(day_ ...

  4. 019.CI4框架CodeIgniter辅助函数类之:Array数组查询

    01. 数组辅助函数,可以方便的查看数组内部的成员,用法如下图所示: <?php namespace App\Controllers; class Hello extends BaseContr ...

  5. 修正png

    这是修正+取MD5的方法 function MD5FileTextPng(filename: AnsiString): AnsiString; var buf: ..MAX_PATH - ] of C ...

  6. Linux下安装JDK及其碰到的问题解决

    1.下载一个linux版本的jdk包 2.新建一个目录,专门用来存放安装包 mkdir    /home/software 3.将jdk包拷贝到/home/software下面,并解压 4.配置jdk ...

  7. 吴裕雄--天生自然java开发常用类库学习笔记:Arrays

    import java.util.* ; public class ArraysDemo{ public static void main(String arg[]){ int temp[] = {3 ...

  8. 吴裕雄 Bootstrap 前端框架开发——Bootstrap 字体图标(Glyphicons):glyphicon glyphicon-play-circle

    <!DOCTYPE html> <html> <head> <meta charset="utf-8"> <meta name ...

  9. P1039 到底买不买

    转跳点:

  10. 球队“食物链”(DFS+剪枝)

    某国的足球联赛中有N支参赛球队,编号从1至N.联赛采用主客场双循环赛制,参赛球队两两之间在双方主场各赛一场. 联赛战罢,结果已经尘埃落定.此时,联赛主席突发奇想,希望从中找出一条包含所有球队的“食物链 ...