题面

50pts

首先当然是二分答案\(mid\),

对于每一个点,以它为圆心的圆,交上攻击轨道;

那么这个点到攻击轨迹的可达范围就是一段圆弧。

怎么求这段圆弧呢?

我们知道圆弧可以用其两端点对于圆心的弧度(角度)来表示。

已知三角形的三边\(R1,R2,Dist\),利用余弦定理,我们可以求出角\(θ\)。

然后我们也可以轻易得知\(C\)对于圆心\(A\)的弧度\(α\)。

那么圆弧两端点的弧度分别为\(α±θ\)。

假设一个圆弧覆盖正多边形的顶点,那么称圆弧和这个顶点可以匹配。

现在就是问:有没有一个攻击轨道内接正多边形偏转角度,

使得圆弧和正多边形的顶点存在完美匹配。

假设现在存在一种可行的方案,

如果这个正多边形没有一个顶点与某一个圆弧的一个端点重合的话,

那么我可以抖动这个正多边形直到其的一个顶点与某一个圆弧的一个端点重合为止,

这个我是想说,如果存在可行的方案,那么一定存在一种可行的方案,使得该正多边形有一个顶点与某个圆弧的一个端点重合。

由于圆弧的端点的个数有\(O(n)\),

所以枚举重合的那个端点,然后就可以通过这个端点确定正多边形的其他顶点。

然后就可以建出二分图,问是否有完美匹配了,完美匹配可以用\(Hungary\)或\(Maxflow\),\(O(n^3)\)解决。

复杂度是\(O(n^4 log n)\)的,可以过50分。

100pts

首先我们注意到,

虽然每个圆弧的端点对应一种偏转角度

这些偏转角度似乎遍布\(360°\),但其实这些偏转角度是可以\(mod 2*π/n\)的,因为正多边形的顶点互相等价。

所以说我们实际的最大偏转角度是不超过\(2*π/n\)。

先把所有偏转角度\(mod 2*π/n\)后,从小到大排序。

对于一个圆弧而言,由于最大的偏转角度不超过相邻两顶点间距

所以在偏转过程中,不能变得能匹配的顶点不超过\(1\)个,能变得不能匹配的顶点不超过\(1\)个。

也就是说,从\(0\)到最大偏转角度对应的残留网络的变化容量的边不超过\(2*n\)条。

所以我们利用网络流的退流,来只对残留网络进行局部修改。

具体而言:

1.删边

比如说一条\((u,v)\)的边容量变为\(0\),那么假如这条边有\(x\)的流量,那么肯定要把最大流减掉\(x\)。

然后要把这条边的双向容量归零。

从\(u\)到\(1\)找出总流量为\(x\)的增广路,把它们都流量翻转;这个可以直接用最大流来实现,限制住初始进入的流量不超过\(x\)即可;

同理,从\(T\)到\(v\)跑一遍最大流。

最后还要从\(1\)到\(T\)检查有没有别的增广路。

2.加边

这个直接加就好了,最后也要检查从\(1\)到\(T\)有没有新增的增广路。

一次退流/加边操作就是\(O(m)\)的,在本题中\(m=n^2(n?)\)

由于每个点只有一次会退流/加边,所以总共是\(O(n^3)\)。

加上预处理的网络流是\(O(n^3)\)的。

加上二分,总共是\(O(n^3 log n)\)。

代码

