围巾的ACM博客

思路：矩阵快速幂的入门题 #include #include #include using namespace std; struct Mat { int a[2][2]; //矩阵大小 }; int n; const int mod = 10000; Mat mul(Mat a,Mat b) { Mat t; memset(t.a,0,sizeof(t.a)); for(i

思路：看k的范围显然不能直接遍历来加，留意到S有一个性质            当k为偶数的时候 Sk = A+A^2+....A^(K/2)+A^(K/2)*(A+A^2+....A^(k/2))            当k为奇数的时候 Sk = A+A^2+....A^(K/2)+A^((K+1)/2) + A^((K+1)/2)*(A+A^2....A^(K/2))          

思路：简单的矩阵快速幂 #include #include #include using namespace std; int n,mod; struct Mat { int a[31][31]; //矩阵大小 }; Mat mul(Mat a,Mat b) { Mat t; memset(t.a,0,sizeof(t.a)); for(int i = 0;i<n;i++)

思路：显然的有f(n)=A^n，那么把g(i)代入后有sum = A^b+A^(k+b)+A^(2k+b)+...A^((n-1)k+b)            提取A^b之后sum = A^b*(A^k^0+A^k^1+A^k^2+.....A^k^(n-1))，将A^k设为B的话那么就是一个简单的二分求和过程，注意最后要加上B^0 坑点：虽然最后结果不会爆long long，可是中途乘法过