dp+spfa优化

最朴素的dp是dp[i][j]表示i->j的最短路,然后把所有pair(i,i)放到队列里跑spfa,但是这样被卡掉了,那么我们要优化一下

问题在于每次我们转移的时候要枚举i和j的邻居,这样会被两个连起来的菊花卡掉,那么我们希望一次只走一步,那么复杂度会大大降低,于是我们设一个状态g[i][j][k],表示当前在i,j,上一条出边的字符为k,这样我们让f和g交替转移,达到了每次只走一步的目标,然后就能过了,原先每次转移最坏O(m*m)?现在大概是O(n*m)的转移,状态数O(n*n*26)

#include<bits/stdc++.h>
using namespace std;
const int N = , inf = 0x3f3f3f3f;
struct edge {
int nxt, to, c;
} e[N * N << ];
struct node {
int u, v, c;
node(int u, int v, int c) : u(u), v(v), c(c) {}
};
int n, m, cnt = , q, last;
int head1[N], head2[N], f[N][N], g[N][N][], vis[N][N][];
void link1(int u, int v, int c)
{
e[++cnt].nxt = head1[u];
head1[u] = cnt;
e[cnt].to = v;
e[cnt].c = c;
}
void link2(int u, int v, int c)
{
e[++cnt].nxt = head2[u];
head2[u] = cnt;
e[cnt].to = v;
e[cnt].c = c;
}
void bfs()
{
queue<node> q;
memset(f, inf, sizeof(f));
memset(g, inf, sizeof(g));
for(int i = ; i <= n; ++i) vis[i][i][] = , f[i][i] = , q.push(node(i, i, ));
for(int i = ; i <= n; ++i)
for(int j = ; j < ; ++j) vis[i][i][j] = , g[i][i][j] = , q.push(node(i, i, j));
while(!q.empty())
{
node o = q.front();
q.pop();
int u = o.u, v = o.v;
vis[u][v][o.c] = ;
if(o.c == )
{
for(int i = head1[v]; i; i = e[i].nxt) if(g[u][e[i].to][e[i].c] > f[u][v] + )
{
g[u][e[i].to][e[i].c] = f[u][v] + ;
if(!vis[u][e[i].to][e[i].c])
{
vis[u][e[i].to][e[i].c] = ;
q.push(node(u, e[i].to, e[i].c));
}
}
}
else
{
for(int i = head2[u]; i; i = e[i].nxt) if(e[i].c == o.c && f[e[i].to][v] > g[u][v][o.c])
{
f[e[i].to][v] = g[u][v][o.c] + ;
if(!vis[e[i].to][v][])
{
vis[e[i].to][v][] = ;
q.push(node(e[i].to, v, ));
}
}
}
}
}
int main()
{
scanf("%d%d", &n, &m);
for(int i = ; i <= m; ++i)
{
int u, v;
char c[];
scanf("%d%d%s", &u, &v, c);
link1(u, v, c[] - 'a');
link2(v, u, c[] - 'a');
}
bfs();
scanf("%d%d", &q, &last);
q--;
while(q--)
{
int x;
scanf("%d", &x);
printf("%d\n", f[last][x] == inf ? - : f[last][x]);
last = x;
}
return ;
}

bzoj1138的更多相关文章

  1. [BZOJ1138][POI2009]Baj 最短回文路

    [BZOJ1138][POI2009]Baj 最短回文路 试题描述 N个点用M条有向边连接,每条边标有一个小写字母. 对于一个长度为D的顶点序列,回答每对相邻顶点Si到Si+1的最短回文路径. 如果没 ...

随机推荐

  1. 售货员的难题(codevs 2596)

    题目描述 Description 某乡有n个村庄(1<n<=15),有一个售货员,他要到各个村庄去售货,各村庄之间的路程s(0<s<1000)是已知的,且A村到B村与B村到A村 ...

  2. 【CF766D】Mahmoud and a Dictionary(并查集)

    题意:有n个单词,给定m个关系,每个关系要么表示单词a与单词b相同,要么表示单词a与单词b相反. 并且“相同”与“相反”有性质:若a与b相同,b与c相同,则a与c相同(从而单词的相同关系是等价关系): ...

  3. JS获取服务器时间并且计算距离当前指定时间差的函数

    项目中遇到了从服务器获取时间,现在记录一下方便以后查询: 1.后台代码:(创建一个date对象并以JSON的形式返回去) // 获取服务器时间 public String getNowServerTi ...

  4. MongoDB基本管理命令操作

    1. 查看所有数据库: show dbs 或: show databases 注意: 该命令不会显示新创建的空数据库,若想显示需要向空数据库插入一些数据. MongoDB中默认的数据库为test,若果 ...

  5. Linux中的进程与线程

    介绍了Linux下fork()创建进程以及使用pthread_create()创建线程的方法 1. 基于进程的斐波那契数列 在下面的代码中,由子进程进行斐波那契数列的输出,父进程要等待子进程输出完毕, ...

  6. 2018 11.1 PION 模拟赛

    期望:250  100+100+50 实际:210   80+100+30 期望:100   实际:80 最后:两个点T了.可能是求逆元的方法太慢了,也可能是闲的又加了一个快速乘的原因. #inclu ...

  7. windows上安装mysql

    安装mysql后 命令行闪退 查看服务 也没有MySQL服务启动 你安装了mysql没有,没有就先安装,安装好mysql以后,在bin目录下有个mysqld.exe,运行这个程序就可以添加mysql服 ...

  8. JUNIT -- springMVC的action进行单元测试

    原文:http://blog.csdn.net/gaopeng0071/article/details/49946575 我开发环境springMVC版本3.0.4 样例代码: package com ...

  9. Intellij IDEA远程调试tomcat

    1.windows系统 文件catalina.bat首行增加下面代码 set CATALINA_OPTS=-server -Xdebug -Xnoagent -Djava.compiler=NONE ...

  10. update语句执行卡死现象原因及解决方案

    https://blog.csdn.net/wpz0713/article/details/51499654 原因分析: 可能在PLSQL Developer执行update时没有commit,ora ...