【HDOJ6598】Harmonious Army（最小割）

题意：有n个人，每个人可以从A,B两种职业中选择一种

有m对两人组，如果两个人都是A能获得p的收益，一个A一个B能获得q的收益，都是B能获得r的收益，其中q=p/4+r/3，保证p%4=0，r%3=0

求最大总收益

n<=5e2，m<=1e4，p,q,r<=4e6

思路：主要是建图

求得一组等效解，答案即为所有边权之和减去最小割

S出发到i的流量总数和i出发到T的流量总数可以累加一下最后再加这两种边，边数可以少一点

 #include<bits/stdc++.h>
 using namespace std;
 typedef long long ll;
 typedef unsigned int uint;
 typedef unsigned long long ull;
 typedef pair<int,int> PII;
 typedef pair<ll,ll> Pll;
 typedef vector<int> VI;
 typedef vector<PII> VII;
 typedef pair<ll,ll>P;
 #define N  100010
 #define M  200010
 #define fi first
 #define se second
 #define MP make_pair
 #define pi acos(-1)
 #define mem(a,b) memset(a,b,sizeof(a))
 #define rep(i,a,b) for(int i=(int)a;i<=(int)b;i++)
 #define per(i,a,b) for(int i=(int)a;i>=(int)b;i--)
 #define lowbit(x) x&(-x)
 #define Rand (rand()*(1<<16)+rand())
 #define id(x) ((x)<=B?(x):m-n/(x)+1)
 #define ls p<<1
 #define rs p<<1|1

 const ll MOD=1e9+,inv2=(MOD+)/;
       double eps=1e-;
       ll INF=1ll<<;
       ll inf=5e13;
       int dx[]={-,,,};
       int dy[]={,,-,};

 int head[N],vet[N],nxt[N],dis[N],gap[N],fan[N],s,S,T,tot;
 ll a[N][];
 double len[N];

 int read()
 {
    int v=,f=;
    char c=getchar();
    while(c<||<c) {if(c=='-') f=-; c=getchar();}
    while(<=c&&c<=) v=(v<<)+v+v+c-,c=getchar();
    return v*f;
 }

 void add(int a,int b,double c)
 {
     nxt[++tot]=head[a];
     vet[tot]=b;
     len[tot]=c;
     head[a]=tot;

     nxt[++tot]=head[b];
     vet[tot]=a;
     len[tot]=;
     head[b]=tot;
 }

 double dfs(int u,double aug)
 {
     if(u==T) return aug;
     int e=head[u],val=s-;
     double flow=;
     while(e)
     {
         int v=vet[e];
         if(len[e])
         {
             if(dis[u]==dis[v]+)
             {
                 double t=dfs(v,min(len[e],aug-flow));
                 len[e]-=t;
                 len[fan[e]]+=t;
                 flow+=t;
                 if(dis[S]>=s) return flow;
                 if(aug==flow) break;
             }
             val=min(val,dis[v]);
         }
         e=nxt[e];
     }
     if(!flow)
     {
         gap[dis[u]]--;
         if(!gap[dis[u]]) dis[S]=s;
         dis[u]=val+;
         gap[dis[u]]++;
     }
     return flow;
 }

 double maxflow()
 {
     rep(i,,s) gap[i]=dis[i]=;
     gap[]=s;
     double ans=;
     while(dis[S]<s) ans+=dfs(S,INF);
     return ans;
 }

 int main()
 {
     rep(i,,N-)
      if(i&) fan[i]=i+;
       else fan[i]=i-;
     int n,m;
     while(scanf("%d%d",&n,&m)!=EOF)
     {
         s=n;
         S=++s; T=++s;
         tot=;
         rep(i,,s) head[i]=;
         rep(i,,n) a[i][]=a[i][]=;
         ll ans=;
         rep(i,,m)
         {
             int x=read(),y=read(),p=read(),q=read(),r=read();
             a[x][]+=p+q;
             a[y][]+=p+q;
             a[x][]+=q+r;
             a[y][]+=q+r;
             add(x,y,(p+r-*q)/2.0);
             add(y,x,(p+r-*q)/2.0);
             ans+=p+q+r;
         }
         rep(i,,n)
         {
             add(S,i,a[i][]/2.0);
             add(i,T,a[i][]/2.0);
         }
         ans=(ll)ans-maxflow();
         printf("%I64d\n",ans);
     }

     return ;
 }