March of the Penguins

poj3498:http://poj.org/problem?id=3498

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

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

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

 #include<iostream>
 #include<cstring>
 #include<algorithm>
 #include<cstdio>
 #include<queue>
 #include<cmath>
 #define INF 100000000
 using namespace std;
 const int N=;
 const int M=;
 struct Node{
    int v;
    int f;
    int next;
 }edge[M];
 int n,m,u,v,cnt,sx,ex;
 int head[N],pre[N];
 int ans[N],top,val[N];
 void init(){
     cnt=;
     memset(head,-,sizeof(head));
 }
 void add(int u,int v,int w){
     edge[cnt].v=v;
     edge[cnt].f=w;
     edge[cnt].next=head[u];
     head[u]=cnt++;
     edge[cnt].f=;
     edge[cnt].v=u;
     edge[cnt].next=head[v];
     head[v]=cnt++;
 }
 bool BFS(){
   memset(pre,,sizeof(pre));
   pre[sx]=;
   queue<int>Q;
   Q.push(sx);
  while(!Q.empty()){
      int d=Q.front();
      Q.pop();
      for(int i=head[d];i!=-;i=edge[i].next    ){
         if(edge[i].f&&!pre[edge[i].v]){
             pre[edge[i].v]=pre[d]+;
             Q.push(edge[i].v);
         }
      }
   }
  return pre[ex]>;
 }
 int dinic(int flow,int ps){
     int f=flow;
      if(ps==ex)return f;
      for(int i=head[ps];i!=-;i=edge[i].next){
         if(edge[i].f&&pre[edge[i].v]==pre[ps]+){
             int a=edge[i].f;
             int t=dinic(min(a,flow),edge[i].v);
               edge[i].f-=t;
               edge[i^].f+=t;
              flow-=t;
              if(flow<=)break;
         }

      }
       if(f-flow<=)pre[ps]=-;
       return f-flow;
 }
 int solve(){
     int sum=;
     while(BFS())
         sum+=dinic(INF,sx);
     return sum;
 }
 double d;
 int xx[N],yy[N],num[N];
 int main(){
     int T,sum;
     scanf("%d",&T);
     while(T--) {
         sum=;
         scanf("%d %lf",&n,&d);
         for(int i=;i<=n;i++){
             scanf("%d%d%d%d",&xx[i],&yy[i],&val[i],&num[i]);
             sum+=val[i];
         }
         top=;
         for(int i=;i<=n;i++){
                 init();
             for(int j=;j<=n;j++){
                 if(j==i)continue;
                add(,j,val[j]);
                add(j,j+n,num[j]);
             }
             for(int j=;j<=n;j++){
                 for(int k=;k<=n;k++){
                     if(k==j)continue;
                     if(j==i)continue;
                     double temp=sqrt((xx[j]-xx[k])*1.0*(xx[j]-xx[k])+(yy[j]-yy[k])*1.0*(yy[j]-yy[k]));
                     if(temp<=d){
                         add(j+n,k,INF);
                     }
                 }
             }
             add(i,*n+,INF);
             sx=,ex=*n+;
            if(solve()==sum-val[i]){
                ans[++top]=i;
             }

         }
         if(top<)printf("-1\n");
         else{
             sort(ans+,ans+top+);
             for(int i=;i<top;i++)
                 printf("%d ",ans[i]-);
             printf("%d\n",ans[top]-);
         }
     }
       return ;
 }