Additions HNUST 1713(矩阵快速幂模板 )

Ada is 5 years old, and she is learning additions. Her father writes some exercises for her.

1+11=?

2+22=?

3+33=?

“It is easy”, she writes the answer correctly. “Try to answer the following questions”,

1+11+111=?

2+22+222+2222=?

Ada scratches her head. “It’s too difficult”.

So, help Ada to compute the value of equation a+aa+aaa+... which have n items. Since the answer may be quite large, you have to module it by an integer m.

输入
Input contains at most 1000 test cases. Each test case contains only one line.

Each line will contain three integers a, n and m. (1<=a<=9, 1<=n<231-1, 1<=m<=100007). Process to end of file.

输出
For each test cases, output the value of (a+aa+aaa+... .)%m。

样例输入
1 1 13
1 2 13
1 3 13
2 2 13
2 3 13
样例输出
1
12
6
11

思路：刚开始看到这个题目想的是纯模拟，但是高达了2^31-1，所以不可以，后来队友告诉我，可以构造出来一个递推的式子

|1 10 1| |fn-1| |fn|
|0 10 1| |xn-1|=|xn|
|0 0 1| |a | |a|

|1 10 1|^n-1 |f1| |fn|
|0 10 1| |x1|=|xn|
|0 0 1| |a | |a|

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
/*1 10 1||fn-1|   |fn|
0 10 1||xn-1|=|xn|
0 0   1||a     |  |a|
*/
int a,b,MOD;
struct mat
{ll a[3][3]; //数组a(一般为存在幂的数组)  3为方阵的长宽
};
mat mat_mul(mat x,mat y)
{mat res;memset(res.a,0,sizeof(res.a));    //初始化for(int i=0;i<3;i++)for(int j=0;j<3;j++)for(int k=0;k<3;k++)res.a[i][j]=(res.a[i][j]+x.a[i][k]*y.a[k][j])%MOD;  //矩阵乘法return res;
}
void mat_pow(int n)
{mat c,res;c.a[0][0]=1,c.a[0][1]=10,c.a[0][2]=1;c.a[1][0]=0,c.a[1][1]=10,c.a[1][2]=1;c.a[2][0]=0,c.a[2][1]=0,c.a[2][2]=1;     //初始化赋值需要幂乘的数组memset(res.a,0,sizeof(res.a)); for(int i=0;i<3;i++) res.a[i][i]=1;      //初始化 a的n次方的数组且对角线均为1while(n){if(n&1) res=mat_mul(res,c);c=mat_mul(c,c);n=n>>1;}                                       //矩阵快速幂的运算过程int f1=a%MOD,x=a;                       //数组b的构成printf("%lld\n",(res.a[0][0]*f1+res.a[0][1]*x+res.a[0][2]*a)%MOD);
}
int main()
{int n;while(scanf("%d%d%d",&a,&n,&MOD)==3){mat_pow(n-1);}return 0;
}

Additions HNUST 1713(矩阵快速幂模板 )相关推荐

POJ3070 矩阵快速幂模板
题目:http://poj.org/problem?id=3070 矩阵快速幂模板.mod写到乘法的定义部分就行了. 别忘了 I ( ) 和 i n i t ( ) 要传引用! #include< ...
快速幂+矩阵快速幂模板
快速..运算快速幂运用位运算代码分析矩阵快速幂题目分析代码拓一..: 快速幂运用位运算强大的位运算把我搞得蒙蒙的理解了之后我表示很喜欢!!! 代码 int power(int ...
51nod 1113 矩阵快速幂模板题
1113 矩阵快速幂基准时间限制:3 秒空间限制:131072 KB 分值: 40 难度:4级算法题收藏关注给出一个N * N的矩阵,其中的元素均为正整数.求这个矩阵的M次方.由于M次方的计 ...
How many ways?? - hdu2157（矩阵快速幂-模板）
分析:求Map^k,刚开始没有用快速幂,TLE了代码如下: =================================================================== ...
(转)矩阵快速幂模板
大佬博客:https://blog.csdn.net/baidu_23081367/article/details/52347256 代码: const int mat_size = 5;//矩阵大小 ...
51nod 1113 矩阵快速幂模板
给出一个N * N的矩阵,其中的元素均为正整数.求这个矩阵的M次方.由于M次方的计算结果太大,只需要输出每个元素Mod (10^9 + 7)的结果. Input 第1行:2个数N和M,中间用空格分隔. ...
43行代码AC_HDU-2604 Queuing（矩阵快速幂，附详细的知识讲解、模板例题）
一道经典的矩阵快速幂模板题. 传送门1-->快速幂基本思想传送门2-->矩阵快速幂讲解(教主传授) 传送门3.1-->HDU-1575(经典矩阵快速幂模板题1) 传送门3.2--& ...
算法分类整理+模板①：矩阵快速幂
一直有一个想法,感觉自己很多基础算法不是很扎实,想要找个机会写一些算法的整理,顺便自己总结一些实用的模板. 最近偶然在训练赛中连续做了2道思维+矩阵快速幂的题目,碰巧有时间,就以矩阵快速幂作为这个系列 ...
快速幂 + 矩阵快速幂
快速幂 1 #include<iostream> 2 #include<algorithm> 3 #include<cstring> 4 #define LL lo ...
蓝桥杯算法提高递推求值（矩阵快速幂）详解
传送门问题描述已知递推公式: F(n, 1)=F(n-1, 2) + 2F(n-3, 1) + 5, F(n, 2)=F(n-1, 1) + 3F(n-3, 1) + 2F(n-3, 2) + 3 ...

Additions HNUST 1713(矩阵快速幂模板 )

Additions HNUST 1713(矩阵快速幂模板 )相关推荐

最新文章

热门文章