zoj 2587(判定が一意かどうか)


まず最大ストリームアルゴリズムで最小分割を求め、ソース点と集約点からそれぞれ遍歴し、遍歴したノードの総数が図中の点の個数に等しい場合、一意であり、そうでなければ一意ではない.
poj 3402と似ていますが、エッジは無方向です.
#include<iostream>
#include<cstdio>
#include<cstring>
#include<set>
#include<cstdio>
#include<string>
using namespace std;
const int maxn=1650;
const int inf=30000000;
struct node
{
    int v,next;
    int val;
}s[10050*4];

int level[maxn];//     
int p[maxn];
int que[200050];//BFS        ,DFS       
int out[maxn];//DFS         
int ind;

inline void insert(int x,int y,int z)
{
    s[ind].v=y;
    s[ind].val=z;
    s[ind].next=p[x];
    p[x]=ind++;
    s[ind].v=x;
    s[ind].val=0;
    s[ind].next=p[y];
    p[y]=ind++;
}

void init()
{
    ind=0;
    memset(p,-1,sizeof(p));
}

int max_flow(int n,int source,int sink)
{
    int ret=0;
    int h=0,r=0;
    while(1)//DFS
    {
        int i;
        for(i=0;i<=n;++i)
         level[i]=0;
         h=0,r=0;
        level[source]=1;
        que[0]=source;
        while(h<=r)//BFS
        {
            int  t=que[h++];
            for(i=p[t];i!=-1;i=s[i].next)
            {
                if(s[i].val&&level[s[i].v]==0)
                {
                    level[s[i].v]=level[t]+1;
                    que[++r]=s[i].v;
                }
            }
        }
        if(level[sink]==0)break;//     
        for(i=0;i<=n;++i)
        out[i]=p[i];
        int  q=-1;
        while(1)
        {
            if(q<0)
            {
                int  cur=out[source];
                for(;cur!=-1;cur=s[cur].next)
                {
                    if(s[cur].val&&out[s[cur].v]!=-1&&level[s[cur].v]==2)
                    {
                        break;
                    }
                }
                if(cur>=0)
                {
                    que[++q]=cur;
                    out[source]=s[cur].next;
                }
                else
                {
                    break;
                }
            }
            int  u=s[que[q]].v;
            if(u==sink)//     
            {
                int  dd=inf;
                int  index=-1;
                for(i=0;i<=q;i++)
                {
                    if(dd>s[que[i]].val)
                    {
                        dd=s[que[i]].val;
                        index=i;
                    }
                }
                ret+=dd;
                for(i=0;i<=q;i++)
                {
                    s[que[i]].val-=dd;
                    s[que[i]^1].val+=dd;
                }
                for(i=0;i<=q;i++)
                {
                    if(s[que[i]].val==0)
                    {
                        q=index-1;
                        break;
                    }
                }
            }
            else
            {
                long cur=out[u];
                for(;cur!=-1;cur=s[cur].next)
                {
                    if(s[cur].val&&out[s[cur].v]!=-1&&level[u]+1==level[s[cur].v])
                    {
                        break;
                    }
                }
                if(cur!=-1)
                {
                    que[++q]=cur;
                    out[u]=s[cur].next;
                }
                else
                {
                    out[u]=-1;
                    q--;
                }
            }
        }
    }
    return ret;
}

int n,m;
int cnts, cntt;
int flag1[maxn], flag2[maxn];

void dfs1(int v)
{
    flag1[v] = 1;
    cnts++;
    for(int i = p[v]; i != 0; i = s[i].next){
        if(flag1[s[i].v]==0 && s[i].val)
            dfs1(s[i].v);
    }
}
void dfs2(int v)
{
    flag2[v] = 1;
    cntt++;
    for(int i = p[v]; i != 0; i = s[i].next){
        if(flag2[s[i].v]==0 && s[i^1].val)
            dfs2(s[i].v);
    }
}
int main()
{
   int u,v,c;
   int start;
   int end;
   while(scanf("%d%d%d%d",&n,&m,&start,&end)!=EOF)
   {
    if(n==0&&m==0&&start==0&&end==0)
       break;
    init();
    for(int i=1;i<=m;i++)
    {
        cin>>u>>v>>c;
        insert(u,v,c);
        insert(v,u,c);
    }
        int maxx=max_flow(n+3,start,end);;
        memset(flag1,0,sizeof(flag1));
        memset(flag2,0,sizeof(flag2));
        cnts = cntt = 0;
        dfs1(start);
        dfs2(end);
        if(cnts+cntt!=n)
             printf("AMBIGUOUS
"); else printf("UNIQUE
"); } return 0; }