Daxia在2016年5月期間去瑞士度蜜月,順便拜訪了Wzc,Wzc給他出了一個問題：

Wzc給Daxia等差數列A(0),告訴Daxia首項a和公差d;

首先讓Daxia求出數列A(0)前n項和,得到新數列A(1);

然后讓Daxia求出數列A(1)前n項和,得到新數列A(2);

接著讓Daxia求出數列A(2)前n項和,得到新數列A(3);

?

規律題，首先寫出 a、a+d、a+2d、a+3d...這個容易寫出

下面一行也容易寫出：a、2a+d、3a+3d....

再下一行，確實難寫，但是通過上面兩行可以得出，dp[i][j] = dp[i-1][j] + dp[i][j-1]

然后可以順利寫出后面的，找到通項公式。

ans = C(m+i-1,m)*a + C(m+i-1,i-2)*d

?

#include <cstdio> #include <cstdlib> #include <cstring> #include <cmath> #include <algorithm> using namespace std; #define inf (0x3f3f3f3f) typedef long long int LL;#include <iostream> #include <sstream> #include <vector> #include <set> #include <map> #include <queue> #include <string> LL a,d,m,i; const int MOD = 1000000007; LL quick_pow (LL a,LL b,LL MOD) {//求解 a^b%MOD的值LL base=a%MOD;LL ans=1; //相乘，所以這里是1while (b){if (b&1){ans=(ans*base)%MOD; //如果這里是很大的數據，就要用quick_mul }base=(base*base)%MOD; //notice//注意這里,每次的base是自己base倍b>>=1;}return ans; } LL C (LL n,LL m,LL MOD) {if (n<m) return 0; //防止sb地在循環，在lucas的時候if (n==m) return 1;LL ans1 = 1;LL ans2 = 1;LL mx=max(n-m,m); //這個也是必要的。能約就約最大的那個LL mi=n-mx;for (int i = 1; i <= mi; ++i){ans1 = ans1*(mx+i)%MOD;ans2 = ans2*i%MOD;}return (ans1*quick_pow(ans2,MOD-2,MOD)%MOD); //這里放到最后進行,不然會很慢 }void work () {if (i==1){printf ("%lld\n",a);return;}LL ans = (C(m+i-1,m,MOD)*a%MOD+C(m+i-1,i-2,MOD)*d%MOD)%MOD;printf ("%lld\n",ans);return ; } int main() { #ifdef localfreopen("data.txt","r",stdin); #endifwhile (scanf("%lld%lld%lld%lld",&a,&d,&m,&i)!=EOF) work();return 0; } View Code

?

轉載于:https://www.cnblogs.com/liuweimingcprogram/p/5785258.html

總結

以上是生活随笔為你收集整理的FZU Problem 2238 Daxia Wzc's problem的全部內容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯，歡迎將生活随笔推薦給好友。