最大流 hdu 4183

/**
hdu 4183
求两点之间来回不重复的两条路就是要求从起点到终点至少有两条不重复的路径，最大流大于等于2即可
**/

#include <iostream>
#include <cstdio>
#include <cstring>
#include <vector>
#include <queue>
using namespace std;

struct edge{int to,cap,rev;};
const int INF=10000000000,max_v=1110;

vector<edge> G[max_v];
int level[max_v];
int iter[max_v];

void add_edge(int from,int to,int cap)
{
   G[from].push_back((edge){to,cap,G[to].size()});
   G[to].push_back((edge){from,0,G[from].size()-1});
}

void bfs(int s)
{
   memset(level,-1,sizeof(level));
   queue<int> que;
   level[s]=0;
   que.push(s);
   while(!que.empty())
   {
      int v=que.front();
      que.pop();
      for(int i=0;i<G[v].size();i++)
      {
         edge& e=G[v][i];
         if(e.cap>0 && level[e.to]<0)
         {
            level[e.to]=level[v]+1;
            que.push(e.to);
         }
      }
   }
}

int dfs(int v,int t,int f)
{
   if(v==t) return f;
   for(int &i = iter[v];i<G[v].size();i++)
   {
      edge &e=G[v][i];
      if(e.cap>0 && level[v]<level[e.to])
      {
         int d=dfs(e.to,t,min(f,e.cap));
         if(d>0)
         {
            e.cap-=d;
            G[e.to][e.rev].cap+=d;
            return d;
         }
      }
   }
   return 0;
}

int max_flow(int s,int t)
{
   int flow=0;
   for(;;)
   {
      bfs(s);
      if(level[t]<0) return flow;
      memset(iter,0,sizeof(iter));
      int f;
      while((f=dfs(s,t,INF))>0)
      {
         flow+=f;
      }
   }
   return flow;
}

struct node
{
   double rate;
   int x,y,r;
}s[305];

int main()
{
   int n,i,ss,tt,j,T;
   scanf("%d",&T);
   while(T--)
   {
      scanf("%d",&n);
      for(i=0;i<700;i++) G[i].clear();
      for(i=0;i<n;i++)
      {
         scanf("%lf%d%d%d",&s[i].rate,&s[i].x,&s[i].y,&s[i].r);
         if(s[i].rate==400) ss=i;
         else if(s[i].rate==789) tt=i;
      }
      bool f=false;
      for(i=0;i<n;i++)
         for(j=0;j<n;j++)
         {
            if(s[j].rate-s[i].rate>0.0001)
            {
               int dis=(s[i].x-s[j].x)*(s[i].x-s[j].x)+(s[i].y-s[j].y)*(s[i].y-s[j].y);
               if(dis<=(s[i].r+s[j].r)*(s[i].r+s[j].r))
               {
                  add_edge(i,j,1);
                  if(s[i].rate==400 && s[j].rate==789) f=true;
               }
            }
         }
      if(max_flow(ss,tt)>=2) f=true;
      if(f) printf("Game is VALID\n");
      else printf("Game is NOT VALID\n");
   }
   return 0;
}