【BZOJ1967】[AHOI2005]穿越磁场(最短路)
【BZOJ1967】[AHOI2005]穿越磁场(最短路)
题面
题解
一个显然的思路是这样的,我们的正方形的边长把整个平面割成了若干块,显然每个联通块都可以看着做一个点,那么接下来只需要把所有的有相邻边的联通块全部连上一条长度为\(1\)的边表示如果要从这个联通块到达另外一个联通块,需要穿过这一条边,那么最终的答案就变成了起点和终点所在联通块的最短路。因为找联通块不好搞,其实离散之后再当成网格图做就好了。
#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<cmath>
#include<algorithm>
#include<queue>
using namespace std;
#define ll long long
#define MAX 250
inline int read()
{
int x=0;bool t=false;char ch=getchar();
while((ch<'0'||ch>'9')&&ch!='-')ch=getchar();
if(ch=='-')t=true,ch=getchar();
while(ch<='9'&&ch>='0')x=x*10+ch-48,ch=getchar();
return t?-x:x;
}
struct Matrix{int x1,x2,y1,y2;}p[MAX<<1];
int n,Sx[MAX<<3],Sy[MAX<<3],tx,ty;
int Bx,By,Ex,Ey;
bool visx[MAX<<1][MAX<<1],visy[MAX<<1][MAX<<1];
struct Line{int v,next,w;}e[MAX*MAX*10];
int h[MAX*MAX<<2],cnt=1;
inline void Add(int u,int v,int w){e[cnt]=(Line){v,h[u],w};h[u]=cnt++;}
int d[4][2]={1,0,-1,0,0,1,0,-1};
int id(int x,int y){return x*ty+y-ty;}
int dis[MAX*MAX<<1];bool vis[MAX*MAX<<1];
void SPFA()
{
queue<int> Q;Q.push(id(Bx,By));
memset(dis,63,sizeof(dis));dis[id(Bx,By)]=0;
while(!Q.empty())
{
int u=Q.front();Q.pop();
for(int i=h[u];i;i=e[i].next)
{
int v=e[i].v;
if(dis[v]>dis[u]+e[i].w)
{
dis[v]=dis[u]+e[i].w;
if(!vis[v])vis[v]=true,Q.push(v);
}
}
vis[u]=false;
}
}
int main()
{
n=read();
for(int i=1;i<=n;++i)
{
int x=read(),y=read(),c=read();
p[i]=(Matrix){x,x+c,y+c,y};
}
Bx=read();By=read();Ex=read();Ey=read();
Sx[++tx]=Bx;Sx[++tx]=Ex;Sy[++ty]=By;Sy[++ty]=Ey;
for(int i=1;i<=n;++i)Sx[++tx]=p[i].x1,Sx[++tx]=p[i].x2;
for(int i=1;i<=n;++i)Sy[++ty]=p[i].y1,Sy[++ty]=p[i].y2;
sort(&Sx[1],&Sx[tx+1]);sort(&Sy[1],&Sy[ty+1]);
tx=unique(&Sx[1],&Sx[tx+1])-Sx-1;ty=unique(&Sy[1],&Sy[ty+1])-Sy-1;
for(int i=2;i<=tx;++i)if(Sx[i]-Sx[i-1]>1)Sx[++tx]=Sx[i]-1;
for(int i=2;i<=ty;++i)if(Sy[i]-Sy[i-1]>1)Sy[++ty]=Sy[i]-1;
Sx[++tx]=-10;Sy[++ty]=-10;Sx[++tx]=1e9;Sy[++ty]=1e9;
sort(&Sx[1],&Sx[tx+1]);sort(&Sy[1],&Sy[ty+1]);
tx=unique(&Sx[1],&Sx[tx+1])-Sx-1;ty=unique(&Sy[1],&Sy[ty+1])-Sy-1;
for(int i=1;i<=n;++i)
{
p[i].x1=lower_bound(&Sx[1],&Sx[tx+1],p[i].x1)-Sx;
p[i].x2=lower_bound(&Sx[1],&Sx[tx+1],p[i].x2)-Sx;
p[i].y1=lower_bound(&Sy[1],&Sy[ty+1],p[i].y1)-Sy;
p[i].y2=lower_bound(&Sy[1],&Sy[ty+1],p[i].y2)-Sy;
for(int j=p[i].x1;j<=p[i].x2;++j)visx[j][p[i].y1-1]=visx[j][p[i].y2-1]=true;
for(int j=p[i].y2;j<=p[i].y1;++j)visy[p[i].x1][j]=visy[p[i].x2][j]=true;
}
Bx=lower_bound(&Sx[1],&Sx[tx+1],Bx)-Sx;By=lower_bound(&Sy[1],&Sy[ty+1],By)-Sy;
Ex=lower_bound(&Sx[1],&Sx[tx+1],Ex)-Sx;Ey=lower_bound(&Sy[1],&Sy[ty+1],Ey)-Sy;
for(int i=1;i<=tx;++i)
for(int j=1;j<=ty;++j)
for(int k=0;k<4;++k)
{
int x=i+d[k][0],y=j+d[k][1];
if(x<1||y<1||x>tx||y>ty)continue;
if(k==0)Add(id(i,j),id(x,y),visy[x][y]);
if(k==1)Add(id(i,j),id(x,y),visy[i][j]);
if(k==2)Add(id(i,j),id(x,y),visx[i][j]);
if(k==3)Add(id(i,j),id(x,y),visx[x][y]);
}
SPFA();
printf("%d\n",dis[id(Ex,Ey)]);
return 0;
}
【BZOJ1967】[AHOI2005]穿越磁场(最短路)的更多相关文章
- bzoj1001--最大流转最短路
http://www.lydsy.com/JudgeOnline/problem.php?id=1001 思路:这应该算是经典的最大流求最小割吧.不过题目中n,m<=1000,用最大流会TLE, ...
- 【USACO 3.2】Sweet Butter(最短路)
题意 一个联通图里给定若干个点,求他们到某点距离之和的最小值. 题解 枚举到的某点,然后优先队列优化的dijkstra求最短路,把给定的点到其的最短路加起来,更新最小值.复杂度是\(O(NElogE) ...
- Sicily 1031: Campus (最短路)
这是一道典型的最短路问题,直接用Dijkstra算法便可求解,主要是需要考虑输入的点是不是在已给出的地图中,具体看代码 #include<bits/stdc++.h> #define MA ...
- 最短路(Floyd)
关于最短的先记下了 Floyd算法: 1.比较精简准确的关于Floyd思想的表达:从任意节点A到任意节点B的最短路径不外乎2种可能,1是直接从A到B,2是从A经过若干个节点X到B.所以,我们假设maz ...
- bzoj1266最短路+最小割
本来写了spfa wa了 看到网上有人写Floyd过了 表示不开心 ̄へ ̄ 改成Floyd试试... 还是wa ヾ(。`Д´。)原来是建图错了(样例怎么过的) 结果T了 于是把Floyd改回spfa 还 ...
- HDU2433 BFS最短路
Travel Time Limit: 10000/2000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others)Total Sub ...
- 最短路(代码来源于kuangbin和百度)
最短路 最短路有多种算法,常见的有一下几种:Dijstra.Floyd.Bellman-Ford,其中Dijstra和Bellman-Ford还有优化:Dijstra可以用优先队列(或者堆)优化,Be ...
- Javascript优化细节:短路表达式
什么是短路表达式? 短路表达式:作为"&&"和"||"操作符的操作数表达式,这些表达式在进行求值时,只要最终的结果已经可以确定是真或假,求值过程 ...
- Python中三目计算符的正确用法及短路逻辑
今天在看别人代码时看到这样一种写法, 感觉是个挺容易踩到的坑, 搞清楚后写出来备忘. 短路逻辑 Python中进行逻辑运算的时候, 默认采用的是一种叫做短路逻辑的运算规则. 名字是很形象的, 下面直接 ...
随机推荐
- 软件设计、DDD概念及落地时的一些零碎思考和记录
DDD理解 DDD体现的是对现实的充分尊重. 1.尊重业务现实,领域专家.领域语言等概念 2.尊重团队现实 3.尊重变化 Application 对某一业务线的整体掌控,流程组装,进度管理,存储时机掌 ...
- scala-数组操作
package com.bigdata import scala.collection.mutable.ArrayBuffer object ArrayO { def main(args: Array ...
- c# speech 文本转语言
using System; using System.Collections.Generic; using System.Linq; using System.Text; using System.W ...
- go语言之行--接口(interface)、反射(reflect)详解
一.interface简介 interface(接口)是golang最重要的特性之一,Interface类型可以定义一组方法,但是这些不需要实现.并且interface不能包含任何变量. 简单的说: ...
- 2015520吴思其 基于《Arm试验箱的国密算法应用》课程设计个人报告
20155200吴思其 基于<Arm试验箱的国密算法应用>课程设计个人报告 课程设计中承担的任务 完成试验箱测试功能4,5,6以及SM3加密实验的实现 测试四 GPIO0按键中断实验 实验 ...
- 【转载】C++文件读写详解(ofstream,ifstream,fstream)
原文:http://blog.csdn.net/kingstar158/article/details/6859379 在看C++编程思想中,每个练习基本都是使用ofstream,ifstream,f ...
- 牛客OI周赛4-提高组-C-战争[并查集]
题意 一个长度为 \(n\) 的序列,每个权值互不相同,给出形如 \(l,r,p\) 的信息表示 \([l,r]\) 区间中最小的数是 \(p\) ,问第几个信息开始出现矛盾. \(n\leq 5 \ ...
- SSIS 更新变量
在Package中声明一个variable,在package运行的过程中,SSIS如何update Variable? 第一种方法:使用 Script Task 来更新Variable的值 1,创建一 ...
- stl源码剖析 详细学习笔记deque(3)
protected: typedef simple_alloc<value_type,Alloc> data_allocator; //用来配置元素的alloc typedef simpl ...
- 用C++实现一个Brainfuck解释器
Brainfuck是一种极小化的计算机语言,只含有8种运算符,由于fuck在英语中是脏话,这种语言有时被称为brainfck或brainf**,甚至被简称为BF.正如它的名字所暗示,brainfuck ...