BZOJ2272: [Usaco2011 Feb]Cowlphabet 奶牛文字

n<=250个大写字母和m<=250个小写字母，给p<=200个合法相邻字母，求用这些合法相邻字母的规则和n+m个字母能合成多少合法串，答案mod 97654321.

什么鬼膜数。。

f(i,j,k)--i个大写字母，j个小写字母，最后一个字母是k，，其中k是小写字母，p是能接在k前面的任意字母，k是大写字母的话同理。

这样复杂度是n*m*26*26?假的！i,j确定后只有p种合法转移，所以就n*m*p。

trick：记得算i=0和j=0的情况！！！

#include<cstring>
#include<cstdlib>
#include<cstdio>
//#include<assert.h>
//#include<time.h>
#include<math.h>
#include<algorithm>
//#include<iostream>
using namespace std;

bool isdigit(char c) {return c>='0' && c<='9';}
int qread()
{
    char c;int s=0,f=1;while (!isdigit(c=getchar())) f=(c=='-'?-1:1);
    do s=s*10+c-'0'; while (isdigit(c=getchar())); return s*f;
}

const int mod=97654321;
int da,xiao,m;
int f[255][255][55];
int id(char c) {return (c>='a' && c<='z')?c-'a':c-'A'+26;}
bool isalpha(char c) {return (c>='a' && c<='z') || (c>='A' && c<='Z');}
struct Edge{int to,next;}edge[255];int first[55],le=2;
void in(int x,int y) {Edge &e=edge[le];e.to=y;e.next=first[x];first[x]=le++;}
int main()
{
    da=qread(),xiao=qread(),m=qread();
    memset(first,0,sizeof(first));
    for (int i=1;i<=m;i++)
    {
        char c1,c2;
        while (!isalpha(c1=getchar()));
        c2=getchar();
        in(id(c2),id(c1));
    }
    for (int i=0;i<26;i++) f[0][1][i]=1,f[1][0][i]=0;
    for (int i=26;i<52;i++) f[0][1][i]=0,f[1][0][i]=1;
    for (int i=0;i<=da;i++)
        for (int j=0;j<=xiao;j++) if (i+j>=2)
        {
            if (j) for (int now=0;now<26;now++)
            {
                f[i][j][now]=0;
                for (int k=first[now];k;k=edge[k].next)
                {
                    const Edge &e=edge[k];
                    f[i][j][now]+=f[i][j-1][e.to],
                    f[i][j][now]-=f[i][j][now]>=mod?mod:0;
                }
            }
            if (i) for (int now=26;now<52;now++)
            {
                f[i][j][now]=0;
                for (int k=first[now];k;k=edge[k].next)
                {
                    const Edge &e=edge[k];
                    f[i][j][now]+=f[i-1][j][e.to],
                    f[i][j][now]-=f[i][j][now]>=mod?mod:0;
                }
            }
        }
    int ans=0;
    for (int i=0;i<52;i++) ans+=f[da][xiao][i],ans-=ans>=mod?mod:0;
    printf("%d
",ans);
    return 0;
}