BZOJ1822[JSOI2010] 冷冻波

luogu链接：https://www.luogu.org/problemnew/show/P4048
bzoj链接：http://www.lydsy.com/JudgeOnline/problem.php?id=1822

冷冻波

题目描述

WJJ喜欢“魔兽争霸”这个游戏。在游戏中，巫妖是一种强大的英雄，它的技能Frozen Nova每次可以杀死一个小精灵。我们认为，巫妖和小精灵都可以看成是平面上的点。

当巫妖和小精灵之间的直线距离不超过R，且巫妖看到小精灵的视线没有被树木阻挡（也就是说，巫妖和小精灵的连线与任何树木都没有公共点）的话，巫妖就可以瞬间杀灭一个小精灵。

在森林里有N个巫妖，每个巫妖释放Frozen Nova之后，都需要等待一段时间，才能再次施放。不同的巫妖有不同的等待时间和施法范围，但相同的是，每次施放都可以杀死一个小精灵。

现在巫妖的头目想知道，若从0时刻开始计算，至少需要花费多少时间，可以杀死所有的小精灵？

输入输出格式

输入格式：

输入文件第一行包含三个整数N、M、K(N,M,K<=200)，分别代表巫妖的数量、小精灵的数量和树木的数量。

接下来N行，每行包含四个整数x, y, r, t，分别代表了每个巫妖的坐标、攻击范围和施法间隔（单位为秒）。

再接下来M行，每行两个整数x, y，分别代表了每个小精灵的坐标。

再接下来K行，每行三个整数x, y, r，分别代表了每个树木的坐标。

输入数据中所有坐标范围绝对值不超过10000，半径和施法间隔不超过20000。

输出格式：

输出一行，为消灭所有小精灵的最短时间（以秒计算）。如果永远无法消灭所有的小精灵，则输出-1。

输入输出样例

输入样例#1：

2 3 1
-100 0 100 3
100 0 100 5
-100 -10
100 10
110 11
5 5 10

输出样例#1：

5

题解

用计算几何求解巫妖和精灵间能否互相到达并以此建图，然后二分时间跑二分图匹配不断更新答案即可AC。

细节

注意先判断-1的情况，二分上界取maxcd*m，即冷却时间最长的巫妖单独杀死所有精灵的情况。更新时间后重新建图，超级源点到每个巫妖的边流量取time/cd[i]+1，即总时间除以冷却时间+1（因为在time=0的时候巫妖可以发射第一发Frozen Nova），其余边流量都为1。  

代码略长大家实现的时候耐心一点，其实就是计算几何+dinic没什么思维难度，主要考代码实现。

另外，洛谷上第三个点博主的裸dinic被卡T了，加个当前弧优化快了500ms，把inline去掉又快了200ms（迷。。。），bzoj无压力过。

代码

#include<bits/stdc++.h>
#define db double
#define R register
using namespace std;
const int M=50005;
const db eps=1e-8;
struct edge{int to,fl;};
struct pt{int x,y,r;};
int sig(db x){return (x>eps)-(x<-eps);}
db dis(pt a,pt b){return sqrt((a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y)*1.0);}
int operator * (pt a,pt b){return a.x*b.y-a.y*b.x;}
pt operator - (pt a,pt b){return (pt){a.x-b.x,a.y-b.y};}
db dis(pt p,pt a,pt b){return abs((a-p)*(b-p)*1.0/dis(a,b));}
pt lich[M],elv[M],tree[M];
edge ed[1000005];
vector<int>x[M];
int d[M],g=0,maxcd,cd[M],s,e,n,m,k,f;
bool vis[M];
char c;
void read(int &r)
{
    r=0;f=1;c=getchar();
    while(!isdigit(c))
    {
        if(c=='-')f=-1;
        c=getchar();
    }
    while(isdigit(c))r=(r<<1)+(r<<3)+c-'0',c=getchar();
    r*=f;
}
void add(int a,int b)
{
    ed[g]=(edge){b,1};
    x[a].push_back(g++);
    ed[g]=(edge){a,0};
    x[b].push_back(g++);
}
int bfs(int s,int e)
{
    R int i;
    memset(d,0,sizeof(d));
    queue<int>dui;
    dui.push(s);
    d[s]=1;
    int f,hh,t;
    while(!dui.empty())
    {
        f=dui.front();
        dui.pop();
        for(i=x[f].size()-1;i>=0;--i)
        {
            hh=x[f][i];
            t=ed[hh].to;
            if(d[t]||!ed[hh].fl)continue;
            d[t]=d[f]+1;
            dui.push(t);
        }
    }
    return d[e];
}
int dfs(int s,int e,int minn)
{
    R int i;
    if(s==e||!minn)return minn;
    int ans=0,hh,fl,to,tmp;
    for(i=x[s].size()-1;i>=0;--i)
    {
        hh=x[s][i];
        fl=ed[hh].fl;
        to=ed[hh].to;
        if(d[to]!=d[s]+1||!fl)continue;
        tmp=dfs(to,e,min(minn-ans,fl));
        if(!tmp)continue;
        ed[hh].fl-=tmp;
        ed[hh^1].fl+=tmp;
        ans+=tmp;
        if(minn==tmp)break;
    }
    return ans;
}
int dinic()
{
    int ans=0;
    while(bfs(e,s))
    ans+=dfs(e,s,1e9);
    return ans;
}
void reb(int t)
{
    R int i;
    int hh;
    for(i=0;i<g;i+=2)
    ed[i].fl=1,ed[i+1].fl=0;
    for(i=x[0].size()-1;i>=0;--i)
    ed[x[0][i]^1].fl=t/cd[ed[x[0][i]].to]+1;
}
void in()
{
    R int i;
    read(n);read(m);read(k);
    e=n+m+1;
    for(i=1;i<=n;++i)
    read(lich[i].x),read(lich[i].y),read(lich[i].r),read(cd[i]),maxcd=max(maxcd,cd[i]);
    for(i=1;i<=m;++i)
    read(elv[i].x),read(elv[i].y);
    for(i=1;i<=k;++i)
    read(tree[i].x),read(tree[i].y),read(tree[i].r);
}
bool check(pt a,pt b)
{
    R int i;
    if(sig(dis(a,b)-a.r)>0) return 0;
    for(i=1;i<=k;++i)
    {
        if(sig(dis(a,tree[i])-tree[i].r)<0)return 0;
        if(sig(dis(b,tree[i])-tree[i].r)<0)return 0;
        if((dis(tree[i],a,b)-tree[i].r)<0)return 0;
    }
    return 1;
}
void build()
{
    R int i,j;
    for(i=1;i<=n;++i)
    add(i,0);
    for(i=1;i<=n;++i)
    for(j=1;j<=m;++j)
    if(check(lich[i],elv[j]))
    {
        add(j+n,i);
        vis[j]=1;
    }
    for(i=1;i<=m;++i)
    add(n+m+1,n+i);
    for(i=1;i<=m;++i)
    if(!vis[i])
    {
        printf("-1");
        exit(0);
    }
}
void ac()
{
    int le=0,ri=m*maxcd,mid,ans,p;
    while(le<=ri)
    {
        mid=(le+ri)>>1;
        reb(mid);
        p=dinic();
        if(p==m)
        {
            ans=mid;
            ri=mid-1;
        }
        else le=mid+1;
    }
    printf("%d",ans);
}
int main()
{
    in();build();ac();
    return 0;
}