[最小割] BZOJ2400: Optimal Marks

论文题。
二进制每位独立算，每个编号就只有 $01$ 两种。可以看作分成两个集合，用最小割模型解。
题目要求在边权和最小的前提下，还要保证编号和最小。这个只需要每次从 $T$ 出发倒着走，能到的点一定是在 $T$ 集合内，其他的都看作是 0 0 <script type="math/tex" id="MathJax-Element-176">0</script>，这样就是最小的。

#include<cstdio>
#include<cctype>
#include<cstring>
#include<queue>
#include<algorithm>
#define _X first
#define _Y second
#define mp(x,y) make_pair(x,y)
const int maxn=1005,maxe=10005;
typedef long long LL;
using namespace std;
inline char gc(){
  static char buf[100000],*p1=buf,*p2=buf;
  return p1==p2&&(p2=(p1=buf)+fread(buf,1,100000,stdin),p1==p2)?EOF:*p1++;
}
inline int getint(){
  char ch=gc(); int res=0,ff=1;
  while(!isdigit(ch)) ch=='-'?ff=-1:0, ch=gc();
  while(isdigit(ch)) res=(res<<3)+(res<<1)+ch-'0', ch=gc();
  return res*ff;
}
struct Edge{
  int from,to,cap,flow;
  Edge(int t1=0,int t2=0,int t3=0,int t4=0){ from=t1;to=t2;cap=t3;flow=t4; }
} Es[maxe];
int fir[maxn],nxt[maxe],tot=1;
void add(int x,int y,int z){
  Es[++tot]=Edge(x,y,z,0);
  nxt[tot]=fir[x]; fir[x]=tot;
  Es[++tot]=Edge(y,x,0,0);
  nxt[tot]=fir[y]; fir[y]=tot;
}
queue<int> que;
int d[maxn],S,T;
bool Bfs(){
  memset(d,63,sizeof(d)); int INF=d[0];
  d[S]=0; que.push(S);
  while(!que.empty()){
    int x=que.front(); que.pop();
    for(int j=fir[x];j;j=nxt[j])
      if(Es[j].flow<Es[j].cap&&d[Es[j].to]==INF){
    d[Es[j].to]=d[x]+1; que.push(Es[j].to);
      }
  }
  return d[T]!=INF;
}
int fir_pos[maxn];
int dfs_find(int x,int flw){
  if(x==T||flw==0) return flw;
  int t,res=0;
  for(int &j=fir_pos[x];j;j=nxt[j])
    if(d[x]+1==d[Es[j].to]&&(t=dfs_find(Es[j].to,min(flw,Es[j].cap-Es[j].flow)))>0){
      res+=t; flw-=t; Es[j].flow+=t; Es[j^1].flow-=t;
      if(flw==0) return res;
    }
  return res;
}
int n,m,a[maxn],N;
int MaxFlow(){
  int res=0;
  while(Bfs()){
    for(int i=1;i<=N;i++) fir_pos[i]=fir[i];
    res+=dfs_find(S,1e9);
  }
  return res;
}
pair<int,int> e[maxn];
int id[maxn],fa[maxn];
int getfa(int x){ return fa[x]==x?x:fa[x]=getfa(fa[x]); }
int merge(int x,int y){
  x=getfa(x); y=getfa(y);
  if(x!=y) fa[x]=y;
}
LL ans,ans2,cnt;
bool vis[maxn];
void dfs(int x){
  vis[x]=true; cnt++;
  for(int j=fir[x];j;j=nxt[j])
    if(Es[j^1].cap>Es[j^1].flow&&!vis[Es[j].to]) dfs(Es[j].to);
}
int main(){
  freopen("bzoj2400.in","r",stdin);
  freopen("bzoj2400.out","w",stdout);
  n=getint(); m=getint();
  for(int i=1;i<=n;i++) a[i]=getint();
  for(int i=1;i<=m;i++) e[i]._X=getint(), e[i]._Y=getint();
  N=n; S=++N; T=++N;
  for(int s=0;s<=30;s++){
    tot=1; memset(fir,0,sizeof(fir));
    for(int i=1;i<=N-2;i++)
      if(a[i]>=0) if((a[i]>>s)&1) add(i,T,1e9); else add(S,i,1e9);
    for(int i=1;i<=m;i++) add(e[i]._X,e[i]._Y,1), add(e[i]._Y,e[i]._X,1); 
    ans+=(LL)(1<<s)*MaxFlow();
    memset(vis,0,sizeof(vis));
    cnt=0; dfs(T); 
    ans2+=(LL)(1<<s)*(cnt-1);
  }
  printf("%lld\n%lld\n",ans,ans2);
  return 0;
}