分组赛时wy大神讲的题,网上都是随机化的题解

我来讲一下正解吧,我们穷举两个点,这两点距离要小于限制

然后我们分别以这两个点为圆心,两点距离为半径画圆

圆圆相交的部分被两点练成线段划分成两部分,不难发现

每个部分内点点之间的距离是小于限制的,很明显想到二分图

对于上半部分与下半部分的两点,如果距离大于限制则连边

然后我们求最大点独立集即可

求最大独立集的方案各种想错,实在太SB

做完最大匹配后首先先把未匹配的点拉出来,然后把他们有边相连的点都标记

然后处理每对匹配,若一个点被标记则选另一个点

这样是不可能出现两个点都标记的匹配,因为这样就会出现一条增广路,与最大匹配矛盾

type node=record
po,next:longint;
end; var e,ee:array[..] of node;
q,q1,q2,c,x,y,cy,cx,p:array[..] of longint;
v,v1,v2:array[..] of boolean;
z,ll,d,t1,t2,i,j,k,l,ans,len,n:longint; procedure add(x,y:longint);
begin
inc(len);
e[len].po:=y;
e[len].next:=p[x];
p[x]:=len;
ee[len].po:=x;
ee[len].next:=q[y];
q[y]:=len;
end; function dfs(x:longint):longint;
var i,y:longint;
begin
i:=p[x];
while i<> do
begin
y:=e[i].po;
if not v[y] then
begin
v[y]:=true;
if (cy[y]=-) or (dfs(cy[y])=) then
begin
cy[y]:=x;
cx[x]:=y;
exit();
end;
end;
i:=e[i].next;
end;
exit();
end; function match:longint;
var i:longint;
begin
fillchar(cy,sizeof(cy),);
fillchar(cx,sizeof(cx),);
match:=t1+t2;
if (t2=) or (t1=) then exit;
for i:= to t1 do
if cx[i]=- then
begin
fillchar(v,sizeof(v),false);
match:=match-dfs(i);
end;
end; function cross(i,j,k:longint):longint;
begin
exit((x[i]-x[k])*(y[j]-y[k])-(x[j]-x[k])*(y[i]-y[k]));
end; function dis(i,j:longint):longint;
begin
exit(sqr(x[i]-x[j])+sqr(y[i]-y[j]));
end; begin
readln(n,d);
d:=d*d;
for i:= to n do
readln(x[i],y[i]);
ans:=;
c[]:=;
for i:= to n- do
for j:=i+ to n do
if dis(i,j)<=d then
begin
ll:=dis(i,j);
t1:=;
t2:=;
fillchar(p,sizeof(p),);
fillchar(q,sizeof(q),);
len:=;
for k:= to n do
if (dis(i,k)<=ll) and (dis(j,k)<=ll) then
begin
if cross(k,j,i)>= then
begin
inc(t1);
q1[t1]:=k;
end
else begin
inc(t2);
q2[t2]:=k;
end;
end;
for k:= to t1 do
for l:= to t2 do
if dis(q1[k],q2[l])>d then add(k,l); l:=match;
if l>ans then
begin
ans:=;
fillchar(v1,sizeof(v1),false);
fillchar(v2,sizeof(v2),false);
for k:= to t1 do
if cx[k]=- then
begin
inc(ans);
c[ans]:=q1[k];
z:=p[k];
while z<> do
begin
v2[e[z].po]:=true;
z:=e[z].next;
end;
end; for k:= to t2 do
if cy[k]=- then
begin
inc(ans);
c[ans]:=q2[k];
z:=q[k];
while z<> do
begin
v1[ee[z].po]:=true;
z:=ee[z].next;
end;
end; for k:= to t1 do
if cx[k]<>- then
begin
inc(ans);
if not v1[k] then c[ans]:=q1[k]
else c[ans]:=q2[cx[k]];
end;
end;
end;
writeln(ans);
for i:= to ans do
write(c[i],' ');
writeln;
end.

bzoj4080的更多相关文章

  1. 【随机化】bzoj4080 [Wf2014]Sensor Network

    #include<cstdio> #include<algorithm> #include<cmath> using namespace std; typedef ...

  2. 【BZOJ4080】【WF2014】Sensor Network [随机化]

    Sensor Network Time Limit: 2 Sec  Memory Limit: 128 MB[Submit][Status][Discuss] Description 魔法炮来到了帝都 ...

  3. 【bzoj4080】[Wf2014]Sensor Network 随机化

    题目描述 魔法炮来到了帝都,除了吃特色菜之外,还准备去尝一尝著名的北京烤鸭.帝都一共有n(1<=1<=100)个烤鸭店,可以看成是二维平面内的点.不过由于魔法炮在吃烤鸭之前没有带钱,所以吃 ...

随机推荐

  1. WPF解析Word为图片

    偶遇需要解析Word为单张图,此做 http://git.oschina.net/jiailiuyan/OfficeDecoder using System; using System.Collect ...

  2. Ui设计哪里有好的素材

    刚看到花瓣网,的确不错,以后得多逛逛了.(不喷广告,只留作笔记)

  3. 文字沟通工具使用SignalR,跨域例子源代码

    其他网站已经有很多关于SignalR的介绍了.这里不多介绍. 安装:Install-Package Microsoft.AspNet.SignalR -Version 1.1.4 参考自:http:/ ...

  4. ios 多任务学习笔记

    一.检测多任务是否支持: - (BOOL) isMultitaskingSupported{ BOOL result = NO; if ([[UIDevice currentDevice] respo ...

  5. EXTJS学习笔记

    由于公司统一决策,决定使用EXTJS作为前台的开发UI框架,所以从今天开始学习EXTJS,对学习过程中遇到的问题做记录. 首先第一个问题:为什么要用ExtJS作为公司开发基础架构?   推荐一篇文章: ...

  6. VS2010字体设置+推荐字体

    字体的设置在工具->选项->环境->字体和颜色. 相信大家在用VS2010的时候都会觉得默认的字体不是很好看,尤其是看的时间长了以后,更是累眼睛,这里推荐一个字体,个人感觉像是加粗加 ...

  7. uva 10994

    一开始想法太简单  错了好多遍 #include <cstdio> #include <cstdlib> #include <cmath> #include < ...

  8. Unity3D WP8发布解决方案名 DefaultPackageName 修改问题

    原地址:http://blog.csdn.net/w337198302/article/details/16960661 在对Unity3D游戏进行Windows phone 8 发布的时候,首先是需 ...

  9. VisualSvn Server介绍

    1 .VisualSvn Server  VisualSvn Server是免费的,而VisualSvn是收费的.VisualSvn是Svn的客户端,和Visual Studio集成在一起,但是不免费 ...

  10. 2010年山东省第一届ACM大学生程序设计竞赛 Balloons (BFS)

    题意 : 找联通块的个数,Saya定义两个相连是 |xa-xb| + |ya-yb| ≤ 1 ,但是Kudo定义的相连是 |xa-xb|≤1 并且 |ya-yb|≤1.输出按照两种方式数的联通块的各数 ...