• poj 1637 Sightseeing tour——最大流+欧拉回路


    题目:http://poj.org/problem?id=1637

    先给无向边随便定向,如果一个点的入度大于出度,就从源点向它连 ( 入度 - 出度 / 2 ) 容量的边,意为需要流出去这么多;流出去1表示改了一条边的方向,会使自己出度-1、入度+1,所以容量要/2;出度大于入度的点类似地连向汇点;无向边按给它定的方向的反方向连上容量为1的边;最后看看能否满流即可。

    #include<cstdio>
    #include<cstring>
    #include<algorithm>
    #include<queue>
    using namespace std;
    const int N=205,M=1005,INF=N*M;
    int n,m,deg[N],hd[N],xnt,cur[N],sm[N],gn[N],gt[N],dfn[N];
    struct Ed{
      int to,nxt,cap;
      Ed(int a=0,int b=0,int c=0):to(a),nxt(b),cap(c) {}
    }ed[M+N<<1];
    queue<int> q;
    int rdn()
    {
      int ret=0;bool fx=1;char ch=getchar();
      while(ch>'9'||ch<'0'){if(ch=='-')fx=0;ch=getchar();}
      while(ch>='0'&&ch<='9')ret=ret*10+ch-'0',ch=getchar();
      return fx?ret:-ret;
    }
    void init()
    {
      xnt=1;memset(deg,0,sizeof deg);//xnt=1!!!
      memset(sm,0,sizeof sm);memset(hd,0,sizeof hd);
      memset(gn,0,sizeof gn);memset(gt,0,sizeof gt);
    }
    void add(int x,int y,int z)
    {
      ed[++xnt]=Ed(y,hd[x],z);hd[x]=xnt;
      ed[++xnt]=Ed(x,hd[y],0);hd[y]=xnt;
    }
    bool bfs()
    {
      memset(dfn,0,sizeof dfn);
      dfn[0]=1;q.push(0);
      while(q.size())
        {
          int k=q.front();q.pop();
          for(int i=hd[k],v;i;i=ed[i].nxt)
        if(!dfn[v=ed[i].to]&&ed[i].cap)
          dfn[v]=dfn[k]+1,q.push(v);
        }
      return dfn[n+1];
    }
    int dinic(int cr,int flow)
    {
      if(cr==n+1)return flow;//
      int use=0;
      for(int& i=cur[cr],v;i;i=ed[i].nxt)
        if(dfn[v=ed[i].to]==dfn[cr]+1&&ed[i].cap)
          {
        int tmp=dinic(v,min(flow-use,ed[i].cap));
        if(!tmp)dfn[v]=n+2;
        use+=tmp;ed[i].cap-=tmp;ed[i^1].cap+=tmp;
        if(use==flow)return use;
          }
      return use;
    }
    int main()
    {
      int T=rdn();
      while(T--)
        {
          n=rdn();m=rdn();init();
          for(int i=1,u,v,fx;i<=m;i++)
        {
          u=rdn();v=rdn();fx=rdn();
          deg[u]++;deg[v]++;sm[v]++;sm[u]--;//in - out
          if(!fx)add(v,u,1);
          else gn[v]++,gt[u]++;
        }
          bool flag=0;int val=0;
          for(int i=1;i<=n;i++)
        {
          if((deg[i]&1)||gn[i]>(deg[i]>>1)||gt[i]>(deg[i]>>1))
            {flag=1;break;}
          if(sm[i]>0)add(0,i,sm[i]>>1),val+=sm[i]>>1;
          else if(sm[i]<0)add(i,n+1,-sm[i]>>1);
        }
          if(flag){puts("impossible");continue;}
          int mxflow=0;
          while(bfs())
        {memcpy(cur,hd,sizeof hd);mxflow+=dinic(0,INF);}//fr:0 !!!
          if(mxflow==val)puts("possible");
          else puts("impossible");
        }
      return 0;
    }
  • 相关阅读:
    BF3,MW3,CF,高端?亲民
    关于#ifdef __cplusplus extern
    lua源码阅读顺序
    (ZZ)如何实现游戏主循环(Game Loop)的详细解析
    D3D学习总结基础篇(二)从古墓丽影的画面设置了解基础概念
    比较两个json是否相等
    IPAD点滴 WebIM
    Remoting与Font对象 WebIM
    使用android隐藏api实现亮度调节
    symbian的HTTP引擎中对302、301事件的处理
  • 原文地址:https://www.cnblogs.com/Narh/p/10117837.html
Copyright © 2020-2023  润新知