BZOJ 1822 [JSOI2010]Frozen Nova 冷冻波

很明显的一道最大流，利用计算几何判断是否被攻击的关系。二分时间以判断攻击次数，每次重新建图。

#include<cstdio>
#include<iostream>
#include<cstring>
#include<vector>
#include<queue>
#include<stack> 
#include<cmath>
using namespace std;
const int maxn = 10005;
const int inf = 0x3f3f3f3f;
struct edge
{
    int u,v,cap,flow;
};
struct P
{
    int x,y,r,t;
};
struct L{
    P a,b;
};
P w[maxn],s[maxn],t[maxn];
vector<int> G[maxn];
vector<edge> f;
void add_edge(int u,int v,int cap)
{
    f.push_back((edge){u,v,cap,0});
    f.push_back((edge){v,u,0,0});
    int m=f.size();
    G[u].push_back(m-2);
    G[v].push_back(m-1);
}
int d[maxn],cur[maxn],mx,x,y,z,n,m,S,T,l,r,k;
bool vis[maxn],map[505][505];
bool bfs()
{
    memset(vis,false,sizeof(vis));
    queue<int> Q;
    Q.push(S);
    vis[S]=true;d[S]=0;
    while(!Q.empty())
    {
        int u=Q.front();
        Q.pop();
        for(int i=0;i<G[u].size();i++)
        {
            edge& e=f[G[u][i]];
            if(e.cap>e.flow&&!vis[e.v])
            {
                d[e.v]=d[u]+1;
                vis[e.v]=true;
                Q.push(e.v);
            }
        }
    }
    return vis[T];
}
int dfs(int u,int a)
{
    int flow=0,h;
    if(u==T||a==0) return a;
    for(int &i=cur[u];i<G[u].size();i++)
    {
        edge& e=f[G[u][i]];
        if(d[e.v]==d[u]+1&&(h=dfs(e.v,min(e.cap-e.flow,a)))>0)
        {
            flow+=h;
            f[G[u][i]^1].flow-=h;
            e.flow+=h;
            a-=h;
            if(a==0) break; 
        }
    }
    return flow;
}
int Dinic(int s,int t)
{
    int flow=0;
    while(bfs()) {
        memset(cur,0,sizeof(cur));
        flow+=dfs(s,inf);
    }
    return flow;
}
inline P operator-(P a,P b)
{
    P t;t.x=a.x-b.x;t.y=a.y-b.y;return t;
}
inline double operator*(P a,P b)
{
    return a.x*b.y-a.y*b.x;
}
inline double dis(P a,P b)
{
    return sqrt((a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y));
}
inline double dot(P a,P b)
{
    return a.x*b.x+a.y*b.y;
}
inline double dis(L l,P p)
{
    if(dot(l.a-p,l.b-p)>0)
        return min(dis(l.a,p),dis(l.b,p));
    return abs((p-l.a)*(l.b-l.a)/dis(l.a,l.b));
}
bool jud(int x,int y)
{
    if(dis(w[x],s[y])>w[x].r)return 0;
    for(int i=1;i<=k;i++)
        if(dis((L){w[x],s[y]},t[i])<t[i].r)return 0;
    return 1;
}
void build(int x)
{
    f.clear();T=n+m+1; 
    for(int i=0;i<=T;i++) G[i].clear();
    for(int i=1;i<=n;i++) add_edge(0,i,x/w[i].t+1);
    for(int i=1;i<=m;i++) add_edge(i+n,T,1);
    for(int i=1;i<=n;i++)
        for(int j=1;j<=m;j++)
            if(map[i][j]) add_edge(i,j+n,1);
}
int main()
{
    scanf("%d%d%d",&n,&m,&k);
    for(int i=1;i<=n;i++) scanf("%d%d%d%d",&w[i].x,&w[i].y,&w[i].r,&w[i].t),mx=max(mx,w[i].t);
    for(int i=1;i<=m;i++) scanf("%d%d",&s[i].x,&s[i].y);
    for(int i=1;i<=k;i++) scanf("%d%d%d",&t[i].x,&t[i].y,&t[i].r);
    for(int i=1;i<=n;i++)
        for(int j=1;j<=m;j++)
        {
            map[i][j]=jud(i,j);
            if(map[i][j]) d[j]=1;           
        }
    for(int i=1;i<=m;i++)
    {
        if(!d[i]){
            puts("-1");return 0;
        }
    }
    l=0,r=mx*m;int mid;
    while(l<=r)
    {
        mid=(l+r)>>1;
        build(mid);
        int ans=Dinic(S,T);
        if(ans==m) r=mid-1;
        else l=mid+1;
    }
    printf("%d\n",l);
    return 0;   
}