poj 3498(满足点容量限制的网络流)

题意：某个冰块上有a只企鹅，总共可以跳出去b只，问是否可能所有的企鹅都跳到某一块冰块上，输出所有的可能的冰块的编号。

由于每个点只能跳出去m只企鹅，所以要拆点假如不拆点，一个点到另一个点可能会跳多于m只企鹅通过拆点后u->u'间的容量来完成题目的要求（对点的一些限制）

建图：i->i+n 容量为m i+n->j容量为INF新建源点s，s->i的容量为i点企鹅的个数然后枚举汇点求最大流就可以判断某个点是否符合条件

Source Code

Problem: 3498  User: 1013101127 
Memory: 508K  Time: 1032MS 
Language: C++  Result: Accepted 

Source Code 
#include<iostream>
#include<cstdio>
#include<cstring>
#include<set>
#include<cstdio>
#include<string>
using namespace std;
const int maxn=505;
//const int maxm=300006;
const int inf=1<<30;
struct node
{
    int v,next;
    int val;
}s[maxn*maxn*2];

int level[maxn];//顶点的层次
int p[maxn];
int que[maxn*10];//BFS中用于遍历的顶点，DFS求增广中记录边
int out[10*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++;
}
int n,f,d;
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;
                }
            }
        }
        //topset=r;//记录原点的集合个数
        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;
                //cout<<ret<<endl;
                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;
}
struct Node
{
    double x,y;
}poi[150];

double dis(Node a,Node b)
{
    return (a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y);
}
double R;
int nn[150],mm[150];
int main()
{
    int cas;
    cin>>cas;
    while(cas--)
    {
        int sum=0;
        int ans=0;
        //cin>>n>>R;
        scanf("%d %lf",&n,&R);
        for(int i=1;i<=n;i++)
        {
             //cin>>poi[i].x>>poi[i].y;
             //cin>>nn[i]>>mm[i];
            scanf("%lf%lf",&poi[i].x,&poi[i].y);
            scanf("%d%d",&nn[i],&mm[i]);
            sum+=nn[i];
        }
        int start=0;
        int end;
        for(int i=1;i<=n;i++)
        {
            start=0;
            end=i;
            init();
            for(int j=1;j<=n;j++)
            {
                insert(start,j,nn[j]);
                insert(j,j+n,mm[j]);
            }
            for(int k=1;k<n;k++)
            {
                for(int t=k+1;t<=n;t++)
                {
                    if(dis(poi[k],poi[t])<=R*R)
                    {
                        insert(n+k,t,inf);
                        insert(n+t,k,inf);
                    }
                }
            }
            int tmp=max_flow(n*2+1,start,end);
            if(tmp==sum)
            {
                ans++;
               // cout<<end-1<<' ';
                printf("%d ",end-1);
            }
        }
        if(!ans)
        cout<<"-1"<<endl;
        else
        cout<<endl;
    }
    return 0;
}