题目

给出一个数组A,经过一次处理,生成一个数组S,数组S中的每个值相当于数组A的累加,比如:A = {1 3 5 6} => S = {1 4 9 15}。如果对生成的数组S再进行一次累加操作,{1 4 9 15} => {1 5 14 29},现在给出数组A,问进行K次操作后的结果。(输出结果 Mod 10^9 + 7)

分析

发现,每次处理相当于将A卷上一个\(I(\forall a_i=1)\)

于是机智的我在wiki又发现狄利克雷卷积满足交换律(我居然才知道)

于是快速幂咯,时间复杂度\(O(nlog^2n)\),常数巨大,在51nod的老爷机上根本过不了。

然后就TLE得一塌糊涂。

于是找规律,发现\(ans_k=\sum_{i+j=k}A_iC_{j+n-1}^{j}\)

然后NTT

但是,\(10^9+7\)并没有原根,所以祭出三模数NTT(如果有人想用高精度,我也没办法)。

因为\((10^9+7)^2*n\approx10^{23}\),所以找三个\(10^9\)左右的模数。

假设

\[ans≡a_0(mod\ m_0)
\]

\[ans≡a_1(mod\ m_1)
\]

\[ans≡a_2(mod\ m_2)
\]

因为\(m_0m_1m_2\)会爆long long

所有,通过CRT(中国剩余定理)合并前两项,于是\(M=m_0*m_1\)

\[ans≡A(mod\ M)
\]

\[ans≡a_2(mod\ m_2)
\]

设\(ans=kM+A\)

因为$$kM+A≡ans≡a_2(mod\ m_2)$$

所以$$k≡(a_2-A)M^{-1}(mod\ m_2)$$

那么根据上面的式子求出k,通过\(kM+A=ans\)就可以求出ans了。

#include <cmath>
#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <algorithm>
#include <queue>
#include <map>
#include <bitset>
#include <set>
const int maxlongint=2147483647;
const long long mo=1e9+7;
const int N=200005;
using namespace std;
long long mod[3]={469762049,998244353,1004535809},M=mod[0]*mod[1],W[3][N];
long long f[3][N],g[3][N],v[3][N],v_[N],jc[3][N],ny[3][N],inv[N];
int n,m,fn;
long long poww(long long x,long long y,int z)
{
long long s=1;
x%=mod[z];
for(;y;x=x*x%mod[z],y>>=1)
if(y&1) s=s*x%mod[z];
return s;
}
long long mul(long long x,long long y,long long mo)
{
long long s=0;
x%=mo;
for(;y;x=(x+x)%mo,y>>=1)
if(y&1) s=(s+x)%mo;
return s;
}
void NTT(long long *f,int type,int z)
{
for(int i=0,p=0;i<fn;i++)
{
if(i<p) swap(f[i],f[p]);
for(int j=fn>>1;(p^=j)<j;j>>=1);
}
for(int i=2;i<=fn;i<<=1)
{
int half=i>>1,pe=fn/i;
for(int j=0;j<half;j++)
{
long long w=!type?W[z][j*pe]:W[z][fn-j*pe];
for(int k=j;k<fn;k+=i)
{
long long x=f[k],y=f[k+half]*w%mod[z];
f[k]=(x+y)%mod[z],f[k+half]=(x-y+mod[z])%mod[z];
}
}
}
if(type)
{
long long ni=poww(fn,mod[z]-2,z);
for(int i=0;i<fn;i++) f[i]=f[i]*ni%mod[z];
}
}
long long CRT(long long a0,long long a1,long long a2)
{
long long n0=poww(mod[1],mod[0]-2,0),n1=poww(mod[0],mod[1]-2,1);
long long A=(mul(a0*mod[1]%M,n0%M,M)+mul(a1*mod[0]%M,n1%M,M))%M,n2=poww(M%mod[2],mod[2]-2,2);
long long k=(a2-A)%mod[2]*n2%mod[2];
k=(k%mod[2]+mod[2])%mod[2];
return (k%mo*(M%mo)%mo+A)%mo;
}
int main()
{
scanf("%d%d",&n,&m);
for(fn=1;fn<=n*2+2;fn<<=1);
inv[0]=inv[1]=1;
for(int i=2;i<=fn;i++) inv[i]=(-(mo/i)*inv[mo%i]%mo+mo)%mo;
for(int i=0;i<n;i++) scanf("%lld",&f[0][i]),f[1][i]=f[2][i]=f[0][i];
for(int j=0;j<3;j++)
{
W[j][0]=1,W[j][1]=poww(3,(mod[j]-1)/fn,j);
for(int i=2;i<=fn;i++) W[j][i]=W[j][i-1]*W[j][1]%mod[j];
for(int i=0;i<n;i++) v[j][i]=1;
g[j][0]=1;
for(int i=1;i<n;i++) g[j][i]=g[j][i-1]*(i+m-1)%mo*inv[i]%mo;
}
for(int j=0;j<3;j++)
{
NTT(f[j],0,j),NTT(g[j],0,j);
for(int i=0;i<fn;i++) f[j][i]=f[j][i]*g[j][i]%mod[j];
NTT(f[j],1,j);
}
for(int i=0;i<n;i++) printf("%lld\n",CRT(f[0][i],f[1][i],f[2][i]));
}

