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 integer T (≤ 200), denoting the number of test cases.

Each case contains two integers N (1 ≤ N ≤ 1015) and K (0 ≤ K ≤ 50) in a single line.

Output

For each case, print the case number and the result.

Sample Input

Output for Sample Input


3

3 1

4 2

3 3


Case 1: 6

Case 2: 30

Case 3: 36

题意:给n和k  计算那串公式的值。

做法:

找出 1^k 怎么推到2^k 再推到n^k的方法,再开一维记录总的值,就ok了。

初始矩阵

1^ 0  1^1 1^2 1^3 .....1^k    总

构造矩阵:

C(0,0) C(0,1) C(0,2) C(0,3)...C(0,k-1)      C(0,k)      0

0          C(1,1) C(1,2) C(1,3)...C(1,k-1)      C(1,k)      0

......

0          0           0         0            C(k-1,k-1)   C(k-1,k)  0

0          0           0         0            0                   C(k,k)     1

0          0           0         0            0                   0             1

C(n,m)=C(n-1,m)+C(n-1,m-1);

#include<stdio.h>
#include<string.h>
#define Matr 55 //矩阵大小,注意能小就小   矩阵从1开始   所以Matr 要+1   最大可以100
#define ll unsigned int
#define LL long long
struct mat//矩阵结构体,a表示内容,size大小 矩阵从1开始   但size不用加一
{
	ll a[Matr][Matr];
	mat()//构造函数
	{
		memset(a,0,sizeof(a));
	}
};
int Size ; 

mat multi(mat m1,mat m2)//两个相等矩阵的乘法,对于稀疏矩阵,有0处不用运算的优化
{
	mat ans=mat();
	for(int i=1;i<=Size;i++)
		for(int j=1;j<=Size;j++)
			if(m1.a[i][j])//稀疏矩阵优化
				for(int k=1;k<=Size;k++)
					ans.a[i][k]=(ans.a[i][k]+m1.a[i][j]*m2.a[j][k]); //i行k列第j项
	return ans;
}

mat quickmulti(mat m,LL n)//二分快速幂
{
	mat ans=mat();
	int i;
	for(i=1;i<=Size;i++)ans.a[i][i]=1;
	while(n)
	{
		if(n&1)ans=multi(m,ans);//奇乘偶子乘 挺好记的.
		m=multi(m,m);
		n>>=1;
	}
	return ans;
}

void print(mat m)//输出矩阵信息,debug用
{
	int i,j;
	printf("%d\n",Size);
	for(i=1;i<=Size;i++)
	{
		for(j=1;j<=Size;j++)
			printf("%u ",m.a[i][j]);
		printf("\n");
	}
}
int main()
{
	LL n;
	int t;
	int cas=1;
	int k;
	scanf("%d",&t);
	while(t--)
	{

		scanf("%lld%d",&n,&k);
		mat gouzao=mat(),chu=mat();

		printf("Case %d: ",cas++);
		Size=k+2;
		for(int i=1;i<=k+1;i++)
		{
			chu.a[1][i]=1;
		}
		for(int j=1;j<=k+1;j++)
		{
			for(int i=1;i<=j;i++)
			{
				if(i==1||i==j)
				{
					gouzao.a[i][j]=1;
					continue;
				}
				else
				{
					gouzao.a[i][j]=gouzao.a[i][j-1]+gouzao.a[i-1][j-1];

				}
			}
		}
		gouzao.a[k+1][k+2]=1;
		gouzao.a[k+2][k+2]=1;
		/*printf("chu\n");
		print(chu);
		printf("gouzao\n");
		print(gouzao);*/

		printf("%u\n",multi(chu,quickmulti(gouzao,n)).a[1][k+2]);
	}
	return 0;
} 

版权声明:本文为博主原创文章,未经博主允许不得转载。

时间: 2024-10-08 17:26:34

LightOJ 1132 - Summing up Powers 矩阵快速幂+排列组合的相关文章

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

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