[bzoj3630][JLOI2014]镜面通道_计算几何_网络流_最小割
镜面通道 bzoj-3630 JLOI-2014
题目大意:题目链接。
注释:略。
想法:
我们发现,只要上下界没有被完全封死,我们就一定有一条合法的光路。
所以只需要将上界和下界拆开即可。
拆点,把每个点分为入点和出点,入点向出点连1的边。
元件之间如果连通就连$inf$。
和上界连通就对源点连$inf$,和下界连通就和汇点连$inf$。
最后求最小割即为答案。
Code:
// luogu-judger-enable-o2
#include<iostream>
#include<cstdio>
#include<cstring>
#define N 303
#define M 370000
#define inf 0x7fffffff
#define ll long long
using namespace std;
int to[M],hd[M],lk[N<<1],len[M],cnt=1,T;
void add(int u,int v,int w)
{
to[++cnt]=v,hd[cnt]=lk[u],len[cnt]=w,lk[u]=cnt;
to[++cnt]=u,hd[cnt]=lk[v],lk[v]=cnt;
}
struct pt
{
ll x,y;
void read()
{scanf("%lld%lld",&x,&y);}
}tp1,tp2;
inline ll Dis(pt a,pt b)
{return (a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y);}
struct cir
{pt o;int r;}Cir[N];
struct rec
{pt x1,x2;}Rec[N],ss,tt;
short n,loc[N];
int X,Y,ans;
ll d,xx,yy;
bool judge(cir a,rec b)
{
d=a.r,xx=a.o.x,yy=a.o.y;
if(xx>=b.x1.x-d&&xx<=b.x2.x+d&&yy>=b.x1.y&&yy<=b.x2.y
||xx>=b.x1.x&&xx<=b.x2.x&&yy>=b.x1.y-d&&yy<=b.x2.y+d)return 1;
tp1.x=b.x1.x,tp2.y=b.x2.y,tp2.x=b.x2.x,tp2.y=b.x1.y;d*=d;
return Dis(a.o,tp1)<=d||Dis(a.o,tp2)<=d||Dis(a.o,b.x1)<=d||Dis(a.o,b.x2)<=d;
}
bool check(int a,int b)
{
if(loc[a]>loc[b])swap(a,b);
if(loc[b]<2)return((Cir[a].r+Cir[b].r)*(Cir[a].r+Cir[b].r)>=Dis(Cir[a].o,Cir[b].o));
if(loc[a]>1)return (Rec[a].x1.x<=Rec[b].x1.x&&Rec[b].x1.x<=Rec[a].x2.x||
Rec[a].x1.x<=Rec[b].x2.x&&Rec[b].x2.x<=Rec[a].x2.x)&&
(Rec[a].x1.y<=Rec[b].x1.y&&Rec[b].x1.y<=Rec[a].x2.y||
Rec[a].x1.y<=Rec[b].x2.y&&Rec[b].x2.y<=Rec[a].x2.y);
return judge(Cir[a],Rec[b]);
}
int k,q[N<<1],h,t,dis[N<<1],cur[N<<1];
bool bfs()
{
for(int i=0;i<=T;i++)
dis[i]=0,cur[i]=lk[i];
h=0,t=dis[0]=1;
while(h<t)
{
k=q[h++];
for(int i=lk[k];i;i=hd[i])
if(len[i]&&!dis[to[i]])
dis[q[t++]=to[i]]=dis[k]+1;
}
return dis[T];
}
int dfs(int x,int f)
{
if(x==T||!f)return f;
int r=0,cst;
for(int &i=cur[x];i;i=hd[i])
if(dis[to[i]]==dis[x]+1)
{
cst=dfs(to[i],f<len[i]?f:len[i]);
f-=cst,len[i]-=cst;
r+=cst,len[i^1]+=cst;
if(!f)break;
}
if(!r)dis[x]=0;
return r;
}
int main()
{
scanf("%d%d",&X,&Y);
ss={{0,0},{X,0}};tt={{0,Y},{X,Y}};
scanf("%d",&n);T=n<<1|1;
for(int i=1;i<=n;i++)
{
add(i,i+n,1);
scanf("%d",loc+i);
if(loc[i]<2)
{
Cir[i].o.read(),scanf("%lld",&Cir[i].r);
if(judge(Cir[i],ss))add(0,i,inf);
if(judge(Cir[i],tt))add(i+n,T,inf);
}
else
{
Rec[i].x1.read(),Rec[i].x2.read();
if(Rec[i].x1.y<=0&&Rec[i].x1.x<=X&&Rec[i].x2.x>=0)add(0,i,inf);
if(Rec[i].x2.y>=Y&&Rec[i].x1.x<=X&&Rec[i].x2.x>=0)add(i+n,T,inf);
}
for(int j=1;j<i;j++)
if(check(j,i))
add(j+n,i,inf),add(i+n,j,inf);
}
while(bfs())ans+=dfs(0,inf);
printf("%d",ans);
}
小结:好题好题。就是特判有点恶心。
[bzoj3630][JLOI2014]镜面通道_计算几何_网络流_最小割的更多相关文章
- 【bzoj3630】[JLOI2014]镜面通道 对偶图+计算几何+网络流最小割
题目描述 在一个二维平面上,有一个镜面通道,由镜面AC,BD组成,AC,BD长度相等,且都平行于x轴,B位于(0,0).通道中有n个外表面为镜面的光学元件,光学元件α为圆形,光学元件β为矩形(这些元件 ...
- BZOJ3630 : [JLOI2014]镜面通道
从左边不能到达右边当且仅当存在一条与上下底边相连的分割线将它们分开 设下底边为S,上底边为T,每个元件作为点,有公共部分的两个点互相连边 最后拆点求最小割 #include<cstdio> ...
- 【BZOJ3630】[JLOI2014]镜面通道 几何+最小割
[BZOJ3630][JLOI2014]镜面通道 Description 在一个二维平面上,有一个镜面通道,由镜面AC,BD组成,AC,BD长度相等,且都平行于x轴,B位于(0,0).通道中有n个外表 ...
- 【BZOJ2229】[ZJOI2011]最小割(网络流,最小割树)
[BZOJ2229][ZJOI2011]最小割(网络流,最小割树) 题面 BZOJ 洛谷 题解 戳这里 那么实现过程就是任选两点跑最小割更新答案,然后把点集划分为和\(S\)联通以及与\(T\)联通. ...
- 【CF331E】Biologist(网络流,最小割)
[CF331E]Biologist(网络流,最小割) 题面 洛谷 翻译: 有一个长度为\(n\)的\(01\)串,将第\(i\)个位置变为另外一个数字的代价是\(v_i\). 有\(m\)个要求 每个 ...
- 【BZOJ1391】Order(网络流,最小割)
[BZOJ1391]Order(网络流,最小割) 题面 BZOJ权限题... 良心洛谷 题目描述 有N个工作,M种机器,每种机器你可以租或者买过来. 每个工作包括若干道工序,每道工序需要某种机器来完成 ...
- 【BZOJ3144】切糕(网络流,最小割)
[BZOJ3144]切糕(网络流,最小割) 题面 BZOJ 题解 这样的类型很有趣 先不考虑相邻距离差不能超过\(D\)的限制 我们考虑答案,显然就是在每个位置选一个最小的高度割就行了 化成最小割的模 ...
- [TJOI2013]攻击装置(网络流,最小割)
前言 网络流被hbx吊起来打 Solution 考虑一下这个走法是不是和象棋中马的走法一模一样(废话) 那么显然我每一次移动是走三次,如果将棋盘二分图染色一下,不就是每一次只能走到另一个颜色的吗? 然 ...
- cogs750栅格网络流(最小割)
750. 栅格网络流 ★★☆ 输入文件:flowa.in 输出文件:flowa.out 简单对比时间限制:1 s 内存限制:128 MB [问题描述] Bob 觉得一般图的最大流问题太 ...
随机推荐
- AJPFX总结关于Java中过滤出字母、数字和中文的正则表达式
1.Java中过滤出字母.数字和中文的正则表达式 (1)过滤出字母的正则表达式 [^(A-Za-z)] (2) 过滤出 数字 的正则表达式 [^(0-9)] (3) 过滤出 中文 的正则 ...
- AJPFX简述java语言现状和发展
作为一种最流行的网络编程语言之一,java语言在当今信息化社会中发挥了 重要的作用.Java语言具有面向对象.跨平台.安全性.多线程等特点,这使得java成为许多应用系统的理想开发语言.java应用在 ...
- 5.1点击4个按钮显示相应的div
事件:onclick 属性:display,className 用到for语句,index标记,this当前事件 先清空后附加 <!DOCTYPE html><html>< ...
- SQL——视图、事务、锁、存储过程
https://www.bilibili.com/video/av15496406/?p=57 https://blog.csdn.net/u013630349/article/details/750 ...
- ERROR 1 (HY000): Can't create/write to file '/tmp/#sql_830_0.MYI' (Errcode: 13)
mysql操作时,出现报错. 执行describe 命令时, 临时文件目录没有创建或者无写入权限:于是: cd /var/lib/mysql/ #进入mysql数据目录 mkdir tmp #创建需要 ...
- Android(java)学习笔记189:ContentProvider使用(银行数据库创建和增删改查的案例)
1. Android的四大组件: (1)Activity 用户交互的UI界面 (2)Service 后台运行的服务 (3)BroadcastReceiver 广播接收者 (4)ContentPro ...
- jekyll 将纯文本转化为静态网站和博客 静态网站生成器
jekyll 将纯文本转化为静态网站和博客 静态网站生成器 这个貌似对windows 支持不是很好~ 但是有支持,官方说不建议使用
- Discuz!代码
我如何使用Discuz!代码 Discuz!代码 效果 [b]粗体文字 Abc[/b] 粗体文字 Abc [i]斜体文字 Abc[/i] 斜体文字 Abc [u]下划线文字 Abc[/u] 下划线 ...
- java线程池 多线程 搜索包含关键字的文件路径
package org.jimmy.searchfile20180807.main; public class ThreadMain implements Runnable{ private int ...
- html嵌入pdf && html嵌入多媒体文件,word,flash,pdf,音视频
<object classid="clsid:CA8A9780-280D-11CF-A24D-444553540000" width="1000" hei ...