BZOJ2647 : [Neerc2011]Journey
$|x|+|y|=\max(x+y,x-y,-x+y,-x-y)$,设$f[i][j]$表示在$(0,0)$,朝向方向$j$,执行第$i$条指令后的信息:
$cir$:是否陷入循环
$d$:朝向
$x,y$:坐标
$v_0$:$\max(x+y)$
$v_1$:$\max(x-y)$
$v_2$:$\max(-x+y)$
$v_3$:$\max(-x-y)$
然后记忆化搜索,当陷入循环时不再执行后面的指令。
用栈维护所有正在计算的状态,若发现循环,则计算循环的$x$和$y$之和,若为$0$则是循环,否则会走到无限远处。
时间复杂度$O(n^3)$。
#include<cstdio>
#include<string>
#include<iostream>
#include<cstdlib>
using namespace std;
const int N=110,Base=100000000,MAXL=32;
int n,i,j,a[N],q[N*4][2],in[N][4];string b[N][N];bool vis[N][4],cir[N][4];
inline int max(int a,int b){return a>b?a:b;}
struct Num{
int a[MAXL],len,fu;
void init(){len=1,fu=a[1]=0;}
Num(){init();}
bool iszero(){return len==1&&!a[1];}
Num operator+(const Num&b){
Num c;
c.len=max(len,b.len)+2;
int i;
for(i=1;i<=c.len;i++)c.a[i]=0;
if(fu==b.fu){
for(i=1;i<=len;i++)c.a[i]=a[i];
for(i=1;i<=b.len;i++)c.a[i]+=b.a[i];
for(i=1;i<=c.len;i++)if(c.a[i]>=Base)c.a[i+1]++,c.a[i]-=Base;
while(c.len>1&&!c.a[c.len])c.len--;
c.fu=fu;
}else{
bool flag=0;
if(len==b.len){
for(i=len;i;i--)if(a[i]!=b.a[i]){
if(a[i]>b.a[i])flag=1;
break;
}
}else{
if(len>b.len)flag=1;
}
if(flag){
for(i=1;i<=len;i++)c.a[i]=a[i];
for(i=1;i<=b.len;i++)c.a[i]-=b.a[i];
for(i=1;i<=c.len;i++)if(c.a[i]<0)c.a[i+1]--,c.a[i]+=Base;
while(c.len>1&&!c.a[c.len])c.len--;
c.fu=fu;
}else{
for(i=1;i<=b.len;i++)c.a[i]=b.a[i];
for(i=1;i<=len;i++)c.a[i]-=a[i];
for(i=1;i<=c.len;i++)if(c.a[i]<0)c.a[i+1]--,c.a[i]+=Base;
while(c.len>1&&!c.a[c.len])c.len--;
c.fu=b.fu;
}
}
if(c.iszero())c.init();
return c;
}
void operator+=(const Num&b){*this=*this+b;}
Num operator-(const Num&b){
Num c;
c.len=max(len,b.len)+2;
int i;
for(i=1;i<=c.len;i++)c.a[i]=0;
if(fu!=b.fu){
for(i=1;i<=len;i++)c.a[i]=a[i];
for(i=1;i<=b.len;i++)c.a[i]+=b.a[i];
for(i=1;i<=c.len;i++)if(c.a[i]>=Base)c.a[i+1]++,c.a[i]-=Base;
while(c.len>1&&!c.a[c.len])c.len--;
c.fu=fu;
}else{
bool flag=0;
if(len==b.len){
for(i=len;i;i--)if(a[i]!=b.a[i]){
if(a[i]>b.a[i])flag=1;
break;
}
}else{
if(len>b.len)flag=1;
}
if(flag){
for(i=1;i<=len;i++)c.a[i]=a[i];
for(i=1;i<=b.len;i++)c.a[i]-=b.a[i];
for(i=1;i<=c.len;i++)if(c.a[i]<0)c.a[i+1]--,c.a[i]+=Base;
while(c.len>1&&!c.a[c.len])c.len--;
c.fu=fu;
}else{
for(i=1;i<=b.len;i++)c.a[i]=b.a[i];
for(i=1;i<=len;i++)c.a[i]-=a[i];
for(i=1;i<=c.len;i++)if(c.a[i]<0)c.a[i+1]--,c.a[i]+=Base;
while(c.len>1&&!c.a[c.len])c.len--;
c.fu=b.fu^1;
}
}
if(c.iszero())c.init();
return c;
}
void operator-=(const Num&b){*this=*this-b;}
bool operator<(const Num&b){
if(fu!=b.fu)return fu;
if(fu){
if(len!=b.len)return len>b.len;
for(int i=len;i;i--)if(a[i]!=b.a[i])return a[i]>b.a[i];
}else{
if(len!=b.len)return len<b.len;
for(int i=len;i;i--)if(a[i]!=b.a[i])return a[i]<b.a[i];
}
return 0;
}
void write(){
if(fu)putchar('-');
printf("%d",a[len]);
for(int i=len-1;i;i--)printf("%08d",a[i]);
}
void set(int x){
fu=0;
if(x<0)x=-x,fu=1;
len=1;
a[1]=x;
}
}zero,one,x,y;
inline void up(Num&a,const Num&b){if(a<b)a=b;}
struct E{
int d;Num x,y,v0,v1,v2,v3;
void left(){d=(d+1)&3;}
void right(){d=(d+3)&3;}
void init(int _d){
d=_d;
x.set(0);
y.set(0);
v0.set(0);
v1.set(0);
v2.set(0);
v3.set(0);
}
void go(){
if(d==0)x+=one;
if(d==1)y+=one;
if(d==2)x-=one;
if(d==3)y-=one;
up(v0,x+y);
up(v1,x-y);
up(v2,y-x);
up(v3,zero-x-y);
}
void merge(const E&b){
d=b.d;
up(v0,x+y+b.v0);
up(v1,x-y+b.v1);
up(v2,y-x+b.v2);
up(v3,zero-x-y+b.v3);
x+=b.x;
y+=b.y;
}
void write(){
up(v0,v1);
up(v0,v2);
up(v0,v3);
v0.write();
}
}f[N][4];
inline int cal(const string&s){
int t=0;
for(int i=1;i<s.size();i++)t=t*10+s[i]-'0';
return t;
}
inline void check(int l,int r){
for(x.set(0),y.set(0);l<=r;l++)x+=f[q[l][0]][q[l][1]].x,y+=f[q[l][0]][q[l][1]].y;
if(x.iszero()&&y.iszero())return;
puts("Infinity");
exit(0);
}
void cal(int x,int y,int d){
if(vis[x][y])return;
vis[x][y]=1;
q[in[x][y]=d][0]=x;
q[d][1]=y;
f[x][y].init(y);
for(int i=0;i<a[x];i++){
if(b[x][i]=="GO")f[x][y].go();
else if(b[x][i]=="LEFT")f[x][y].left();
else if(b[x][i]=="RIGHT")f[x][y].right();
else{
int A=cal(b[x][i]),B=f[x][y].d;
if(in[A][B]){
check(in[A][B],d);
cir[x][y]=1;
}else{
cal(A,B,d+1);
f[x][y].merge(f[A][B]);
if(cir[A][B])cir[x][y]=1;
}
if(cir[x][y])break;
}
}
in[x][y]=0;
}
int main(){
zero.set(0);
one.set(1);
cin>>n;
for(i=1;i<=n;i++){
cin>>a[i];
for(j=0;j<a[i];j++)cin>>b[i][j];
}
cal(1,0,1);
f[1][0].write();
return 0;
}
BZOJ2647 : [Neerc2011]Journey的更多相关文章
- bzoj AC倒序
Search GO 说明:输入题号直接进入相应题目,如需搜索含数字的题目,请在关键词前加单引号 Problem ID Title Source AC Submit Y 1000 A+B Problem ...
- 2632: [neerc2011]Gcd guessing game
2632: [neerc2011]Gcd guessing game Time Limit: 10 Sec Memory Limit: 128 MBSubmit: 144 Solved: 84[S ...
- CF721C. Journey[DP DAG]
C. Journey time limit per test 3 seconds memory limit per test 256 megabytes input standard input ou ...
- POJ2488A Knight's Journey[DFS]
A Knight's Journey Time Limit: 1000MS Memory Limit: 65536K Total Submissions: 41936 Accepted: 14 ...
- CF #374 (Div. 2) C. Journey dp
1.CF #374 (Div. 2) C. Journey 2.总结:好题,这一道题,WA,MLE,TLE,RE,各种姿势都来了一遍.. 3.题意:有向无环图,找出第1个点到第n个点的一条路径 ...
- POJ2488-A Knight's Journey(DFS+回溯)
题目链接:http://poj.org/problem?id=2488 A Knight's Journey Time Limit: 1000MS Memory Limit: 65536K Tot ...
- codeforces 721C C. Journey(dp)
题目链接: C. Journey time limit per test 3 seconds memory limit per test 256 megabytes input standard in ...
- A Knight's Journey 分类: POJ 搜索 2015-08-08 07:32 2人阅读 评论(0) 收藏
A Knight's Journey Time Limit: 1000MS Memory Limit: 65536K Total Submissions: 35564 Accepted: 12119 ...
- HDOJ-三部曲一(搜索、数学)- A Knight's Journey
A Knight's Journey Time Limit : 2000/1000ms (Java/Other) Memory Limit : 131072/65536K (Java/Other) ...
随机推荐
- 在windows下Apache安装配置
安装,从官网下载,安装即可. 配置遇到一些问题: 1. the requested operation has failed 这是因为安装后的文件目录没有没有写的权限.通过安全设置安装目录的所有 ...
- Entity Framework介绍
1.Entity Framework介绍 下图显示EF整体架构.现在我们来看看架构的各个组件: EDM(Entity Data Model): EDM由三个主要部分组成:概念模型,映射和存储模型.映射 ...
- [转]pyCharm最新2018激活码
https://blog.csdn.net/u014044812/article/details/78727496 因公司的需求,需要做一个爬取最近上映的电影.列车号.航班号.机场.车站等信息,所以需 ...
- linux 锁定重要文件 更改重要命令
锁定重要文件如下: chattr +i /etc/passwd /etc/shadow /etc/groupp /etc/gshadow 解锁文件 chattr -i /etc/passwd /etc ...
- Coolpy开源项目简介
1.Coolpy初识 基于ARDUINO的迷你操作系统.只需一个ARDUINO主板+Ethernet Shield即可运行.成本低,Coolpy主件以目前淘宝价只需要76元人民币. 2.Coolpy能 ...
- spark操作Kudu之读 - 使用DataFrame API
虽然我们可以通过上面显示的KuduContext执行大量操作,但我们还可以直接从默认数据源本身调用读/写API. 要设置读取,我们需要为Kudu表指定选项,命名我们要读取的表以及为表提供服务的Kudu ...
- P1057 传球游戏 dp
题目描述 上体育课的时候,小蛮的老师经常带着同学们一起做游戏.这次,老师带着同学们一起做传球游戏. 游戏规则是这样的:nn个同学站成一个圆圈,其中的一个同学手里拿着一个球,当老师吹哨子时开始传球,每个 ...
- day 69 orm操作之表关系,多对多,多对一(wusir总结官网的API)
对象 关系 模型 wusir博客地址orm官网API总结 django官网orm-API orm概要: ORM 跨表查询 class Book(models.Model): title = mod ...
- KMP算法2
给定一个主串s,一个子串sub,将主串中的所有子串替换成replaceStr,并将最终结果输出来. #include<stdio.h> #include<string.h> # ...
- Django之Models(二)
Django之Models(二) 创建一对多的关系 一个出版商可以有多本出版的书 一本书只有一个出版商 创建一对多的语法: 字段名= models.ForeignKey(关联表(类名),on_dele ...