#include<bits/stdc++.h>
#define ll long long
#define db double
#define fo(i,x,y) for(int i=x;i<=y;i++)
#define fd(i,x,y) for(int i=x;i>=y;i--)
using namespace std; const int maxn=407,maxm=maxn*maxn;
const db pi=acos(-1); int n,m,Last=1;
db Step; struct P{
int x,y;
}a[maxn]; void Init(){
scanf("%d%d",&n,&m);
fd(i,n,1) scanf("%d%d",&a[i].x,&a[i].y);
//random_shuffle(a+1,a+n+1);
Step=2*pi/n;
} int N;
db q[maxn*2],L[maxn],R[maxn];
void newQuery(db x){
while (x>=-pi+Step) x-=Step;
while (x<-pi) x+=Step;
q[++N]=x;
}
int fi[maxn],ne[maxm],la[maxm],va[maxm],tot,T;
void add(int a,int b,int c){
tot++;
ne[tot]=fi[a];
la[tot]=b;
va[tot]=c;
fi[a]=tot;
}
void Add(int a,int b,int c){
add(a,b,c);
add(b,a,0);
}
int bz[maxn],id;
int sap(int v,int flow){
bz[v]=id;
if (v==T) return flow;
for(int k=fi[v];k;k=ne[k])
if (bz[la[k]]<id && va[k]){
int i=sap(la[k],va[k]&flow);
if (i){
va[k]-=i;
va[k^1]+=i;
return i;
}
}
return 0;
}
int go_flow(int St,int Next,int Lim=0x7fffffff){
int res=0;
T=Next;
while (1){
id++;
int temp=sap(St,1);
if (temp){
res+=temp;
if (res==Lim) break;
}else break;
}
return res;
}
int Road[maxn][maxn],res;
int cnt[maxn];
bool check(db ang){
fo(j,1,n){
if (cnt[j]>=2) continue;
db y=ang;
fo(i,1,n){
if (L[j]>R[j] && (L[j]<=y && y<=R[j]+pi*2 || L[j]-pi*2<=y && y<=R[j]) || L[j]<=R[j] && L[j]<=y && y<=R[j]){
if (!Road[j][i]){
Add(1+j,1+n+i,1);
Road[j][i]=tot;
//cnt[j]++;
/*int z=go_flow(1+n+i,n+n+2)&go_flow(1,1+j);
res+=z;
va[tot]^=z;
va[tot^1]^=z;*/
if (cnt[j]){
cnt[j]++;
break;
}
}
}else{
if (Road[j][i]){
int k=Road[j][i];
if (va[k]){
va[k]=va[k^1]=0;
go_flow(n+n+2,1+n+i,1);
go_flow(1+j,1,1);
res--;
}
va[k^1]=va[k]=0;
Road[j][i]=0;
cnt[j]++;
}
}
y+=Step;
}
}
res+=go_flow(1,n+n+2);
return res>=n;
}
bool judge(db rad){
tot=1;
res=0;
memset(fi,0,sizeof fi);
memset(cnt,0,sizeof cnt);
fo(i,1,n) Add(1,1+i,1);
fo(i,1,n) Add(1+n+i,n+n+2,1);
memset(Road,0,sizeof Road); N=0;
fo(i,1,n){
db Dist=sqrt(a[i].x*a[i].x+a[i].y*a[i].y);
if (rad+m<Dist || Dist<m-rad || a[i].x==0 && a[i].y==0 && rad<m) return false;
if (rad<m+Dist){
db Alpha=atan2(a[i].y,a[i].x),Beta=acos((m*m+Dist*Dist-rad*rad)/(2.0*m*Dist));
L[i]=Alpha-Beta;
R[i]=Alpha+Beta;
if (L[i]<-pi) L[i]+=2*pi;
if (R[i]>=pi) R[i]-=2*pi;
}else{
L[i]=-pi;
R[i]=pi;
}
newQuery(L[i]);
newQuery(R[i]);
}
sort(q+1,q+N+1);
if (n<=10){
fo(i,1,N) if ((i==Last || q[i]-q[i-1]>1e-7) && check(q[i])){
Last=i;
//printf("%d\n",i);
return true;
}
return false;
}
/*fd(i,N,1) if ((i==N || q[i+1]-q[i]>1e-7) && check(q[i])){
printf("%d\n",i);
return true;
}*/
fo(i,Last,N) if ((i==Last || q[i]-q[i-1]>1e-7) && check(q[i])){
Last=i;
//printf("%d\n",i);
return true;
}
fo(i,1,Last-1) if ((i==1 || q[i]-q[i-1]>1e-7) && check(q[i]+Step)){
//Last=i;
//printf("%d\n",i);
return true;
}
return false;
}
db Ans;
void Solve(){
db l=0,r=200;
while (r-l>1e-7){
db mid=(l+r)/2;
if (judge(mid)) r=mid;
else l=mid;
}
Ans=l;
} void Print(){
printf("%.8lf\n",Ans);
} int main(){
freopen("3.in","r",stdin);
freopen("3.out","w",stdout);
Init();
Solve();
Print();
return 0;
}

