BZOJ4519: [Cqoi2016]不同的最小割

清空时没清空好

最小割树=分治+最小割
在当前集合随便找两个点，算他们之间的最小割
然后找当前最小割分割的两个集合（包含不在当前分治集合的点），就从起点遍历，每次只走未满流的边即可
然后更新这两个点集之间的最小割（包含不在当前分治集合的点，这些点之间的最小割也要更新），然后把当前集合按照最小割分成两个集合，分治下去

code：

#include<set>#include<map>#include<deque>#include<queue>#include<stack>#include<cmath>#include<ctime>#include<bitset>#include<string>#include<vector>#include<cstdio>#include<cstdlib>#include<cstring>#include<climits>#include<complex>#include<iostream>#include<algorithm>#define ll long long#define inf 1e15using namespace std;const int maxn = 900;const int maxm = 110000;int n,m;int E[maxm][2];struct edge{    int y,nex;    ll c;    edge(){}    edge(const int _y,const ll _c,const int _nex){y=_y;c=_c;nex=_nex;}}a[maxm],e[maxm]; int len,fir[maxn];inline void ins(const int x,const int y,const ll c){    a[++len]=edge(y,c,fir[x]);fir[x]=len;    a[++len]=edge(x,c,fir[y]);fir[y]=len;}int st,ed;int h[maxn];queue<int>q;bool bfs(){    for(int i=1;i<=n;i++) h[i]=0;    h[st]=1; q.push(st);    while(!q.empty())    {        const int x=q.front(); q.pop();        for(int k=fir[x];k;k=a[k].nex)        {            const int y=a[k].y;            if(a[k].c&&!h[y]) h[y]=h[x]+1,q.push(y);        }    }    return h[ed]!=0;}ll dfs(const int x,const ll flow){    if(x==ed) return flow;    ll delta=0;    for(int k=fir[x];k;k=a[k].nex)    {        const int y=a[k].y;        if(a[k].c&&h[y]==h[x]+1)        {            ll minc=dfs(y,min(a[k].c,flow-delta));            a[k].c-=minc; a[k^1].c+=minc;            delta+=minc;        }        if(delta==flow) return delta;    }    if(!delta) h[x]=0;    return delta;}void clear(){    for(int i=1;i<=m;i++)    {        int ii=E[i][0],c=E[i][1]-a[ii].c;        a[ii].c+=c,a[ii^1].c-=c;    }}ll Flow(){    //clear();    ll re=0;    while(bfs()) re+=dfs(st,(ll)inf);    return re;}map<ll,int>reg;int ans;bool v[maxn];void mark(){    for(int i=1;i<=n;i++) v[i]=false;    q.push(st); v[st]=true;    while(!q.empty())    {        const int x=q.front(); q.pop();        for(int k=fir[x];k;k=a[k].nex)        {            const int y=a[k].y;            if(a[k].c&&!v[y]) v[y]=true,q.push(y);        }    }}int p[maxn],temp[maxn];void solve(int L,int R){    if(L==R) return;    st=p[L],ed=p[R];    for(int i=2;i<=len;i+=2) a[i].c=a[i^1].c=(a[i].c+a[i^1].c)/2ll;    ll tk=Flow();    if(reg.count(tk)==0) reg[tk]=1,ans++;    mark();    int l=L,r=R;    for(int i=L;i<=R;i++)    {        if(v[p[i]]) temp[l++]=p[i];        else temp[r--]=p[i];    }    for(int i=L;i<=R;i++) p[i]=temp[i];    solve(L,l-1);     solve(r+1,R);}int main(){    len=1; //    scanf("%d%d",&n,&m);    for(int i=1;i<=m;i++)    {        int x,y,c; scanf("%d%d%d",&x,&y,&c);        E[i][0]=len+1;        E[i][1]=c;        ins(x,y,c);    }    for(int i=1;i<=n;i++) p[i]=i;    solve(1,n);    printf("%d\n",ans);    return 0;}