zoukankan      html  css  js  c++  java
  • HDU 2157 How many ways??(矩阵快速幂)

    How many ways??

    Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Others)
    Total Submission(s): 1912    Accepted Submission(s): 692


    Problem Description
    春天到了, HDU校园里开满了花, 姹紫嫣红, 非常美丽. 葱头是个爱花的人, 看着校花校草竞相开放, 漫步校园, 心情也变得舒畅. 为了多看看这迷人的校园, 葱头决定, 每次上课都走不同的路线去教室, 但是由于时间问题, 每次只能经过k个地方, 比方说, 这次葱头决定经过2个地方, 那他可以先去问鼎广场看看喷泉, 再去教室, 也可以先到体育场跑几圈, 再到教室. 他非常想知道, 从A 点恰好经过k个点到达B点的方案数, 当然这个数有可能非常大, 所以你只要输出它模上1000的余数就可以了. 你能帮帮他么?? 你可决定了葱头一天能看多少校花哦
     
    Input
    输入数据有多组, 每组的第一行是2个整数 n, m(0 < n <= 20, m <= 100) 表示校园内共有n个点, 为了方便起见, 点从0到n-1编号,接着有m行, 每行有两个整数 s, t (0<=s,t<n) 表示从s点能到t点, 注意图是有向的.接着的一行是两个整数T,表示有T组询问(1<=T<=100),
    接下来的T行, 每行有三个整数 A, B, k, 表示问你从A 点到 B点恰好经过k个点的方案数 (k < 20), 可以走重复边。如果不存在这样的走法, 则输出0
    当n, m都为0的时候输入结束
     
    Output
    计算每次询问的方案数, 由于走法很多, 输出其对1000取模的结果
     
    Sample Input
    4 4 0 1 0 2 1 3 2 3 2 0 3 2 0 3 3 3 6 0 1 1 0 0 2 2 0 1 2 2 1 2 1 2 1 0 1 3 0 0
     
    Sample Output
    2 0 1 3
     

    解题思路:

      这道题就是模板题,图论中有一个结论就是说,我现在有一个n个顶点构成的邻接矩阵,那么我们要得到某个顶点救过k条路径到达某个顶点的方案数目,其实就是所说的矩阵的快速幂了。输出经过k次幂后的邻接矩阵。

     1 # include<iostream>
     2 # include<cstdio>
     3 # include<cstring>
     4 
     5 using namespace std;
     6 
     7 int n,m;
     8 
     9 struct matrix
    10 {
    11     int p[25][25];
    12 };
    13 
    14 struct matrix multi( matrix a,matrix b )//连个矩阵相乘
    15 {
    16     matrix c;
    17     for( int i=0;i<n;i++ )
    18     {
    19         for( int j=0;j<n;j++ )
    20         {
    21             c.p[i][j]=0;
    22             for( int k=0;k<n;k++ )
    23                 c.p[i][j]=(c.p[i][j]+a.p[i][k]*b.p[k][j])%1000;
    24         }
    25     }
    26 
    27     return c;
    28 }
    29 
    30 matrix My_power( matrix a, int k )
    31 {
    32 
    33     matrix b;
    34     for( int i=0;i<n;i++)
    35         for( int j=0;j<n;j++)
    36             if( i==j )
    37                 b.p[i][j]=1;
    38             else
    39                 b.p[i][j]=0;
    40     while( k )
    41     {
    42         if(k%2==1)
    43             b=multi(b,a);
    44         k=k/2;
    45         a=multi(a,a);
    46     }
    47     return b;
    48 }
    49 
    50 int main(void)
    51 {
    52     int S,E,k;
    53     matrix a,b;
    54     while(scanf("%d%d",&n,&m)!=EOF&&(n!=0||m!=0))
    55     {
    56         memset(a.p,0,sizeof(a.p));
    57         for( int i=0;i<m;i++)
    58         {
    59             int x1,x2;
    60             scanf("%d%d",&x1,&x2);
    61             a.p[x1][x2]=1;
    62         }
    63         int T; scanf("%d",&T);
    64         while(T--)
    65         {
    66             scanf("%d%d%d",&S,&E,&k);
    67             b = My_power(a,k);
    68             while( b.p[S][E] < 0 )//以后碰到取模的情况记得添加
    69                 b.p[S][E] += 1000;
    70             printf("%d
    ",b.p[S][E]);
    71         }
    72     }
    73     return 0;
    74 }

    代码:

  • 相关阅读:
    Java hello/hi的简单的网络聊天程序
    案例分析:设计模式与代码的结构特性
    网络相关的命令工具研究报告
    如何提高程序员的键盘使用效率
    分析一套源代码的代码规范和风格并讨论如何改进优化代码
    用例建模Use Case Modeling
    结合工程实践选题调研分析同类软件产品
    领域建模
    SecureCRT 连接虚拟机Linux
    用css固定div层在页面顶部和底部(兼容IE6)
  • 原文地址:https://www.cnblogs.com/wikioibai/p/4508971.html
Copyright © 2011-2022 走看看