hdu 2254 奥运 矩阵求路径方案数

#include <cstdio>
#include <cstring>
#include <cmath>
#include <map>
#include <string>
#include <iostream>
#include <algorithm>
using namespace std;
#define LL long long
#define maxn 35
map<string,int>mm;
const int mod=2008;
struct matrix{
    int f[maxn][maxn];
};
matrix mul(matrix a,matrix b,int n)
{
    matrix c;
    memset(c.f,0,sizeof(c.f));
    int i,j,k;

    for(k=1;k<=n;k++)
    {
        for(i=1;i<=n;i++)
        {
            if(!a.f[i][k])continue;
            for(j=1;j<=n;j++)
            {
                if(!b.f[k][j])continue;
                c.f[i][j]=(c.f[i][j]+a.f[i][k]*b.f[k][j])%mod;
            }
        }
    }
    return c;
}
matrix e,g[11111];
int main()
{
    int n;
    while(scanf("%d",&n)!=EOF)
    {
        int i,j,k,m=0,t1,t2;
        char a[22],b[22],v1[22],v2[22];
        mm.clear();

        memset(e.f,0,sizeof(e));
        for(i=0;i<n;i++)
        {
            scanf("%s%s",a,b);
            if(mm[a]==0)
                mm[a]=++m;
            if(mm[b]==0)
                mm[b]=++m;
            e.f[mm[a]][mm[b]]++;
        }
        memset(g[0].f,0,sizeof(g[0].f));
        for(i=1;i<=m;i++)
            g[0].f[i][i]=1;
        for(i=1;i<=10005;i++)
            g[i]=mul(g[i-1],e,m);
        scanf("%d",&k);
        while(k--)
        {
            scanf("%s%s%d%d",v1,v2,&t1,&t2);
            if(mm[v1]==0||mm[v2]==0||(t1==0&&t2==0))
            {
                printf("0\n");
                continue;
            }
            if(t1>t2)
                swap(t1,t2);
            int ans=0;
            for(i=t1;i<=t2;i++)
            {
                if(i==0)continue;//当i=0时，必定为0，即使v1=v2且存在e.f[v1][v2]>0
                ans=(ans+g[i].f[mm[v1]][mm[v2]])%mod;
            }
            printf("%d\n",ans);
        }
    }
    return 0;
}
/*
    矩阵的基本：
    矩阵f中f[i][j]表示i到j可以走的条数
    那么f^2中,f^2[i][j]表示i到j要经过天的方案数
    可以认作是状态的转移。。
    
    特例：
    Input:
    1
    1 1
    1 
    1 1 0 1
    Ouput:
    1
*/