bzoj 4519 [Cqoi2016]不同的最小割

用Gusfield算法建出Gomory-Hu树,之后统计有多少不同的最小割就可以了~

不了解Gusfield算法与Gomory-Hu树的话可以查查相关资料~

代码:

#include
#include
#include
#include
using namespace std;
typedef long long LL;
const LL INF=0x7ffffffffffffLL;//long long!!! 
struct front_star{LL ne,to,v;}a[100005];
LL h[100005],cnt=1,q[1000005],f[1005][1005],N,M,tot,dis[100005],Gap[100005];
LL vst[100005],prt[100005];
void Addedge(LL x,LL y,LL v)
{   a[++cnt].to=y;a[cnt].ne=h[x];a[cnt].v=v;h[x]=cnt;
    a[++cnt].to=x;a[cnt].ne=h[y];a[cnt].v=0;h[y]=cnt;
}
LL DFS(LL S,LL x,LL maxf,LL T)
{   if(x==T)return maxf;
    LL y,k,ret=0,dlt;
    for(k=h[x];k;k=a[k].ne)
      {y=a[k].to;
       if(a[k].v&&dis[x]==dis[y]+1)
         {dlt=DFS(S,y,min(maxf,a[k].v),T);
          a[k].v-=dlt;a[k^1].v+=dlt;
          maxf-=dlt;ret+=dlt;
          if(!maxf||dis[S]>=tot)return ret;
         }
      }
    if(!(--Gap[dis[x]]))dis[S]=tot;
    Gap[++dis[x]]++;
    return ret;
}
LL Sap(LL S,LL T)
{   LL i,ret=0;
    for(i=0;i<=N;i++)dis[i]=Gap[i]=0;
    Gap[0]=tot=N;
    while(dis[S]


#include
#include
#include
#include
using namespace std;
typedef long long LL;
const LL INF=0x7ffffffffffffLL;//long long!!! 
struct front_star{LL ne,to,v;}a[100005];
LL h[100005],cnt=1,q[1000005],f[1005][1005],N,M,tot,dis[100005],Gap[100005];
LL vst[100005],prt[100005];
void Addedge(LL x,LL y,LL v)
{   a[++cnt].to=y;a[cnt].ne=h[x];a[cnt].v=v;h[x]=cnt;
    a[++cnt].to=x;a[cnt].ne=h[y];a[cnt].v=0;h[y]=cnt;
}
LL DFS(LL S,LL x,LL maxf,LL T)
{   if(x==T)return maxf;
    LL y,k,ret=0,dlt;
    for(k=h[x];k;k=a[k].ne)
      {y=a[k].to;
       if(a[k].v&&dis[x]==dis[y]+1)
         {dlt=DFS(S,y,min(maxf,a[k].v),T);
          a[k].v-=dlt;a[k^1].v+=dlt;
          maxf-=dlt;ret+=dlt;
          if(!maxf||dis[S]>=tot)return ret;
         }
      }
    if(!(--Gap[dis[x]]))dis[S]=tot;
    Gap[++dis[x]]++;
    return ret;
}
LL Sap(LL S,LL T)
{   LL i,ret=0;
    for(i=0;i<=N;i++)dis[i]=Gap[i]=0;
    Gap[0]=tot=N;
    while(dis[S]


你可能感兴趣的:(网络流,各省省选)