BZOJ 1822[JSOI2010]Frozen Nova 冷冻波

网络流+二分。

n^3枚举判断每个巫妖可以攻击的精灵，向其连1的边，每个精灵向汇点连1的边。

二分答案，修改源点流向每个巫妖的cap，跑最大流看是否等于精灵数。

恩，看起来没什么毛病。

然后狂WA不止。调了一晚上。拍了大半晚上，发现网上找的拿来拍的程序是个WA的。。。我还能说些什么呢。。

这时候才发现我应该算点到线段的距离而不是直线。保持微笑。

原来这题还有一个计算几何的tag?

算点到直线的距离d,点到线段两点的距离，短的为l,长的为r。

勾股定理算出tmp=sqrt(r*r-d*d);若是tmp小于线段长度，则返回d,否则返回l;

//Twenty
#include<cstdio>
#include<cstdlib>
#include<iostream>
#include<algorithm>
#include<cmath>
#include<cstring>
#include<queue>
#include<vector>
#define INF 0x3f3f3f3f
typedef long long LL;
const int maxn=**+;
using namespace std;
int ans,s,t,n,u,v,w,ecnt=,fir[maxn],d[maxn],cur[maxn],c[maxn],p[maxn],ed[maxn];
struct wuyao{
    int x,y,r,t;
}wy[maxn],jl[maxn],sm[maxn];
struct edge {
    int from,to,cap,flow,nxt;
    edge(){}
    edge(int from,int to,int cap,int flow,int nxt):from(from),to(to),cap(cap),flow(flow),nxt(nxt){}
}e[maxn];
void add(int u,int v,int w) {
    e[++ecnt]=edge(u,v,w,,fir[u]);
    e[++ecnt]=edge(v,u,,,fir[v]);
    fir[u]=ecnt-; fir[v]=ecnt;
}
queue<int>que;
void bfs(int s,int t) {
    for(int i=;i<=n;i++) d[i]=n;
    d[t]=;
    que.push(t);
    while(!que.empty()) {
        int x=que.front() ;que.pop();
        for(int i=fir[x];i;i=e[i].nxt)
        if(d[e[i].to]==n&&e[i].flow==e[i].cap) {
            d[e[i].to]=d[x]+;
            que.push(e[i].to);
        }
    }
}
int cal(int s,int t) {
    int fl=INF;
    for(int x=t;x!=s;x=e[p[x]].from)
        fl=min(fl,e[p[x]].cap-e[p[x]].flow);
    for(int x=t;x!=s;x=e[p[x]].from) {
        e[p[x]].flow+=fl;
        e[p[x]^].flow-=fl;
    }
    return fl;
}
int maxflow(int s,int t) {
    bfs(s,t);
    int res=;
    for(int i=;i<=n;i++) cur[i]=fir[i],c[d[i]]++;
    for(int x=s;d[x]<n;) {
        if(x==t) {
            res+=cal(s,t);
            x=s;
        }
        int ok=;
        for(int &i=cur[x];i;i=e[i].nxt)
            if(d[e[i].to]+==d[x]&&e[i].cap>e[i].flow){
                p[x=e[i].to]=i;
                ok=; break;
            }
        if(!ok) {
            cur[x]=fir[x]; int M=n;
            for(int i=cur[x];i;i=e[i].nxt)
                if(e[i].cap>e[i].flow)
                    M=min(M,d[e[i].to]+);
            if(!(--c[d[x]])) break;
            c[d[x]=M]++;
            if(x!=s) x=e[p[x]].from;
        }
    }
    return res;
}

double dis(int x,int y,int xx,int yy) {
    return sqrt((double)(x-xx)*(x-xx)+(double)(y-yy)*(y-yy));
}

double Dis(int x,int y,double A,double B,double C) {
    return fabs((A*x+B*y+C))/sqrt(A*A+B*B);
}

int wys,jls,sms;
double yyj;
int check(int ti) {
    for(int i=;i<=ecnt;i++) e[i].flow=;
    for(int i=;i<=wys;i++) {
        int fl=ti/wy[i].t+;
        e[ed[i]].cap=fl;
    }
    return (maxflow(s,t)==jls);
}

int main()
{
    //freopen(".in","r",stdin);
    //freopen(".out","w",stdout);

    scanf("%d%d%d",&wys,&jls,&sms);
    n=wys+jls+; s=n-; t=n;
    for(int i=;i<=wys;i++) {
        scanf("%d%d%d%d",&wy[i].x,&wy[i].y,&wy[i].r,&wy[i].t);
        add(s,i,); ed[i]=ecnt-;
    }
    for(int i=;i<=jls;i++) {
        scanf("%d%d",&jl[i].x,&jl[i].y);
        add(wys+i,t,);
    }
    for(int i=;i<=sms;i++)
        scanf("%d%d%d",&sm[i].x,&sm[i].y,&sm[i].r);    

    for(int i=;i<=wys;i++) {
        for(int j=;j<=jls;j++)
            if((yyj=dis(wy[i].x,wy[i].y,jl[j].x,jl[j].y))<=(double)wy[i].r) {
                double A=(wy[i].y-jl[j].y),B=jl[j].x-wy[i].x,C=wy[i].x*jl[j].y-wy[i].y*jl[j].x;
                  if(i==) {
                   int debug=;
                  }
                if(!sms)
                    add(i,wys+j,);
                for(int k=;k<=sms;k++) {
                    double tmp;
                    double ddx=Dis(sm[k].x,sm[k].y,A,B,C);
                    double zb=dis(wy[i].x,wy[i].y,sm[k].x,sm[k].y),yb=dis(jl[j].x,jl[j].y,sm[k].x,sm[k].y);
                    if(zb<yb) swap(zb,yb);
                    double woc=sqrt(zb*zb-ddx*ddx);
                    tmp=woc<=yyj?ddx:yb;
                    if(tmp<=(double)sm[k].r) break;
                    if(k==sms)
                        add(i,wys+j,);
                }
            }
        }

    int l=,r=4e6+;
    if(!check(r)) ans=-;
    else {
        while(l<=r) {
            int mid=(l+r)>>;
            if(check(mid)) ans=mid,r=mid-;
            else l=mid+;
        }
    }
    printf("%d\n",ans);
    return ;
}

题面