题目描述:
定义无向图中的一条边的值为:这条边连接的两个点的值的异或值。
定义一个无向图的值为:这个无向图所有边的值的和。
给你一个有n个结点m条边的无向图。其中的一些点的值是给定的,而其余的点的值由你决定(但要求均为非负数),使得这个无向图的值最小。在无向图的值最小的前提下,使得无向图中所有点的值的和最小。
题目分析:
边的每个二进制位的贡献是独立的,所以可以按位处理
一条边同样是一个二元关系,两个点不同会产生1的代价
定义割完之后与S相连的点为1,与T相连的点为0,图就这样建:
- 当前位确定为1的点,从S向它连 ∞ \infty ∞的边
- 当前为确定为0的点,从T向它连 ∞ \infty ∞的边
- 如果原图有边(x,y),那么x,y之间连一条1的双向边
跑出来的最小割乘上对应的位贡献,累加就是第一问的答案
如何保证点权和最小?就是既可以选1又可以选0的时候优先选0。
从S到T的路径上有很多边都是满流的,最小割可以取其中的任意一条或多条,那么从S开始dfs,遇到满流边就停止,可以保证最小割是靠近S集合的,也就使得T集合的点尽量多(在T集合表示为0)
或者更为稳妥的做法是给代价加上一个第二关键字
具体做法是把原来的代价*10000,然后如果选1就附加上1的代价,连边就是从每个点向T连一条容量为1的边,割掉这条边的代价为1
然后第一问的答案就是最小割/10000,选1的点的个数就是最小割%10000
Code:
#include<cstdio>
#include<cctype>
#include<cstring>
#include<queue>
#include<algorithm>
#define LL long long
#define maxn 505
#define maxm 10005
using namespace std;
char cb[1<<15],*cs,*ct;
#define getc() (cs==ct&&(ct=(cs=cb)+fread(cb,1,1<<15,stdin),cs==ct)?0:*cs++)
inline void read(int &a){
char c;bool f=0;
while(!isdigit(c=getc())) if(c=='-') f=1;
for(a=c-'0';isdigit(c=getc());a=a*10+c-'0');
if(f) a=-a;
}
const int inf = 0x3f3f3f3f;
int n,m,S,T;
int fir[maxn],cur[maxn],dis[maxn],nxt[maxm],to[maxm],tot=1,cap[maxm],sum;
inline void line(int x,int y,int z,int rz=0){
nxt[++tot]=fir[x],fir[x]=tot,to[tot]=y,cap[tot]=z;
nxt[++tot]=fir[y],fir[y]=tot,to[tot]=x,cap[tot]=rz;
}
queue<int>q;
bool bfs()
{
memset(dis,0,(T+1)<<2);
dis[T]=1,q.push(T);
while(!q.empty()){
int u=q.front();q.pop();
for(int i=fir[u];i;i=nxt[i]) if(cap[i^1]&&!dis[to[i]]){
dis[to[i]]=dis[u]+1;
q.push(to[i]);
}
}
return dis[S];
}
int dfs(int u,int lim)
{
if(u==T) return lim;
int need=lim,delta;
for(int &i=cur[u];i;i=nxt[i]) if(cap[i]&&dis[u]==dis[to[i]]+1){
delta=dfs(to[i],min(cap[i],need));
cap[i]-=delta,cap[i^1]+=delta;
if(!(need-=delta)) break;
}
return lim-need;
}
int Dinic(){
int flow=0;
while(bfs()) memcpy(cur,fir,(T+1)<<2),flow+=dfs(S,inf);
return flow;
}
LL ans1,ans2;bool vis[maxn];
void calc(int u,int val){
if(u) ans2+=val;vis[u]=1;
for(int i=fir[u];i;i=nxt[i]) if(cap[i]&&!vis[to[i]]) calc(to[i],val);
}
int x[maxm],y[maxm],a[maxn];
int main()
{
read(n),read(m);S=0,T=n+1;
for(int i=1;i<=n;i++) read(a[i]);
for(int i=1;i<=m;i++) read(x[i]),read(y[i]);
for(int bit=0;bit<30;bit++){
memset(fir,0,(T+1)<<2),tot=1;
for(int i=1;i<=m;i++) line(x[i],y[i],1,1);
for(int i=1;i<=n;i++){
if(a[i]<0) continue;
if(a[i]&(1<<bit)) line(S,i,inf);
else line(i,T,inf);
}
ans1+=1ll*Dinic()*(1<<bit);
memset(vis,0,(T+1)),calc(S,1<<bit);
}
printf("%lld\n%lld",ans1,ans2);
}