Codeforces 223APartial Sums 数论+组合数学

标签：os   io   art   for   re   c   

题意很简单，求不是那么好求的，k很大 要操作很多次，所以不可能直接来的，印象中解决操作比较多无非线段树 循环节 矩阵 组合数等等吧，这道题目 也就只能多画画什么 的了 

就以第一个案例为主吧 ，

3 

1 2 3

k我们依据画的次数来自己定好了

下面的每个数表示这个位置的 数由最初的 数组num[]中多少个数加起来得到的

当k为0的时候呢，就是

1 1 1

k为1的时候呢

1 2 3

k为2的时候呢

1 3 6

那么k为3的时候

1 4 10

这里看一下 从数组下标0开始，那么其实就是 C(i + k,i)

觉得不够的话呢 可以再多写几个，发现就是这么回事啊 ，跟组合数有联系了，那么肯定不可能每一次都求啊 ，所以可以试着搞一个矩阵，这样就可以做了

做的时候组合数直接来超时了，可以先用一个c数组预处理出所有的组合数答案，求答案运用C(n,m) == C(n,n - m)这样省时间这样也是700+ms比较慢，当然若是再把乘法逆元给 先预处理存到数组里会更加省时间 而且会省很多，

#include<iostream>
#include<cstdio>
#include<list>
#include<algorithm>
#include<cstring>
#include<string>
#include<queue>
#include<stack>
#include<map>
#include<vector>
#include<cmath>
#include<memory.h>
#include<set>

#define ll long long

#define eps 1e-8

const int inf = 0xfffffff;

const ll INF = 1ll<<61;

using namespace std;

//vector<pair<int,int> > G;
//typedef pair<int,int > P;
//vector<pair<int,int> > ::iterator iter;
//
//map<ll,int >mp;
//map<ll,int >::iterator p;

#define MOD 1000000007

ll num[2000 + 5];

int n,k;

ll exgcd(ll a, ll b, ll &x, ll &y)  {
	if(!b) {
		x = 1; y = 0;
		return a;
	}
	ll r = exgcd(b, a%b, y, x);
	y -= a/b*x;
	return r;
}

ll inv(ll a, ll m)
{
	ll x,y,gcd = exgcd(a, m, x, y);
	if(x < 0)  
		x += m;
	return x;
}

ll C(ll n,ll m) {
	ll ans = 1;
	for(int i=1;i<=m;i++)
		ans = ((ans * inv(i,MOD))%MOD * (n - i + 1))%MOD;
	return ans;
}

ll c[2000 + 5];

int main() {
	while(scanf("%d %d",&n,&k) == 2) {
		for(int i=0;i<n;i++)
			scanf("%d",&num[i]);
		if(k == 0) {
			for(int i=0;i<n;i++)
				printf("%I64d%c",num[i],i == n - 1?'\n':' ');
			continue;
		}
		for(int i=0;i<n;i++)
			c[i] = C(i + k - 1,i);
		for(int i=0;i<n;i++) {
			ll ans = 0ll;
			for(int j=0;j<=i;j++) {
				//ll tmp = C(i - j + k - 1,i - j);
				//ll tt = num[j];
				ans = (ans + num[j] * c[i - j])%MOD;
			}
			printf("%I64d%c",ans,i == n - 1?'\n':' ');
		}
	}
	return 0;
}

Codeforces 223APartial Sums 数论+组合数学,布布扣,bubuko.com

Codeforces 223APartial Sums 数论+组合数学

标签：os   io   art   for   re   c   

原文地址：http://blog.csdn.net/yitiaodacaidog/article/details/38071257