【JSOI2018】绝地反击的更多相关文章

  1. 【BZOJ5316】[JSOI2018]绝地反击(网络流,计算几何,二分)

    [BZOJ5316][JSOI2018]绝地反击(网络流,计算几何,二分) 题面 BZOJ 洛谷 题解 很明显需要二分一个答案. 那么每个点可以确定的范围就是以当前点为圆心,二分出来的答案为半径画一个 ...

  2. BZOJ5316 : [Jsoi2018]绝地反击

    若$R=0$,那么显然答案为离原点最远的点到原点的距离. 否则若所有点都在原点,那么显然答案为$R$. 否则考虑二分答案$mid$,检查$mid$是否可行. 那么每个点根据对应圆交,可以覆盖圆上的一部 ...

  3. LGP4518[JSOI2018]绝地反击

    题解: 只要确定了每艘飞船的就位位置,就可以用二分+网络流求得答案: 定义偏转角度$a$为离$x$正半轴逆时针最近的边的弧度,$a \in [0,\frac{2\pi}{n})$ 二分一个值,对于一个 ...

  4. 洛谷P4518 [JSOI2018]绝地反击(计算几何+二分图+退流)

    题面 传送门 题解 调了咱一个上午-- 首先考虑二分答案,那么每个点能够到达的范围是一个圆,这个圆与目标圆的交就是可行的区间,这个区间可以用极角来表示 首先,如果我们知道这个正\(n\)边形的转角,也 ...

  5. yyb省选前的一些计划

    突然意识到有一些题目的计划,才可以减少大量查水表或者找题目的时间. 所以我决定这样子处理. 按照这个链接慢慢做. 当然不可能只做省选题了. 需要适时候夹杂一些其他的题目. 比如\(agc/arc/cf ...

  6. LOJ 2548 「JSOI2018」绝地反击 ——二分图匹配+网络流手动退流

    题目:https://loj.ac/problem/2548 如果知道正多边形的顶点,就是二分答案.二分图匹配.于是写了个暴力枚举多边形顶点的,还很愚蠢地把第一个顶点枚举到 2*pi ,其实只要 \( ...

  7. 【LOJ】#2548. 「JSOI2018」绝地反击

    题解 卡常卡不动,我自闭了,特判交上去过了 事实上90pts= = 我们考虑二分长度,每个点能覆盖圆的是一段圆弧 然后问能不能匹配出一个正多边形来 考虑抖动多边形,多边形的一个端点一定和圆弧重合 如果 ...

  8. JSOI2018简要题解

    来自FallDream的博客,未经允许,请勿转载,谢谢. 有幸拜读到贵省的题目,题的质量还不错,而且相比zjoi可做多了,简单发一下题解吧. 还有就是,怎么markdown在博客园上的代码这么丑啊 「 ...

  9. JSOI2018 简要题解

    潜入行动 复杂度分析题. 定义状态fi,j,0/1,0/1f_{i,j,0/1,0/1}fi,j,0/1,0/1​表示以iii为根子树放jjj个机器iii这个放不放,iii这个是否已放来进行dpdpd ...

随机推荐

  1. react 高阶组件之小学版

    高阶组件  多么高大上的概念,一般用来实现组件逻辑的抽象和复用,在很多三方库(redux)中都被使用到,但是开发普通有任务项目时,如果能合理使用高阶组件,也会显著的提高代码质量. 我们今天就用最简单的 ...

  2. 资源-Java:Java资源列表

    ylbtech-资源-Java:Java资源列表 1. 开发软件返回顶部 1.Eclipse https://www.eclipse.org/ 2.IntelliJ IDEA https://www. ...

  3. iOS开发线程之NSThread

    1.初始化 - (instancetype)init API_AVAILABLE(macos(10.5), ios(2.0), watchos(2.0), tvos(9.0)) NS_DESIGNAT ...

  4. QT 获取文件的路径、打开文件的弹出对话框

    1.打开获取文件夹路径的对话框 QString filePath = QFileDialog::getExistingDirectory(this, "请选择文件保存路径...", ...

  5. springboot启动器:spring-boot-starter

    今天想要导入thymeleaf的依赖,但是又不想从其他博复制粘贴,于是去spring官方文档找一找 在idea新建的springbootweb项目中,有一个HELP.md文件,里面包含spring w ...

  6. <爬虫>黑板爬虫闯关01

    import requests from lxml import etree import time ''' 黑板爬虫闯关 网址:http://www.heibanke.com/lesson/craw ...

  7. 基于jdk8的格式化时间方法

    背景 jdk8之前,java使用Date表示时间,在做时间的格式化时,通常使用SimpleDateFormat,但是SimpleDateFormat是非线程安全的,在写代码时通常要将之定义为局部变量或 ...

  8. Leetcode148. Sort List排序链表

    在 O(n log n) 时间复杂度和常数级空间复杂度下,对链表进行排序. 示例 1: 输入: 4->2->1->3 输出: 1->2->3->4 示例 2: 输入 ...

  9. Android开发 Html工具类详解

    前言 在一些需求富文本显示或者编辑的开发情况下,数据都是用html的格式来保存文本信息的.而google是有提供解析html的工具类那就是Html.有了Html可以让TextView也支持富文本(其实 ...

  10. 8年前诞生于淘宝,细数阿里云RPA 的前世今生!

    9月10日,踏入55岁的马云正式卸任阿里巴巴董事局主席一职,由阿里巴巴集团CEO张勇接任.公寓创业.西湖论剑.美国敲钟,从成立到登顶中国最值钱的公司,阿里巴巴只用了20年. 阿里云RPA,2011年诞 ...