码迷，mamicode.com
首页 > 其他好文 > 详细

HDOJ 4686 Arc of Dream 矩阵快速幂

时间：2014-11-05 06:58:32      阅读：347      评论：0      收藏：0      [点我收藏+]

标签：des   style   blog   http   io   color   ar   os   java   

矩阵快速幂:


根据关系够建矩阵 , 快速幂解决.

bubuko.com,布布扣

Arc of Dream

Time Limit: 2000/2000 MS (Java/Others)    Memory Limit: 65535/65535 K (Java/Others)
Total Submission(s): 2164    Accepted Submission(s): 680


Problem Description
An Arc of Dream is a curve defined by following function:
bubuko.com,布布扣

where
a0 = A0
ai = ai-1*AX+AY
b0 = B0
bi = bi-1*BX+BY
What is the value of AoD(N) modulo 1,000,000,007?
 


Input
There are multiple test cases. Process to the End of File.
Each test case contains 7 nonnegative integers as follows:
N
A0 AX AY
B0 BX BY
N is no more than 1018, and all the other integers are no more than 2×109.
 


Output
For each test case, output AoD(N) modulo 1,000,000,007.
 


Sample Input
1 1 2 3 4 5 6 2 1 2 3 4 5 6 3 1 2 3 4 5 6
 


Sample Output
4 134 1902
 


Author
Zejun Wu (watashi)
 


Source
2013 Multi-University Training Contest 9
 


[cpp] view plaincopybubuko.com,布布扣bubuko.com,布布扣
    1. #include <iostream>  
    2. #include <cstdio>  
    3. #include <cstring>  
    4. #include <algorithm>  
    5. #include <cmath>  
    6.   
    7. using namespace std;  
    8.   
    9. typedef long long int LL;  
    10.   
    11. const LL mod=1000000007LL;  
    12.   
    13. struct Matrix  
    14. {  
    15.     int x,y;  
    16.     LL m[6][6];  
    17.     Matrix() {x=y=5;memset(m,0,sizeof(m));}  
    18.     void one()  
    19.     {  
    20.         for(int i=0;i<5;i++) m[i][i]=1LL;  
    21.     }  
    22.     void show()  
    23.     {  
    24.         cout<<x<<" * "<<y<<endl;  
    25.         for(int i=0;i<x;i++)  
    26.         {  
    27.             for(int j=0;j<y;j++)  
    28.                 cout<<m[i][j]<<",";  
    29.             cout<<endl;  
    30.         }  
    31.     }  
    32. };  
    33.   
    34. Matrix Mul(Matrix& a,Matrix& b)  
    35. {  
    36.     Matrix ret;  
    37.     ret.x=a.x; ret.y=b.y;  
    38.     for(int i=0;i<a.x;i++)  
    39.     {  
    40.         for(int j=0;j<b.y;j++)  
    41.         {  
    42.             LL temp=0;  
    43.             for(int k=0;k<b.y;k++)  
    44.             {  
    45.                 temp=(temp+(a.m[i][k]*b.m[k][j])%mod)%mod;  
    46.             }  
    47.             ret.m[i][j]=temp%mod;  
    48.         }  
    49.     }  
    50.     return ret;  
    51. }  
    52.   
    53. Matrix quickPow(Matrix m,LL x)  
    54. {  
    55.     Matrix e;  
    56.     e.one();  
    57.     while(x)  
    58.     {  
    59.         if(x&1LL) e=Mul(e,m);  
    60.         m=Mul(m,m);  
    61.         x/=2LL;  
    62.     }  
    63.     return e;  
    64. }  
    65.   
    66. LL n,A0,B0,AX,AY,BX,BY;  
    67.   
    68. Matrix init_matrix()  
    69. {  
    70.     Matrix ret;  
    71.     ret.m[0][0]=1;  
    72.     ret.m[1][0]=AY; ret.m[1][1]=AX;  
    73.     ret.m[2][0]=BY; ret.m[2][2]=BX;  
    74.     ret.m[3][0]=(BY*AY)%mod; ret.m[3][1]=(AX*BY)%mod;  
    75.     ret.m[3][2]=(BX*AY)%mod; ret.m[3][3]=(AX*BX)%mod;  
    76.     ret.m[4][3]=1LL; ret.m[4][4]=1LL;  
    77.     return ret;  
    78. }  
    79.   
    80. Matrix Beg()  
    81. {  
    82.     Matrix beg;  
    83.     beg.m[0][0]=1;  
    84.     beg.m[1][0]=A0;  
    85.     beg.m[2][0]=B0;  
    86.     beg.m[3][0]=A0*B0%mod;  
    87.     return beg;  
    88. }  
    89.   
    90. int main()  
    91. {  
    92.     while(cin>>n)  
    93.     {  
    94.         cin>>A0>>AX>>AY>>B0>>BX>>BY;  
    95.         A0=A0%mod; AX=AX%mod; AY=AY%mod;  
    96.         B0=B0%mod; BX=BX%mod; BY=BY%mod;  
    97.         Matrix m=init_matrix();  
    98.         m=quickPow(m,n);  
    99.         Matrix beg=Beg();  
    100.         LL ans=0;  
    101.         for(int i=0;i<5;i++)  
    102.             ans=(ans+beg.m[i][0]*m.m[4][i]%mod)%mod;  
    103.         cout<<ans<<endl;  
    104.     }  
    105.     return 0;  
    106. } 

HDOJ 4686 Arc of Dream 矩阵快速幂

标签：des   style   blog   http   io   color   ar   os   java   

原文地址：http://www.cnblogs.com/yuyanbian/p/4075307.html
踩
(0)
赞
(0)
   
举报
评论 一句话评论（0）
登录后才能评论！
分享档案
更多>
2021年07月29日 (22)
2021年07月28日 (40)
2021年07月27日 (32)
2021年07月26日 (79)
2021年07月23日 (29)
2021年07月22日 (30)
2021年07月21日 (42)
2021年07月20日 (16)
2021年07月19日 (90)
2021年07月16日 (35)
周排行
mamicode.com排行更多图片
更多
友情链接
兰亭集智   国之画   百度统计   站长统计   阿里云   chrome插件   新版天听网
关于我们 - 联系我们 - 留言反馈
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码！