8.3.4 奥运

奥运

Time Limit: 1000/1000 MS (Java/Others) Memory Limit: 65536/65536 K (Java/Others)

Total Submission(s): 82 Accepted Submission(s): 40

Problem Description

北京迎来了第一个奥运会，我们的欢呼声响彻中国大地，所以今年的奥运金牌 day day up!
比尔盖兹坐上鸟巢里，手里摇着小纸扇，看的不亦乐乎，被俺们健儿的顽强拼搏的精神深深的感动了。反正我的钱也多的没地方放了，他对自己说，我自己也来举办一个奥运会，看谁的更火。不过他的奥运会很特别：
1 参加人员必须是中国人；
2 至少会加法运算（因为要计算本人获得的金牌数）
他 知道中国有很多的名胜古迹，他知道自己在t1 到 t2天内不可能把所有的地方都玩遍，所以他决定指定两个地方v1,v2，如果参赛员能计算出在t1到t2天（包括t1,t2）内从v1到v2共有多少种走 法(每条道路走需要花一天的时间，且不能在某个城市停留,且t1=0时的走法数为0)，那么他就会获得相应数量的金牌，城市的总数<=30,两个城 市间可以有多条道路
，每条都视为是不同的。

 

Input

本题多个case，每个case：
输入一个数字n表示有n条道路 0<n<10000
接下来n行每行读入两个数字 p1，p2 表示城市p1到p2有道路，并不表示p2到p1有道路 (0<=p1,p2<2^32)
输入一个数字k表示有k个参赛人员
接下来k行，每行读入四个数据v1,v2,t1,t2 (0<=t1,t2<10000)

 

Output

            对于每组数据中的每个参赛人员输出一个整数表示他获得的金牌数（mod 2008）

 

Sample Input

6
1 2
1 3
2 3
3 2
3 1
2 1
3
1 2 0 0
1 2 1 100
4 8 3 50

 

Sample Output

0
1506
0

思路：用map和邻接矩阵

路径条数就是把邻接矩阵乘起来

想想为什么？

我们怎么计算路径条数？

假设计算到a的路径条数=sigma({集合b：所有能到达a的点}的路径条数*(集合b的每一个对应元素到a)有多少条路)

#include <cmath>
#include <cstdio>
#include <algorithm>
#include <cstring>
#include <string>
#include <cstdlib>
#include <map>
using namespace std;

typedef long long ll;
typedef double dd;
const int maxn=10010;
const int mod=2008;
map<int,int> mp;
int n,v1,v2,V1,V2,query,T,t1,t2,m,x,y,ans;
struct matrix
{
    int m[31][31];
} c,r[maxn],base;

void close()
{
exit(0);
}

void print(matrix a)
{
    for (int i=1;i<=n;i++)
    {
        for (int j=1;j<=n;j++)
            printf("%d ",a.m[i][j]);
        puts("");
    }
}

matrix mul(matrix a,matrix b)
{
    memset(c.m,0,sizeof(c.m));
    for (int i=1;i<=n;i++)
        for (int j=1;j<=n;j++)
            for (int k=1;k<=n;k++)
            {
                c.m[i][j]+=a.m[i][k]*b.m[k][j];
                c.m[i][j] %= mod;
            }
    return c;
}

void init()
{

}

int main ()
{
    while (scanf("%d",&T)!=EOF)
    {
        n=0;
        mp.clear();
        memset(base.m,0,sizeof(base.m));
        while (T--)
        {
            scanf("%d %d",&x,&y);
            if (mp[x]==0)
            {
                n++;
                mp[x]=n;
            }
            if (mp[y]==0)
            {
                n++;
                mp[y]=n;
            }
            base.m[mp[x]][mp[y]]++;
            // X -> Y
        }
        m=10000; //Need to be changed!
        r[1]=base;
        for (int i=2;i<=m;i++)
            r[i]=mul(r[i-1],base);
        scanf("%d",&query);
        while (query--)
        {
            ans=0;
            scanf("%d %d %d %d",&v1,&v2,&t1,&t2);
            if (t1>t2 || t2==0)
            {
                ans=0;
            }
            else
            {
                if (mp[v1]==0 || mp[v2]==0)
                    ans=0;
                else
                {
                    V1=mp[v1];
                    V2=mp[v2];
                    for (int i=max(1,t1);i<=t2;i++)
                    {

                        ans+=r[i].m[V1][V2];
                        ans %= mod;
                    }
                }
            }
            printf("%d
",ans);
        }
    }
    return 0;
}