第一次听说斯坦纳树这种东西

先dfs预处理出来dis[i][j][k]表示格子(i,j)向k方向转移能到哪,记忆话搜索预处理,注意如果有环的话特判一下

设f[i][j][x][y]表示复合机器人i-j在(x,y)生成需要推得步数,用spfa转移,因为时间比较紧所以优化spfa,把能转移的放进一个队列(准确的来讲是栈)按f排序,spfa转移到的放到另一个,每次取队首小的更新答案

#include<iostream>
#include<cstdio>
#include<cstring>
#include<queue>
using namespace std;
const int N=505,M=600003,dx[]={0,1,0,-1},dy[]={1,0,-1,0};
int n,h,w,f[10][10][N][N],v[N][N][5],ti,s[M+5],va[N*N],top,inf;
bool mk[N][N];
char c[N][N];
struct qwe
{
int x,y;
qwe(int X=0,int Y=0)
{
x=X,y=Y;
}
}dis[N][N][5],q[N*N],p[M+5];
qwe dfs(int x,int y,int k)
{
if(v[x][y][k]==ti)
return qwe(-1,0);
v[x][y][k]=ti;
if(dis[x][y][k].x!=0)
return dis[x][y][k];
int t=k;
if(c[x][y]=='C')
k=(k+1)%4;
if(c[x][y]=='A')
k=(k+3)%4;
int xx=x+dx[k],yy=y+dy[k];
if(xx>w||xx<1||yy>h||yy<1||c[xx][yy]=='x')
return dis[x][y][t]=qwe(x,y);
return dis[x][y][t]=dfs(xx,yy,k);
}
void spfa(int l,int r)
{
memset(s,0,sizeof(s));
memset(mk,0,sizeof(mk));
int mn=inf,mx=-inf;
for(int i=1;i<=top;i++)
{
s[va[i]]++;
if(va[i]<mn)
mn=va[i];
if(va[i]>mx)
mx=va[i];
mk[q[i].x][q[i].y]=1;
}
for(int i=mn+1;i<=mx;i++)
s[i]+=s[i-1];
for(int i=1;i<=top;i++)
p[s[va[i]]--]=q[i];
for(int i=1;i<=top;i++)
q[i]=p[top-i+1];
int ll=0,rr=0;
while(top||ll!=rr)
{
int now=ll%M+1;
qwe u;
if(ll==rr||(top&&f[l][r][q[top].x][q[top].y]<f[l][r][p[now].x][p[now].y]))
u=q[top--];
else
u=p[ll=now];
int x=u.x,y=u.y;
mk[x][y]=0;
for(int i=0;i<4;i++)
{
int tx=dis[x][y][i].x,ty=dis[x][y][i].y;
if(tx!=-1&&f[l][r][x][y]+1<f[l][r][tx][ty])
{
f[l][r][tx][ty]=f[l][r][x][y]+1;
if(!mk[tx][ty])
{
mk[tx][ty]=1;
p[rr=rr%M+1]=qwe(tx,ty);
}
}
}
}
}
int main()
{
scanf("%d%d%d",&n,&h,&w);
for(int i=1;i<=w;i++)
scanf("%s",c[i]+1);
memset(f,127/3,sizeof(f));
inf=f[0][0][0][0];cerr<<inf<<endl;
for(int i=1;i<=w;i++)
for(int j=1;j<=h;j++)
if(c[i][j]!='x')
for(int k=0;k<4;k++)
ti++,dis[i][j][k]=dfs(i,j,k);
for(int i=1;i<=w;i++)
for(int j=1;j<=h;j++)
if(c[i][j]>='0'&&c[i][j]<='9')
f[c[i][j]-48][c[i][j]-48][i][j]=0;
for(int l=1;l<=n;l++)
for(int x=1;x<=n-l+1;x++)
{
int y=x+l-1;top=0;
for(int i=1;i<=w;i++)
for(int j=1;j<=h;j++)
{
for(int k=x;k<y;k++)
if(f[x][y][i][j]>f[x][k][i][j]+f[k+1][y][i][j])
f[x][y][i][j]=f[x][k][i][j]+f[k+1][y][i][j];
if(f[x][y][i][j]!=inf)
{
q[++top]=qwe(i,j);
va[top]=f[x][y][i][j];
}
}
spfa(x,y);
}
int ans=inf;
for(int i=1;i<=w;i++)
for(int j=1;j<=h;j++)
if(f[1][n][i][j]<ans)
ans=f[1][n][i][j];
printf("%d\n",ans==inf?-1:ans);
return 0;
}

