bzoj 2304: [Apio2011]寻路【spfa】
我是智障*3,读优写错了调了半天没发现= =
虽然是个最短路却有网络流一般的神建图啊。
首先发现在拐角处转弯是最优的,于是先离散化,然后矩形的四个顶点向距离它最近的上下左右点连边,跑spfa即可。
就是难写啊,还要判断无解:st在矩形里;dis[t]=inf
#include<iostream>
#include<cstdio>
#include<queue>
#include<cstring>
#include<algorithm>
#include<map>
#include<cstdlib>
using namespace std;
const int N=2005;
const long long inf=1e15;
int n,m,sx,sy,tx,ty,a[N],b[N],c[N],d[N],hax,hay,x[N],y[N],mk[N][N],tot,s,t,fx[N],fy[N],h[24005],cnt;
long long dis[24005];
bool v[N][N],fl[N],vis[24005];
map<int,int>mx,my;
struct edg
{
int ne,to;
long long va;
}e[200005];
struct qwe
{
int x,y,id;
qwe(int X=0,int Y=0,int ID=0)
{
x=X,y=Y,id=ID;
}
}p[24005];
bool cmpx(const qwe &a,const qwe &b)
{
return a.x<b.x||(a.x==b.x&&a.y<b.y);
}
bool cmpy(const qwe &a,const qwe &b)
{
return a.y<b.y||(a.y==b.y&&a.x<b.x);
}
int read()
{
int r=0,f=1;
char p=getchar();
while(p>'9'||p<'0')
{
if(p=='-')
f=-1;
p=getchar();
}
while(p>='0'&&p<='9')
{
r=r*10+p-48;
p=getchar();
}
return r*f;
}
inline long long jl(qwe a,qwe b)
{
return abs(fx[a.x]-fx[b.x])+abs(fy[a.y]-fy[b.y]);
}
void zou(int x,int y,int dx,int dy)
{
x+=dx,y+=dy;
while(!mk[x][y])
{
if(x<1||x>hax||y<1||y>hay)
return;
x+=dx,y+=dy;
}
v[x][y]=1;
}
inline void add(int u,int v,long long w)
{
cnt++;
e[cnt].ne=h[u];
e[cnt].to=v;
e[cnt].va=w;
h[u]=cnt;
}
inline void ins(int u,int v,long long w)
{//cerr<<u<<" "<<v<<" "<<w<<endl;
add(u,v,w);
add(v,u,w);
}
int main()
{
// int T=read();
int T;
scanf("%d",&T);
while(T--)
{
mx.clear(),my.clear();
memset(h,0,sizeof(h));
memset(v,0,sizeof(v));
memset(mk,0,sizeof(mk));
memset(fl,0,sizeof(fl));
memset(vis,0,sizeof(vis));
cnt=0;tot=0,hax=0,hay=0;
// sx=read(),sy=read(),tx=read(),ty=read(),n=read();
scanf("%d%d%d%d%d",&sx,&sy,&tx,&ty,&n);
for(int i=1;i<=n;i++)
{
// a[i]=read(),b[i]=read(),c[i]=read(),d[i]=read();
scanf("%d%d%d%d",&a[i],&b[i],&c[i],&d[i]);
if(a[i]>c[i])
swap(a[i],c[i]);
if(b[i]>d[i])
swap(b[i],d[i]);
x[i]=a[i],x[i+n]=c[i],y[i]=b[i],y[i+n]=d[i];
}
x[2*n+1]=sx,x[2*n+2]=tx,y[2*n+1]=sy,y[2*n+2]=ty;
sort(x+1,x+1+2*n+2);
sort(y+1,y+1+2*n+2);
for(int i=1;i<=2*n+2;i++)
if(i==1||x[i]!=x[i-1])
mx[x[i]]=++hax,fx[hax]=x[i];
for(int i=1;i<=2*n+2;i++)
if(i==1||y[i]!=y[i-1])
my[y[i]]=++hay,fy[hay]=y[i];//cerr<<y[i]<<" "<<my[y[i]]<<endl;
sx=mx[sx],sy=my[sy],tx=mx[tx],ty=my[ty];
mk[sx][sy]=mk[tx][ty]=-1;
for(int i=1;i<=n;i++)
for(int j=1;j<=n;j++)
if(i!=j&&a[i]>=a[j]&&b[i]>=b[j]&&c[i]<=c[j]&&d[i]<=d[j])
{
fl[i]=1;
break;
}
for(int i=1;i<=n;i++)
if(!fl[i])
{
a[i]=mx[a[i]],b[i]=my[b[i]],c[i]=mx[c[i]],d[i]=my[d[i]];
//cerr<<a[i]<<" "<<b[i]<<" "<<c[i]<<" "<<d[i]<<endl;
for(int j=a[i];j<=c[i];j++)
mk[j][b[i]]=1,mk[j][d[i]]=2;
for(int j=b[i];j<=d[i];j++)
mk[a[i]][j]=1,mk[c[i]][j]=2;
mk[a[i]][b[i]]=mk[a[i]][d[i]]=mk[c[i]][b[i]]=mk[c[i]][d[i]]=-1;
}
for(int i=1;i<=n;++i)
if((sx>a[i]&&sy>b[i]&&sx<c[i]&&sy<d[i])||(tx>a[i]&&ty>b[i]&&tx<c[i]&&ty<d[i]))
{
puts("No Path");
return 0;
}
for(int i=1;i<=n;i++)
if(!fl[i])
{
zou(a[i],b[i],-1,0),zou(a[i],b[i],0,-1);
zou(a[i],d[i],-1,0),zou(a[i],d[i],0,1);
zou(c[i],b[i],1,0),zou(c[i],b[i],0,-1);
zou(c[i],d[i],1,0),zou(c[i],d[i],0,1);
v[a[i]][b[i]]=v[a[i]][d[i]]=v[c[i]][b[i]]=v[c[i]][d[i]]=1;
}
v[sx][sy]=v[tx][ty]=1;
zou(sx,sy,-1,0),zou(sx,sy,1,0),zou(sx,sy,0,-1),zou(sx,sy,0,1);
zou(tx,ty,-1,0),zou(tx,ty,1,0),zou(tx,ty,0,-1),zou(tx,ty,0,1);
for(int i=1;i<=hax;i++)
for(int j=1;j<=hay;j++)
if(v[i][j])
{
tot++;
p[tot]=qwe(i,j,tot);
if(i==sx&&j==sy)
s=tot;
if(i==tx&&j==ty)
t=tot;
}//cerr<<s<<" "<<t<<endl;
sort(p+1,p+1+tot,cmpx);
for(int i=1,j;i<=tot;i=j)
for(j=i+1;j<=tot&&p[j].x==p[i].x;j++)
if(mk[p[j].x][p[j].y]==-1||mk[p[j-1].x][p[j-1].y]==-1||mk[p[j-1].x][p[j-1].y]>=mk[p[j].x][p[j].y])
ins(p[j-1].id,p[j].id,jl(p[j-1],p[j]));
sort(p+1,p+1+tot,cmpy);
for(int i=1,j;i<=tot;i=j)
for(j=i+1;j<=tot&&p[j].y==p[i].y;j++)
if(mk[p[j].x][p[j].y]==-1||mk[p[j-1].x][p[j-1].y]==-1||mk[p[j-1].x][p[j-1].y]>=mk[p[j].x][p[j].y])
ins(p[j-1].id,p[j].id,jl(p[j-1],p[j]));
queue<int>q;
for(int i=1;i<=tot;i++)
dis[i]=inf;
dis[s]=0;
vis[s]=1;
q.push(s);
while(!q.empty())
{
int u=q.front();
q.pop();
vis[u]=0;
for(int i=h[u];i;i=e[i].ne)
if(dis[e[i].to]>dis[u]+e[i].va)
{
dis[e[i].to]=dis[u]+e[i].va;
if(!vis[e[i].to])
{
vis[e[i].to]=1;
q.push(e[i].to);
}
}
}
if(dis[t]==inf)
puts("No Path");
else
printf("%lld\n",dis[t]);
}
return 0;
}
bzoj 2304: [Apio2011]寻路【spfa】的更多相关文章
- bzoj 2304 [Apio2011]寻路 Dij+模拟+恶心建图
[Apio2011]寻路 Time Limit: 20 Sec Memory Limit: 256 MBSubmit: 624 Solved: 193[Submit][Status][Discus ...
- BZOJ 4898 [APIO2017] 商旅 | SPFA判负环 分数规划
BZOJ 4898 [APIO2017] 商旅 | SPFA判负环 分数规划 更清真的题面链接:https://files.cnblogs.com/files/winmt/merchant%28zh_ ...
- bzoj 3875 骑士游戏 - spfa - 动态规划
Description [故事背景] 长期的宅男生活中,JYY又挖掘出了一款RPG游戏.在这个游戏中JYY会 扮演一个英勇的骑士,用他手中的长剑去杀死入侵村庄的怪兽. [问题描述] 在这个游戏中,J ...
- BZOJ 2763: [JLOI2011]飞行路线 spfa dp
题目链接: http://www.lydsy.com/JudgeOnline/problem.php?id=2763 题解: d[x][kk]表示从s到x用了kk次免费机会的最少花费. 代码: #in ...
- BZOJ 2303: [Apio2011]方格染色 题解
题目大意: 有n*m的方格,中间的数要么是1,要么是0,要求任意2*2的方格中的数异或和为1.已知一部分格子中的数,求合法的填数的方案数. 思路: 由题意得:a[i][j]^a[i][j+1]^a[i ...
- bzoj 2303: [Apio2011]方格染色
传送门 Description Sam和他的妹妹Sara有一个包含n × m个方格的表格.她们想要将其的每个方格都染成红色或蓝色.出于个人喜好,他们想要表格中每个2 × 2的方形区域都包含奇数个(1 ...
- BZOJ 2303: [Apio2011]方格染色 [并查集 数学!]
题意: $n*m:n,m \le 10^6$的网格,每个$2 \times 2$的方格必须有1个或3个涂成红色,其余涂成蓝色 有一些方格已经有颜色 求方案数 太神了!!!花我三节课 首先想了一下只有两 ...
- bzoj 2303: [Apio2011]方格染色【并查集】
画图可知,每一行的状态转移到下一行只有两种:奇数列不变,偶数列^1:偶数列不变,奇数列^1 所以同一行相邻的变革染色格子要放到同一个并查集里,表示这个联通块里的列是联动的 最后统计下联通块数(不包括第 ...
- bzoj AC倒序
Search GO 说明:输入题号直接进入相应题目,如需搜索含数字的题目,请在关键词前加单引号 Problem ID Title Source AC Submit Y 1000 A+B Problem ...
随机推荐
- BZOJ 1123 tarjan
题目链接 题意:一张无向图,把第$i$个点关联的所有边去掉,求无向图中有多少个点对不连通. 题解: 如果割的不是割点,那么总答案是$2\times (n-1)$. 如果是割点,要分别考虑每个子树的贡献 ...
- Fibonacci--poj3070(矩阵快速幂)
http://poj.org/problem?id=3070 Description In the Fibonacci integer sequence, F0 = 0, F1 = 1, and Fn ...
- UVA 12035 War Map
让我先把微笑送给出题人 这个题最基础的一个想法:先找出一个度数和为总度数和的1/2的点集,然后判断这个点集和这个点集的补集能否形成二分图.但是就算我们把判断的复杂度看成O(1),这个算法的复杂度也是 ...
- 阿里oss上传图片react组件alioss-react,vue组件alioss-vue (不用我先收藏着,后端看下前端处理方法)
1.介绍 最近开发了一个项目,其中需要一个上传图片到阿里云的 oss 上面,就是上传图片到阿里云的 oss 上面. 因为之前开发过 vue 的阿里云 oss 上传,所以直接复制粘 vue 的组件. 因 ...
- Nginx: 解决connect() to xxxx failed (13: Permission denied) while connecting to upstream的问题
一句话:setsebool httpd_can_network_connect true
- 怎样使用SSH连接OpenStack上的云主机
转载请注明出处.否则将追究法律责任http://blog.csdn.net/xingjiarong/article/details/47021815 在上一篇博客中我介绍了怎样在OpenStack中创 ...
- discuz的php7版本号
php7的安装 wget http://am1.php.net/get/php-7.0.4.tar.gz/from/this/mirror tar zvxf php-7.0.4.tar.gz cd p ...
- dva/dynamic
1.安装: yarn add dva 2.引入: import dynamic from 'dva/dynamic'; * dva路由跳转 * dynamic(app, model, componen ...
- Web安全漏洞及攻击
背景介绍 先说一个在互联网上常见,但是普通人又不太理解的东西--“验证码”. 验证码(CAPTCHA)是“Completely Automated Public Turing test to tell ...
- mina客户端与服务端通信的易错点
使用mina进行项目开发时,如果客户端与服务端不在同一个项目下,需要关注一下两点: 第一.服务端与客户端的编码解码器一致 第二.过程中所用到的实体类的包名需要一致