【JSOI2018】绝地反击
题面
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】绝地反击的更多相关文章
- 【BZOJ5316】[JSOI2018]绝地反击(网络流,计算几何,二分)
[BZOJ5316][JSOI2018]绝地反击(网络流,计算几何,二分) 题面 BZOJ 洛谷 题解 很明显需要二分一个答案. 那么每个点可以确定的范围就是以当前点为圆心,二分出来的答案为半径画一个 ...
- BZOJ5316 : [Jsoi2018]绝地反击
若$R=0$,那么显然答案为离原点最远的点到原点的距离. 否则若所有点都在原点,那么显然答案为$R$. 否则考虑二分答案$mid$,检查$mid$是否可行. 那么每个点根据对应圆交,可以覆盖圆上的一部 ...
- LGP4518[JSOI2018]绝地反击
题解: 只要确定了每艘飞船的就位位置,就可以用二分+网络流求得答案: 定义偏转角度$a$为离$x$正半轴逆时针最近的边的弧度,$a \in [0,\frac{2\pi}{n})$ 二分一个值,对于一个 ...
- 洛谷P4518 [JSOI2018]绝地反击(计算几何+二分图+退流)
题面 传送门 题解 调了咱一个上午-- 首先考虑二分答案,那么每个点能够到达的范围是一个圆,这个圆与目标圆的交就是可行的区间,这个区间可以用极角来表示 首先,如果我们知道这个正\(n\)边形的转角,也 ...
- yyb省选前的一些计划
突然意识到有一些题目的计划,才可以减少大量查水表或者找题目的时间. 所以我决定这样子处理. 按照这个链接慢慢做. 当然不可能只做省选题了. 需要适时候夹杂一些其他的题目. 比如\(agc/arc/cf ...
- LOJ 2548 「JSOI2018」绝地反击 ——二分图匹配+网络流手动退流
题目:https://loj.ac/problem/2548 如果知道正多边形的顶点,就是二分答案.二分图匹配.于是写了个暴力枚举多边形顶点的,还很愚蠢地把第一个顶点枚举到 2*pi ,其实只要 \( ...
- 【LOJ】#2548. 「JSOI2018」绝地反击
题解 卡常卡不动,我自闭了,特判交上去过了 事实上90pts= = 我们考虑二分长度,每个点能覆盖圆的是一段圆弧 然后问能不能匹配出一个正多边形来 考虑抖动多边形,多边形的一个端点一定和圆弧重合 如果 ...
- JSOI2018简要题解
来自FallDream的博客,未经允许,请勿转载,谢谢. 有幸拜读到贵省的题目,题的质量还不错,而且相比zjoi可做多了,简单发一下题解吧. 还有就是,怎么markdown在博客园上的代码这么丑啊 「 ...
- JSOI2018 简要题解
潜入行动 复杂度分析题. 定义状态fi,j,0/1,0/1f_{i,j,0/1,0/1}fi,j,0/1,0/1表示以iii为根子树放jjj个机器iii这个放不放,iii这个是否已放来进行dpdpd ...
随机推荐
- 运用shtml类型文件,实现项目页面的分割。
学过动态网页,如asp.php的人知道如何去引用网站头部.底部文件包含.例如在php中有一个方法:include();这个方法能在当前文档中引入外部文件,从而方便网站的开发和维护,然而html静态文件 ...
- linux下常见的包安装方式
linux下常见的包安装方式 一.总结 一句话总结: rpm包安装 tar.gz源代码包安装 yum方式安装rpm包 bin文件安装 1.yum是什么? 安装所有依赖的软件包 Yum(全称为 Yell ...
- Jenkins 自动部署
一.安装插件[系统管理 → 插件管理 ] 为了通过SSH上传war包,我们需要安装Publish Over SSH 插件. 二.添加SSH 服务器[系统管理→系统设置] 参数说明: Name:ss ...
- 本地项目通过 git 同步到 github
1. github创建仓库并克隆仓库地址 2. 在本地通过git命令:git clone <仓库地址> 生成github仓库文件夹 3. 将本地项目复制到该文件夹 4. 通过git命令:g ...
- 使用virtualenv发布Python程序
客户环境不能上网,开始想把所有依赖包下载下来,进入客户环境进行安装.但为了避免出差,部署工作交给其他同事了,我想还是需要更简单的方式. 实验了一下virtualenv是可以的 1. 创建一个新的环境( ...
- 注解@Qualifier@Primary
在Controller中需要注入service那么我的这个server有两个实现类如何区分开这两个impl呢 根据注入资源的注解不同实现的方式有一点小小的区别 在Controller中使用 @Auto ...
- java、jsp导出excel功能备份
问题踩坑: ajax请求不能下载文件 必须这样: <a href="/media">点击下载Excel</a> 或者 location.href = '/m ...
- Python全栈开发:json与pickle
#!/usr/bin/env python # -*- coding;utf-8 -*- """ 正解(序列化):将Python数据类型转换成json或者pickle格式 ...
- JavaWeb中请求转发和请求重定向的区别
针对于JavaWeb中请求与重定向的一个cheatsheep: 1.转发 1)完成一次转发,用户浏览器发送一次请求 2)转发之后,浏览器URL地址栏不改变(服务器帮忙完成) 3)请求域中数据不丢失 4 ...
- ThinkPHP实现了ActiveRecords模式的ORM模型
ThinkPHP实现了ActiveRecords模式的ORM模型,采用了非标准的ORM模型:表映射到类,记录映射到对象.最大的特点就是使用方便和便于理解(因为采用了对象化),提供了开发的最佳体验,从而 ...