LightOJ - 1132 Summing up Powers 矩阵快速幂

题目大意:求(1^K + 2^K + 3K + … + N^K) % 2^32

解题思路:

借用别人的图

可以先打表,求出Cnm,用杨辉三角可以快速得到

#include<cstdio>
typedef unsigned long long ll;
const int N = 55;
const ll mod = (1LL << 32);
struct Matrix{
    ll mat[N][N];
}A, B, tmp;

ll n, num[N];
ll C[N][N];
int K;

void init2() {
    C[0][0] = 1;
    for(int i = 1; i <= 50; i++) {
        C[i][i] = C[i][0] = 1;
        for(int j = 1; j < i; j++)
            C[i][j] = (C[i-1][j-1] + C[i-1][j]) % mod;
    }
}

Matrix matMul(const Matrix &x, const Matrix &y) {
    for(int i = 0; i < K + 2; i++)
        for(int j = 0; j < K + 2; j++) {
            tmp.mat[i][j] = 0;
            for(int k = 0; k < K + 2; k++) {
                tmp.mat[i][j] = (tmp.mat[i][j] + x.mat[i][k] * y.mat[k][j]) % mod;
            }
        }
    return tmp;
}
void solve() {
    while(n) {
        if(n & 1)
            B = matMul(B,A);
        A = matMul(A,A);
        n >>= 1;
    }
}

void init() {
    for(int i = 0; i < K + 2; i++)
        for(int j = 0; j < K + 2; j++) {
            B.mat[i][j] = A.mat[i][j] = 0;
            if(i == j)
                B.mat[i][j] = 1;
        }

    A.mat[0][0] = 1;
    for(int i = 1; i < K + 2; i++)
        A.mat[i][0] = A.mat[i][1] = C[K][i-1];
    for(int i = 2; i < K + 2; i++)
        for(int j = i; j < K + 2; j++) {
            A.mat[j][i] = C[K-i+1][j-i];
        }
}

int main() {
    int test, cas = 1;
    scanf("%d", &test);
    init2();
    while(test--) {
        scanf("%lld%d", &n, &K);
        printf("Case %d: ", cas++);
        ll ans = ( (n % mod) * ( (n + 1) % mod) / 2) % mod;
        if(K == 1) {
            printf("%lld\n", ans);
            continue;
        }
        init();

        n--;
        solve();
        ans = 0;
        for(int i = 0; i < K + 2; i++)
            ans = (ans + B.mat[i][0] ) % mod;
        printf("%lld\n", ans);
    }
    return 0;
}
时间: 2024-10-10 00:38:03

LightOJ - 1132 Summing up Powers 矩阵快速幂的相关文章

LightOJ 1132 - Summing up Powers 矩阵快速幂+排列组合

链接:http://lightoj.com/volume_showproblem.php?problem=1132 1132 - Summing up Powers PDF (English) Statistics Forum Time Limit: 2 second(s) Memory Limit: 32 MB Given N and K, you have to find (1K + 2K + 3K + ... + NK) % 232 Input Input starts with an i

LightOJ 1132 Summing up Powers(矩阵快速幂+二项式定理)

LightOJ 1132 题意: 给出N(1≤N≤1015),K(0≤K≤50): 计算:(1K+2K+3K+...+NK)%232. 思路: 根据二项式定理,我们可以得到: (n+1)k=C0knk+C1knk?1+C2knk?2+...+Ck?1kn+1 于是可以构造变化矩阵: ???????????C0k00...01C1kC0k?10...00C2kC1k?1C0k?2...00C3kC2k?1C1k?2...00..................111...11???????????

LightOJ 1132 Summing up Powers:矩阵快速幂 + 二项式定理

题目链接:http://lightoj.com/volume_showproblem.php?problem=1132 题意: 给定n.k,求(1K + 2K + 3K + ... + NK) % 232. 题解: 设sum(i) = 1K + 2K + 3K + ... + iK 所以要从sum(1)一直推到sum(n). 所以要找出sum(i)和sum(i+1)之间的关系: 好了可以造矩阵了. (n = 6时) 矩阵表示(大小为 1 * (k+2)): 初始矩阵start: 也就是: 特殊矩

LightOJ 1070 - Algebraic Problem 推导+矩阵快速幂

http://www.lightoj.com/volume_showproblem.php?problem=1070 思路:\({(a+b)}^n =(a+b){(a+b)}^{n-1} \) \((ab)C_{n}^{r}a^{n-r}b{r} = C_{n+2}^{r}a^{n-r+2}b{r} - a^{n+2} - b^{n+2} \) 综上\( f(n) = (a+b)f(n-1)-(ab)f(n-2) \) /** @Date : 2016-12-19-19.53 * @Author

LightOJ 1070 Algebraic Problem (推导+矩阵快速幂)

题目链接:LightOJ 1070 Algebraic Problem 题意:已知a+b和ab的值求a^n+b^n.结果模2^64. 思路: 1.找递推式 得到递推式之后就是矩阵快速幂了 注意:模2^64,定义成unsigned long long 类型,因为无符号类型超过最大范围的数与该数%最大范围 的效果是一样的. AC代码: #include<stdio.h> #include<string.h> #define LL unsigned long long struct Ma

矩阵快速幂刷题系列

来源自http://blog.csdn.net/chenguolinblog/article/details/10309423 hdu 1575 Tr A Time Limit: 1000/1000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Others)Total Submission(s): 5587    Accepted Submission(s): 4200 Problem Description A为一个方阵,则Tr

HDU 1757 A Simple Math Problem (矩阵快速幂)

[题目链接]:click here~~ [题目大意]: If x < 10 f(x) = x. If x >= 10 f(x) = a0 * f(x-1) + a1 * f(x-2) + a2 * f(x-3) + -- + a9 * f(x-10); 问f(k)%m的值. [思路]:矩阵快速幂,具体思路看代码吧,注意一些细节. 代码: #include<bits/stdc++.h> using namespace std; typedef long long LL; const

Codeforces Round #291 (Div. 2) E - Darth Vader and Tree (DP+矩阵快速幂)

这题想了好长时间,果断没思路..于是搜了一下题解.一看题解上的"快速幂"这俩字,不对..这仨字..犹如醍醐灌顶啊...因为x的范围是10^9,所以当时想的时候果断把dp递推这一方法抛弃了.我怎么就没想到矩阵快速幂呢.......还是太弱了..sad..100*100*100*log(10^9)的复杂度刚刚好. 于是,想到了矩阵快速幂后,一切就变得简单了.就可以把距离<=x的所有距离的点数都通过DP推出来,然后一个快速幂就解决了. 首先DP递推式很容易想到.递推代码如下: for(

POJ 3233 - Matrix Power Series ( 矩阵快速幂 + 二分)

POJ 3233 - Matrix Power Series ( 矩阵快速幂 + 二分) #include <cstdio> #include <cstring> #include <algorithm> using namespace std; typedef long long LL; #define MAX_SIZE 30 #define CLR( a, b ) memset( a, b, sizeof(a) ) int MOD = 0; int n, k; st