bzoj5037 [Jsoi2014]电信网络

bzoj5037 [Jsoi2014]电信网络

​​ 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#include#include#include#includeusing 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(){ queueq;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;}

版权声明：本文内容由网络用户投稿，版权归原作者所有，本站不拥有其著作权，亦不承担相应法律责任。如果您发现本站中有涉嫌抄袭或描述失实的内容，请联系我们jiasou666@gmail.com 处理，核实后本网站将在24小时内删除侵权内容。