bzoj5037 [Jsoi2014]电信网络

http://www.elijahqi.win/archives/3167
Description

JYY创建的电信公司，垄断着整个JSOI王国的电信网络。JYY在JSOI王国里建造了很多的通信基站。目前所有的基站
都是使用2G网络系统的。而现在3G时代已经到来了，JYY在思考，要不要把一些基站升级成3G网络的呢？JSOI王国
可以被看作为一个无穷大的二维平面，JYY一共建造了N个通信基站，第i个基站的坐标是(Xi,Yi)。每个基站有一个
通信范围Ri。第i号基站会向所有到其距离不超过Ri的基站发送信息。每个基站升级到3G网络都会有一个收益Si，
这个收益可能是正数（比如基站附近有个大城市，用户很多，赚的流量费也就很多了），也可能是负数（比如基站
周围市场不佳，收益不能填补升级基站本身的投资）。此外，由于原有的使用2G网络系统的基站无法解析从升级成
3G网络系统的基站所发来的信息（但是升级之后的基站是可以解析未升级基站发来的信息的），所以，JYY必须使
得，在升级工作全部完成之后，所有使用3G网络的基站，其通信范围内的基站，也都是使用3G网络的。由于基站数
量很多，你可以帮助JYY计算一下，他通过升级基站，最多能获得的收益是多少吗？

Input

第一行一个整数N；
接下来N行，每行4个整数，Xi，Yi，Ri，Si，表示处在(Xi,Yi)的基站的通信
范围是Ri，升级可以获得的收益是Si。
数据满足任意两个基站的坐标不同。
1≤N≤500，1≤Ri≤20000，|Xi|,|Yi|,|Si|≤10^4。

Output

输出一行一个整数，表示可以获得的最大收益。

Sample Input

5
0 1 7 10
0 -1 7 10
5 0 1 -15
10 0 6 10
15 1 2 -20
Sample Output

5
【样例说明】
我们可以将前三座基站升级成 3G 网络，以获得最佳收益。
HINT

Source

Round1

分析出来这是一个选这个另外一个必须选但是有负代价 那么就是最大权闭合子图

源向收益正的连边 收益负的向汇连边 如果在覆盖范围内就连边 跑最大流用总正费用减最大流即可

#include<queue>
#include<cstdio>
#include<cctype>
#include<cstring>
#include<algorithm>
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(!isdigit(ch)) {if (ch=='-') f=-1;ch=gc();}
    while(isdigit(ch)) x=x*10+ch-'0',ch=gc();
    return x*f;
}
const int inf=0x3f3f3f3f;
const int N=550;
struct node{
    int y,next,z;
}data[600000];
int num=1,h[N],level[N],cur[N],x[N],y[N],r[N],s[N],ans,n,T,dis[N][N];
inline void insert1(int x,int y,int z){
    data[++num].y=y;data[num].next=h[x];h[x]=num;data[num].z=z;
    data[++num].y=x;data[num].next=h[y];h[y]=num;data[num].z=0;
}
inline bool bfs(){
    queue<int>q;memset(level,0,sizeof(level));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;if (y==T) return 1;q.push(y);
        }
    }return 0;
}
inline int dfs(int x,int s){
    if (x==T) return s;int ss=s;
    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(z,s));if(!xx) level[y]=0;
            s-=xx;data[i].z-=xx;data[i^1].z+=xx;if(!s) return ss;
        }
    }return ss-s;
}
inline int sqr(int x){return x*x;}
inline int calc(int i,int j){
    return sqr(x[i]-x[j])+sqr(y[i]-y[j]);
}
int main(){
    freopen("bzoj5037.in","r",stdin);
    n=read();T=n+1;
    for (int i=1;i<=n;++i) x[i]=read(),y[i]=read(),r[i]=read(),s[i]=read(),ans+=s[i]>0?s[i]:0;
    for (int i=1;i<=n;++i)
        for (int j=i+1;j<=n;++j) dis[i][j]=dis[j][i]=calc(i,j);
    for (int i=1;i<=n;++i) if (s[i]>0) insert1(0,i,s[i]);else insert1(i,T,-s[i]);
    for (int i=1;i<=n;++i){static int tmp;tmp=sqr(r[i]);
        for (int j=1;j<=n;++j) {
            if (i==j) continue;
            if (dis[i][j]<=tmp) insert1(i,j,inf);
        }
    }while(bfs()) memcpy(cur,h,sizeof(h)),ans-=dfs(0,inf);
    printf("%d\n",ans);
    return 0;
}