hdu 4975 最大流及其唯一性判定（有向图环判断算法升级）

来源：互联网发布：linux查看所有进程编辑：程序博客网时间：2024/04/25 01:51

就当时最大流再次复习吧。。动手敲一下。。。经典解法不想说了。。这题主要是坑时间，10个提交7个tle。

环的判断，曾经用简单dfs方法，这次的就tle了！别人说要用很屌的dinic，我感觉自己dinic不可能超时，坚信是判断环慢了，于是学习了新断环的方法：删除点/边！从某点进去，若该点的所有边都遍历过还是无功而返，那么该店以后不用再进入了（这么简单的道感觉自己应该要想到啊！愚蠢啊！）开始时用只删除边，还是tle！nb!于是自己删点又删边，一下到156ms,前5了！

#include<cstdio>#include<iostream>#include<queue>#include<cstring>#include<string>using namespace std;const int maxv=1200;const int maxe=2*501*501+2000;const int inf=0x3f3f3f3f;int n,m;int allsumn=0,allsumm=0;int nume=0;int e[maxe][3];int head[maxv];bool flag;void inline adde(int i,int j,int c){    e[nume][0]=j;e[nume][1]=head[i];head[i]=nume;    e[nume++][2]=c;    e[nume][0]=i;e[nume][1]=head[j];head[j]=nume;    e[nume++][2]=0;}int lev[maxv];int vis[maxv];int ss=0;int tt=0;bool bfs(){    memset(lev,0,sizeof(lev));    memset(vis,0,sizeof(vis));    queue<int>q;    q.push(ss);    vis[ss]=1;    while(!q.empty())    {        int cur=q.front();        q.pop();        for(int i=head[cur];i!=-1;i=e[i][1])        {            int v=e[i][0];             if(e[i][2]>0&&!vis[v])             {                 lev[v]=lev[cur]+1;               // if(v==tt)return 1;          //这句不加，速度更快                 q.push(v);                 vis[v]=1;             }        }    }    return vis[tt];}int dfs(int u,int minf){    if(u==tt||minf==0)return minf;    int sumf=0,f;    for(int i=head[u];i!=-1&&minf;i=e[i][1])    {        int v=e[i][0];        if(lev[v]==lev[u]+1&&e[i][2]>0)        {            f=dfs(v,e[i][2]<minf?e[i][2]:minf);            sumf+=f;            e[i][2]-=f;e[i^1][2]+=f;            minf-=f;        }    }    if(!sumf)lev[u]=-1;                 return sumf;}int dinic(){    int sum=0;    while(bfs())      {        sum+=dfs(ss,inf);      }    return sum;}int hasv[maxv];void init(){    allsumm=allsumn=nume=0;    scanf("%d%d",&n,&m);    ss=n+m;tt=n+m+1;   for(int i=0;i<=tt;i++)   {       head[i]=-1;    hasv[i]=0;   }}void read_build(){    int temps=0;     for(int i=0;i<n;i++)          {              scanf("%d",&temps);              allsumn+=temps;              adde(ss,i,temps);              for(int j=0;j<m;j++)                adde(i,j+n,9);          }     for(int j=0;j<m;j++)         {            scanf("%d",&temps);             allsumm+=temps;             adde(j+n,tt,temps);         }}bool dfs_getother_ans(int u,int fa){    if(hasv[u])return 0;      for(int i=head[u];i!=-1;i=e[i][1])    {        int v=e[i][0];        if(v==fa||e[i][2]<=0||v==ss||v==tt)continue;        if(!vis[v])        {            vis[v]=1;            if(dfs_getother_ans(v,u))return 1;            vis[v]=0;                           //有向图的环出来的时候标记回来啊！        }        else        {            return 1;        }        e[i][2]=0;   //删除边！-50ms    }    hasv[u]=1;            //删除点，关键！-900ms    return 0;}int main(){     int T;    scanf("%d",&T);    for(int ii=1;ii<=T;ii++)    {       init();       read_build();       printf("Case #%d: ",ii);       if(allsumm!=allsumn)          {printf("So naive!\n");continue;}       int ans=dinic();       if(ans!=allsumm)        {            printf("So naive!\n");        }       else       {            flag=0;         for(int i=0;i<n;i++)         {              vis[i]=1;                 if(dfs_getother_ans(i,-1))                 {                     flag=1;                     break;                 }              vis[i]=0;         }         if(flag)           printf("So young!\n");         else         {             printf("So simple!\n");         }        }    }    return 0;}

0 0