bzoj 1138: [POI2009]Baj 最短回文路

额，，貌似网上的题解都说超时之类的。

然而我这个辣鸡在做的时候不知道在想什么，连超时的都不会。

超时的大概是这样的，f[x][y]表示x到y的最短回文路，然后更新的话就是 f[x][y]更新到 f[a][b] 当x->a,y->b且边的颜色是一样的。

然后yy了一下为什么会超时呢。。。。

然后想到了一个sb的情况。。两个菊花图连起来。。。这样的复杂度就呵呵呵了。。每次用f[x][y]更新的话都要枚举一个颜色，然后把另一边的颜色找出来，这样显然是要爆炸的。

为了避免这种sb情况，就引入了一个叫题解的神奇的东西。

题解多设了一个东西，g[x][y][col]表示x到y且除了和y相连的那条颜色为col的边之外是回文最短路的长度。

然后有了g，就可以在一遍更新过g，f在另一边直接拿来用了。神奇2333

#include<iostream>
#include<iostream>
#include<algorithm>
#include<cstdio>
#include<cstring>
#define inf 0x3f3f3f3f
using namespace std;

const int maxn=100005;
const int QWQ=505;

struct edge{
    int to,from,next;
}e[maxn<<1],E[maxn<<1];
int head[QWQ],cnt,HEAD[QWQ],CNT;
void insert(int x, int y)
{
    e[++cnt].next=head[x]; e[cnt].to=y; e[cnt].from=x; head[x]=cnt;
    E[++CNT].next=HEAD[y]; E[CNT].to=x; E[CNT].from=y; HEAD[y]=CNT;
}

int l,r=1;
int n,m,map[QWQ][QWQ],f[QWQ][QWQ],g[QWQ][QWQ][30];
struct node{
    int from,to,col;
}q[QWQ*QWQ*30];
void bfs()
{
    while (l<r)
    {
        node u=q[l++];
        int x=u.from,y=u.to;
        if (!u.col)
        {
            for (int i=head[y];i;i=e[i].next)
            {
                if (f[x][y]+1<g[x][e[i].to][map[y][e[i].to]])
                {
                    g[x][e[i].to][map[y][e[i].to]]=f[x][y]+1;
                    q[r++]=(node){x,e[i].to,map[y][e[i].to]};
                }
            }
        }
        else
        {
            for (int i=HEAD[x];i;i=E[i].next)
            {
                if (g[x][y][map[E[i].to][x]]+1<f[E[i].to][y])
                {
                    f[E[i].to][y]=g[x][y][map[E[i].to][x]]+1;
                    q[r++]=(node){E[i].to,y,0};
                }
            }
        }
    }
}
int main()
{
    scanf("%d%d",&n,&m);
    memset(f,0x3f,sizeof(f)); memset(g,0x3f,sizeof(g));
    for (int i=1; i<=m; i++)
    {
        int x,y; char Orz[2];
        scanf("%d%d%s",&x,&y,Orz);
        insert(x,y);
        map[x][y]=Orz[0]-'a'+1;
        f[x][y]=1;
        q[r++]=(node){x,y,0};
    }
    for (int i=1; i<=n; i++) f[i][i]=0,q[r++]=(node){i,i,0};
    bfs();
    int T; scanf("%d",&T);
    int x,y; scanf("%d",&x);
    for (int i=2; i<=T; i++)
    {
        scanf("%d",&y);
        if (f[x][y]!=inf) cout<<f[x][y]<<endl;
            else puts("-1");
        x=y;
    }
    return 0;
}