51nod 1172 Partial Sums V2的更多相关文章

  1. 51nod 1172 Partial Sums V2 卡精度的任意模数FFT

    卡精度的任意模数fft模板题……这道题随便写个表就能看出规律来(或者说考虑一下实际意义),反正拿到这题之后,很快就会发现他是任意模数fft模板题.然后我就去网上抄了一下板子……我打的是最土的任意模数f ...

  2. 51nod1172 Partial Sums V2

    推一下式子发现是裸的FFT,$ans[k]=\sum_{i}\sum_{j}[i+j=k]a[i]*C_{m-1+j}^{j}$ 比较坑爹的就是这个模数,于是我们上任意模数的FFT 任意模数的FFT目 ...

  3. 51nod 1161 Partial Sums

    给出一个数组A,经过一次处理,生成一个数组S,数组S中的每个值相当于数组A的累加,比如:A = {1 3 5 6} => S = {1 4 9 15}.如果对生成的数组S再进行一次累加操作,{1 ...

  4. 51nod1161 Partial Sums

    开始想的是O(n2logk)的算法但是显然会tle.看了解题报告然后就打表找起规律来.嘛是组合数嘛.时间复杂度是O(nlogn+n2)的 #include<cstdio> #include ...

  5. Non-negative Partial Sums(单调队列)

    Non-negative Partial Sums Time Limit: 6000/3000 MS (Java/Others)    Memory Limit: 65536/32768 K (Jav ...

  6. hdu 4193 Non-negative Partial Sums 单调队列。

    Non-negative Partial Sums Time Limit: 6000/3000 MS (Java/Others)    Memory Limit: 65536/32768 K (Jav ...

  7. TOJ 1721 Partial Sums

    Description Given a series of n numbers a1, a2, ..., an, the partial sum of the numbers is defined a ...

  8. 【计数】cf223C. Partial Sums

    考试时候遇到这种题只会找规律 You've got an array a, consisting of n integers. The array elements are indexed from ...

  9. CodeForces 223C Partial Sums 多次前缀和

    Partial Sums 题解: 一个数列多次前缀和之后, 对于第i个数来说他的答案就是 ; i <= n; ++i){ ; j <= i; ++j){ b[i] = (b[i] + 1l ...

随机推荐

  1. Hadoop部署(伪分布式系统)

    hadoop安装 #修改主机名 hostnamectl set-hostname hadoop #修改hosts vim /etc/hosts #追加到末尾 10.0.0.11 hadoop 安装必备 ...

  2. PAT B1006 换个格式输出整数 (15)

    AC代码 #include <cstdio> const int max_n = 3; char radix[max_n] = {' ', 'S', 'B'}; int ans[max_n ...

  3. Mysql 三大特性详解

    Mysql 三大特性详解 Mysql Innodb后台线程 工作方式 首先Mysql进程模型是单进程多线程的.所以我们通过ps查找mysqld进程是只有一个. 体系架构 InnoDB存储引擎的架构如下 ...

  4. X86逆向1:软件破解入门课【课件下载】

    从本节课开始,我将带领小白入门学习软件破解的相关内容,大佬绕过,以后将会定期更新从最基本的破解知识点开始学习,由简单到复杂循序渐进,难度会逐步提高. 为了防止版权方面的争议,我将自行编写一些破解案例来 ...

  5. Java反射的理解(六)-- 通过反射了解集合泛型的本质

    Java反射的理解(六)-- 通过反射了解集合泛型的本质 上述写了那么多,我们可能会有个疑问,为什么要用反射,步骤比我们常规的加载类操作复杂多了,别急,这个问题我最后才解答,我们先来了解集合泛型的本质 ...

  6. js ajax return false了,仍然会往下执行

    function checkMust(){ var flag=false; $.getJSON("../Ajax/Carton/Ajax_TMSOrder_Create.ashx?r=&qu ...

  7. 搭建自己的框架WedeNet(五)

    WedeNet2018.WedeWcfServices-WCF服务层:结构如下: 就是定义了服务契约接口和服务类,以OrderServices为例,如下: using System; using Sy ...

  8. 安装多个ORACLE导致多个Oracle HOME的情况!

    Oracle由于版本的不同,在注册表中产生的注册表信息也有所不同,但主要的键值信息还是一样的,例如Oracle10g比oracle9i在注册表中表现的更为“简洁”,在未知的情况下,获取Oracle10 ...

  9. 谷歌浏览器禁用JS步骤

    前奏: 当你想访问一个页面,因为不符合访问条件,而被JS阻拦:或者你打开的页面特效太多,干扰视线:亦或者JS的某个功能禁止了你的某些行为: 是不是很难受, 好办,咱给他禁了不就好啦,大胆地屏蔽它,不管 ...

  10. vue-$watch属性方法

    特性 https://www.cnblogs.com/widgetbox/p/8954162.html https://segmentfault.com/a/1190000012948175?utm_ ...