hdu4975--dinic算法（1）

好像还能用isap和矩阵dp的，不过我没咋闹明白isap，分层能看懂后面就完全看不懂了，等今天晚上看懂了在写一个用isap做这道题的方法。嗯，立flag；
题意：给出 行数m和列数n，和各行各列的和，求这个矩阵可能的排列，不同的结果不同的输出，我都写到注释上面了。
方法：dinic，就是ek算法的增强版，先以源点为0，然后对整个网络进行分层，一层一层的分，分完层在进行dfs，跑一个while，不断的对残余网络分层，搜索，知道分层分不了了，也就是说到不了终点了就完事。
思路是不难理解的，见图的思路很值得学习，
注意点1 有一个弧优化，就是在进行dfs的时候，这样dfs就可以用while了，试想dfs的对树的搜索，如果用这个方法，就可以在达到目的之前对 左，中，右进行遍历。不然只会返回一个结果，
2 在对 残残残残残余网络之后见图 在判断 路径的大小，这样就会很省时间

#include <queue>#include <vector>#include <stdio.h>#include <string.h>#include <stdlib.h>#include <iostream>#include <algorithm>#include <functional>const int maxn=2000;using namespace std;typedef long long ll;struct edge{int to,cap,rev;edge(int a,int b,int c){to=a,cap=b,rev=c;}//分别是目的地，记录反向边的坐标，和容量;};vector<edge>G[maxn];vector<int>G1[maxn];bool vis[maxn];int cengci[maxn];int n,m;int itr[maxn];void addedge(int from,int to,int value){   G[from].push_back((edge){to,value,G[to].size()});    G[to].push_back((edge){from,0,G[from].size()-1});}void add_edge(int from,int to){   G1[from].push_back(to);}void bfs(int s){    queue<int>q;     q.push(s);     memset(cengci,-1,sizeof(cengci));     cengci[s]=0;//这和上面一行都是分层必要的初始化;     while(!q.empty())     {   int u=q.front();     q.pop();        for(int i=0;i<G[u].size();i++)         {   edge &s=G[u][i];             if(cengci[s.to]<0&&s.cap>0)               {cengci[s.to]=cengci[u]+1;//分层;               q.push(s.to);}         }     }}int dfs(int s,int t,int f){   //int flow;   if(s==t) return f;//dfs结束的条件   for(int &i=itr[s];i<G[s].size();i++)   {  edge &w=G[s][i];       if(w.cap>0&&cengci[w.to]>cengci[s])       {  int d=dfs(w.to,t,min(w.cap,f));            if(d>0)            {   w.cap-=d;                G[w.to][w.rev].cap+=d;//反向边，自己画一个二维数组就理解了。g就是二位数组。                return d;            }       }   }return 0;}int max_flow(int s,int t){  int INF=0x3f3f3f3f;int flow=0;int cas=1;   while(1)   {   bfs(s);     memset(itr,0,sizeof(itr));      if(cengci[t]<0)        {  //printf("%d!!! %d\n",t,cengci[t]);            //printf("%d!!!!\n",cas++);            return flow;}      int l;//=dfs(s,t,INF);       while(l=dfs(s,t,INF))//多次进行dfs，而不是一次分层然后一次dfs，         flow+=l;   }}bool judge_dfs(int x,int val){   vis[x]=1;    for(int i=0;i<G1[x].size();i++)    {   int t=G1[x][i];        if(t==val) continue;//不找反向边。反向边一定能成为环所以你在逗我？？        if(vis[t]) return 1;//确定已经找到环了，因为vis被预定就是说明是环了        if(judge_dfs(t,x)) return 1;//继续dfs，把后继节点的t变成左边的x，这样在下一个dfs中如果t==val就是说明        //他的后继节点是前驱结点，所以就是环咯。    }    vis[x]=0;    return 0;}bool judge(){   memset(vis,0,sizeof(vis));   for(int i=0;i<=m+n+1;i++)   {   for(int j=0;j<G[i].size();j++)         {edge &e=G[i][j];         if(e.cap>0)//你看你看，如果在残余网络中边还存在，那么就见图，新建的图已经找不到一条增广路径了555；          add_edge(i,e.to);}   }    for(int i=0;i<=n;i++)       if(judge_dfs(i,-1)) return 1;//第二个参数设定成多少都行，只要不在 点的范围里。。就是钦定的边，就是说，如果判断某个具体的点能不能       //构成环。    return 0;}int main(){    int T;    int a,b,t=1;    scanf("%d",&T);    while(T--)    {  scanf("%d%d",&n,&m);       for(int i=0;i<maxn;i++)       {  G[i].clear();          G1[i].clear();       }       int sum1=0,sum2=0,sum=0;        for(int i=1;i<=n;i++)          {scanf("%d",&a);         addedge(0,i,a);     sum+=a,sum1+=a;//把所有杭上的点给加了     //也就是左边的原点到左边一排的最大流。         }         for(int j=n+1;j<=m+n;j++)         {   scanf("%d",&a);            addedge(j,m+n+1,a);            sum2+=a;//右边的最大流;         }         if(sum1!=sum2)         {printf("Case #%d: So naive!\n",t++);//不能成立，数据有问题啊哈哈          continue;         }          for(int i=1;i<=n;i++)            for(int j=1+n;j<=m+n;j++)          {  addedge(i,j,9);          }         int ans=max_flow(0,m+n+1);          //printf("%d   %d*******\n",ans,sum);        if(ans!=sum)            printf("Case #%d: So naive!\n",t++);//最大流跑不到和，不能成立。        else        {     if(judge())             printf("Case #%d: So young!\n",t++);//好多种类方法。             //在这个图中             //这个图表示的意思 就是各个行和和列和相互制约关系的成立。             //而这个图中有环，就说明了把每个边权加一个少一个都是没关系的。这部就有多种情况了嘛。             else                printf("Case #%d: So simple!\n",t++);//只有一种情况，费了好大劲做出来。-。-好简单啊。        }    }    return 0;}