bzoj 1822 [JSOI2010]Frozen Nova 冷冻波

辣鸡蒟蒻的原blog qwq
http://www.elijahqi.win/2018/01/04/bzoj-1822-jsoi2010frozen-nova-%e5%86%b7%e5%86%bb%e6%b3%a2/
Description
WJJ喜欢“魔兽争霸”这个游戏。在游戏中，巫妖是一种强大的英雄，它的技能Frozen Nova每次可以杀死一个小精灵。我们认为，巫妖和小精灵都可以看成是平面上的点。 当巫妖和小精灵之间的直线距离不超过R，且巫妖看到小精灵的视线没有被树木阻挡（也就是说，巫妖和小精灵的连线与任何树木都没有公共点）的话，巫妖就可以瞬间杀灭一个小精灵。 在森林里有N个巫妖，每个巫妖释放Frozen Nova之后，都需要等待一段时间，才能再次施放。不同的巫妖有不同的等待时间和施法范围，但相同的是，每次施放都可以杀死一个小精灵。 现在巫妖的头目想知道，若从0时刻开始计算，至少需要花费多少时间，可以杀死所有的小精灵？
Input
输入文件第一行包含三个整数N、M、K(N,M,K<=200)，分别代表巫妖的数量、小精灵的数量和树木的数量。 接下来N行，每行包含四个整数x, y, r, t，分别代表了每个巫妖的坐标、攻击范围和施法间隔（单位为秒）。 再接下来M行，每行两个整数x, y，分别代表了每个小精灵的坐标。 再接下来K行，每行三个整数x, y, r，分别代表了每个树木的坐标。 输入数据中所有坐标范围绝对值不超过10000，半径和施法间隔不超过20000。
Output
输出一行，为消灭所有小精灵的最短时间（以秒计算）。如果永远无法消灭所有的小精灵，则输出-1。
Sample Input
2 3 1

-100 0 100 3

100 0 100 5

-100 -10

100 10

110 11

5 5 10

Sample Output
5
HINT

Source
JSOI2010第二轮Contest1
此题好题orz 跪膜icefox 88ms过