bzoj 3205: [Apio2013]机器人【dfs+斯坦纳树+spfa】的更多相关文章

  1. [Bzoj3205][Apio2013]机器人(斯坦纳树)(bfs)

    3205: [Apio2013]机器人 Time Limit: 15 Sec  Memory Limit: 128 MBSubmit: 977  Solved: 230[Submit][Status] ...

  2. BZOJ 3205 [Apio2013]机器人 ——斯坦纳树

    腊鸡题目,实在卡不过去. (改了一下午) 就是裸的斯坦纳树的题目,一方面合并子集,另一方面SPFA迭代求解. 优化了许多地方,甚至基数排序都写了. 还是T到死,不打算改了,就这样吧 #include ...

  3. BZOJ 2595: [Wc2008]游览计划 [DP 状压 斯坦纳树 spfa]【学习笔记】

    传送门 题意:略 论文 <SPFA算法的优化及应用> http://www.cnblogs.com/lazycal/p/bzoj-2595.html 本题的核心就是求斯坦纳树: Stein ...

  4. BZOJ 2595 [Wc2008]游览计划 ——斯坦纳树

    [题目分析] 斯坦纳树=子集DP+SPFA? 用来学习斯坦纳树的模板. 大概就是用二进制来表示树包含的点,然后用跟几点表示树的形态. 更新分为两种,一种是合并两个子集,一种是换根,换根用SPFA迭代即 ...

  5. BZOJ 5180 [Baltic2016]Cities(斯坦纳树)

    斯坦纳树的板子题. 斯坦纳树问题是组合优化问题,与最小生成树相似,是最短网络的一种. 最小生成树是在给定的点集和边中寻求最短网络使所有点连通. 而最小斯坦纳树允许在给定点外增加额外的点,使生成的最短网 ...

  6. bzoj 4006 [JLOI2015]管道连接——斯坦纳树

    题目:https://www.lydsy.com/JudgeOnline/problem.php?id=4006 除了模板,就是记录 ans[ s ] 表示 s 合法的最小代价.合法即保证 s 里同一 ...

  7. [APIO2013]机器人(斯坦纳树)

    题目描述 VRI(Voltron 机器人学会)的工程师建造了 n 个机器人.任意两个兼容的机 器人站在同一个格子时可以合并为一个复合机器人. 我们把机器人用 1 至 n 编号(n ≤ 9).如果两个机 ...

  8. [APIO2013]机器人[搜索、斯坦纳树]

    题意 题目链接 分析 记 g(d,x,y) 表示从 (x,y) 出发,方向为 d 到达的点,这个可以通过记忆化搜索求出,注意如果转移成环(此时向这个方向走没有意义)要特判. 记 f(l,r,x,y) ...

  9. [BZOJ3205][APIO2013]Robot(斯坦纳树)

    3205: [Apio2013]机器人 Time Limit: 15 Sec  Memory Limit: 128 MBSubmit: 1007  Solved: 240[Submit][Status ...

随机推荐

  1. Permutation Sequence(超时,排列问题)

    The set [1,2,3,…,n] contains a total of n! unique permutations. By listing and labeling all of the p ...

  2. 转 Linux里设置环境变量的方法(export PATH)

    1.动态库路径的设置 Linux下调用动态库和windows不一样.linux 可执行程序是靠配置文件去读取路径的,因此有些时候需要设置路径 具体操作如下export LD_LIBRARY_PATH= ...

  3. Antivius for Linux

    http://www.clamav.net/   https://www.avast.com/zh-cn/linux-server-antivirus   http://www.f-prot.com/ ...

  4. python的线程thread笔记

    python的线程是用thread和threading来实现的.其中利用threading会更好,因为thread没有线程保护,当主线程退出了之后,子线程也会被强行退出.threading支持守护线程 ...

  5. LoadRunner中存储表格参数------关联数组

    主要用到 web_reg_save_param_ex函数("Scope=All",), sprintf( CProdNo,"{CProdNo_%d}",i ); ...

  6. Swift String 一些经常用法

    直接上代码 //字符串 //1 推断字符串是否为空 var test1Str="" var test1Str2:String = String(); println("t ...

  7. Python 离线等价类

    离线等价类的概念见离线等价类 最近在清洗数据的时候涉及到要将相似度比较高的文件夹合并,特征比对得到是1:1的对,比如: (a,b),(c,d),(a,c)...,那么合并的时候就涉及到将这些等价的对合 ...

  8. Notepad++ 两个格式化插件

    格式化HTML--Tidy2 本来都可以通过Notepad++中的“插件>Plugin Manager>Show Plugin Manager>Tidy2” 这种方式来安装,不过内地 ...

  9. win7-64bit安装comtypes的问题

    Update 28/12/2014: Please download the latest comtypes 1.1.1 from https://pypi.python.org/pypi/comty ...

  10. 关于mybatis的 insert into select 命令未结束问题

    关于mybatis的 insert into select 命令未结束问题,最后以为是sql写错了,可是,在plsql运行又没问题.最后还是解决这个问题. 是设置问题. ### Cause: java ...