蒟蒻我好菜啊
题意：给定n个巫妖 m个小精灵 每个巫妖每次可以干掉一个小精灵每个巫妖攻击之后有个冷却时间 然后询问最少的时间能够干掉全部的小精灵 那么好 把巫妖放在左边小精灵放在右边 可以感觉出来 因为每个小精灵只可以被干掉一次 所以每个小精灵向汇点连1的边 然后 每个巫妖 能干掉几个小精灵 我去二分一下即可 然后如果这个巫妖可以干掉这些小精灵 那么就他们之间连一条边 至于这个二分 有点意思 还是再次跪膜icefox 我大概想了很久没有想到 还是去膜了他的blog 就是我可以直接枚举我最长的那个人需要花费多少时间然后 每条边容量就是这个(最长时间/冷却时间)+1(因为0s的时候也可以攻击） 另外这个如何去判断巫妖能否攻击小精灵 相当于是他们之间不能有树木那就是看一个圆到线段是否有交点 解析几何我哪会啊orz 然后怎么搞 不能简单的使用点到直线距离公式 因为这是个线段啊
点到直线距离公式：练练latex(划去
$dis=\frac{\left | A*x0+B*y0+C\right|}{\sqrt{A^{2}+B^{2}}}$
这个直线是一般式 x0 y0就是树木的坐标 那么画个图可以知道我求一下点到直线的距离 然后再求一下点到线段每个端点的距离 去线段端点之间最长的距离作为斜边 然后 点到直线的距离为直角边 勾股定理求下另一个直角边 如果另一个直角边比线段距离长那么说明点到线段的最短距离是两端选一个最小的否则就是这个点到之间的距离 然后把这个距离和树木的半径比较一下 如果半径更大说明这个树挡住了精灵和巫妖
注意 有可能存在没有小精灵被干掉的情况提前预判
跪膜icefox 我在这里又gg了

#include<cmath>
#include<queue>
#include<cstdio>
#include<cstring>
#include<algorithm>
#define inf 0x3f3f3f3f
#define N 440
using namespace std;
inline char gc(){
    static char now[1<<16],*S,*T;
    if (T==S){T=(S=now)+fread(now,1,1<<16,stdin);if (T==S) return EOF;}
    return *S++;
}
inline int read(){
    int x=0,f=1;char ch=gc();
    while(ch<'0'||ch>'9') {if (ch=='-') f=-1;ch=gc();}
    while(ch<='9'&&ch>='0'){x=x*10+ch-'0';ch=gc();}
    return x*f;
}
int h[N],level[N],cur[N],T,t[N],num=1,st,n,m,in[N];
struct node{
    int x,y,z,next;
}data[88000];
inline void insert1(int x,int y,int z){
    data[++num].y=y;data[num].z=z;data[num].next=h[x];h[x]=num;data[num].x=x;
    data[++num].y=x;data[num].z=0;data[num].next=h[y];h[y]=num;data[num].x=y;
}
inline bool bfs(){
    memset(level,0,sizeof(level));queue<int>q;level[0]=1;q.push(0);
    while(!q.empty()){
        int x=q.front();q.pop();
        for (int i=h[x];i;i=data[i].next){
            int y=data[i].y,z=data[i].z;
            if (level[y]||!z) continue;level[y]=level[x]+1;q.push(y);if (y==T) return 1;
        }
    }return 0;
}
inline int dfs(int x,int s){
    if (x==T) return s;int ss=s,t=h[x];
    for (int i=cur[x];i;i=data[i].next){
        int y=data[i].y,z=data[i].z;
        if (level[x]+1==level[y]&&z){
            int xx=dfs(y,min(s,z));if(!xx) level[y]=0;t=i;
            s-=xx;data[i].z-=xx;data[i^1].z+=xx;if (!s) {cur[x]=i;return ss;}
        }
    }cur[x]=t;return ss-s;
}
inline double calc1(int x,int y,int x1,int y1){
    return sqrt((x-x1)*(x-x1)+(y-y1)*(y-y1));
}
inline double calc(int x1,int y1,int x2,int y2,int x3,int y3){
    double A=(double)(y1-y2)/(x1-x2);double B=-1;double C=y2-double(y1-y2)/(x1-x2)*x2;
    double dis3=fabs(A*x3+B*y3+C)/(sqrt(A*A+B*B));return dis3;
}
inline bool check(int mid){
    for (int i=2;i<st;++i) if (i&1) data[i].z=0;else data[i].z=1;
    for (int i=st;i<=num;++i){
        int now=i-st+1;now>>=1;now+=1;
        if (i&1) data[i].z=0;else data[i].z=mid/t[now]+1;
    }int ans=0;while(bfs()) memcpy(cur,h,sizeof(h)),ans+=dfs(0,inf);
    if (ans!=m) return 0;else return 1;
}
int x[N],y[N],xx1[N],yy1[N],x2[N],y2[N],r[N],r2[N],k;
int main(){
    freopen("bzoj1822.in","r",stdin);
    n=read();m=read();k=read();int max1=0;T=n+m+1;
    for (int i=1;i<=n;++i) x[i]=read(),y[i]=read(),r[i]=read(),t[i]=read(),max1=max(max1,t[i]);
    for (int i=1;i<=m;++i) xx1[i]=read(),yy1[i]=read();
    for (int i=1;i<=k;++i) x2[i]=read(),y2[i]=read(),r2[i]=read();
    for (int i=1;i<=n;++i){
        for (int j=1;j<=m;++j){
            bool flag=1;double dis=calc1(x[i],y[i],xx1[j],yy1[j]);
            if (dis>r[i]) continue;
            for (int z=1;z<=k;++z){
                double dis3=calc(x[i],y[i],xx1[j],yy1[j],x2[z],y2[z]);
                double dis1=calc1(x[i],y[i],x2[z],y2[z]),dis2=calc1(xx1[j],yy1[j],x2[z],y2[z]);
                if (dis1<dis2) swap(dis1,dis2);double dis4=sqrt(dis1*dis1-dis3*dis3);double min1=0;
                if (dis4<dis) min1=dis3;else min1=dis2;
                if (min1<=r2[z]) {flag=0;break;}
            }
            if (flag) insert1(i,n+j,1),++in[j];
        }
    }int l=0,r=m*max1;
    for (int i=1;i<=m;++i) if (!in[i]){puts("-1");return 0;}
    for (int i=1;i<=m;++i) insert1(n+i,T,1);st=num+1;
    for (int i=1;i<=n;++i) insert1(0,i,1);
    //for (int i=2;i<=num;++i) printf("%d %d %d\n",data[i].x,data[i].y,data[i].z);
    while(l<=r){
        int mid=l+r>>1;
        if (check(mid)) r=mid-1;else l=mid+1;
    }printf("%d",l);
    return 0;
}