题意

已知\(a_{i} = \sum_{j=1}^{i} \{^{i} _{j} \}b_{j}\), 给出\(a_{1} 到 a_{n}\) ;

求\(b_{l} 到 b_{r}\)在\(1e9+7\)的意义下取模的值;

\(1 \le l \le r \le n \le 10^5\)

\(r-l \le 100\)

\(0 \le a_{i} \lt 10^9 + 7\)

题解

  • Part1

  • 斯特林反演(https://www.cnblogs.com/hchhch233/p/10016543.html):

    • \(f(n) = \sum_{i=1}^{n} \{^n_i\} g(n) \Longleftrightarrow g(n) = \sum_{i=1}^{n} (-1)^{n-i} [^n_i]f(i)\)
  • 由于\(r-l\)比较小,可以直接暴力求,只需要快速求出\([^n_1] - [^n_n]\)

  • \([_m^n] = [x^m] x^{\overline n}\) , 其中\(x^{\overline n} = \Pi_{i=1}^n(x+i-1)\)

  • Part 2

  • 直接做是\(n log ^2n\)的,假设在计算\(F_n(x) = x^{\overline n}\)时已经计算了\(F_\frac{n}{2}(x)\),记\(m = \frac{n}{2} , F_{m} = \sum_{i=0}^{m}a_{i} x^{i}\)

  • 只需要计算\(F'_{m} (x) = F_{m}(x+m)\)和\(F_{m}\)相乘;

  • \[F'_m(x) = \sum_{i=0}^{m} a_{i} (x+m)^i \\ = \sum_{i=0}^{m} a _{i}\sum_{j=0}^{i}(^i_j)x^j m^{i-j} \\ = \sum_{i=0}^{m} \sum_{j=0}^{i} a_{i} \frac{i!}{j!(i-j)!}x^{j}m^{i-j} \\ = \sum_{i=0}^{m} \frac{x^i}{i!} \sum_{j=0}^{m-i} a_{i+j}(i+j)! \ \frac{m^{j}}{j!}
    \]

  • 记$A(x)=\sum_{i=0}^{m} a_{i} i! x^{i}, B(x) = \sum_{i=0}{m}\frac{m{i}}{i!} x^{m-i} $

  • 相乘之后取第\(m+i\)项可以算后面的东西

  • Part 3

  • 可是不是\(ntt\)模数,直接取模会爆\(long \ long\)需要用到拆系数\(ntt\);

  • 将\(A(x)*B(x)\)分解成\((k * A1(x) + A2(x)) * (k * B1(x) + B2(x)) , k 一般取\sqrt{n} (2^{15} ) 左右\)

  • 暴力是7次\(fft\),可以优化到需要4次\(fft\)

  • 两次\(fft\)合并成一次的做法,需要计算\(dft(A(x))和dft(B(x))\)

  • 记: $$\begin{align} P(x) = A(x) + iB(x)\ Q(x) = A(x) -iB(x) \end{align}$$

  • 考虑直接做\(P\)为\(P’\),考虑\(Q \ dft\)之后的结果\(Q'\)

    \[\begin{align} Q'[k] &= A(w_{n}^{k}) - iB(w_{n}^{k}) \\ &= \sum_{j=0}^{n-1} (a_{j}-ib_{j})w^{jk}_{n} \\ &= \sum_{j=0}^{n-1}(a_{j}-ib_{j})(cos(\frac{2\pi jk}{n} ) + isin(\frac{2\pi jk}{n}) ) \\ &= \sum_{j=0}^{n-1}(a_{j}cos(\frac{2\pi jk}{n}) + b_{j}sin(\frac{2\pi jk}{n}))- i(b_{j}cos(\frac{2\pi jk}{n}) - a_{j}sin(\frac{2\pi jk}{n})) \\ &= conj \sum_{j=0}^{n-1} (a_{j}cos(\frac{-2\pi jk}{n}) - b_{j}sin(\frac{-2\pi jk}{n} ) ) +i( b_{j}cos(\frac{-2\pi jk}{n}) + a_{j}sin(\frac{-2\pi jk}{n})) \\ &= conj \sum_{j=0}^{n-1} (a_{j}+ib_{j})(cos(\frac{-2\pi jk}{n})+isin(\frac{-2\pi jk}{n})) \\ &= conj \ \sum_{j=0}^{n-1} (a_{j}+ib_{j})w^{-jk}_{n} \\ &= conj \ P'[n-k] \end{align}
    \]

  • 所以由\(P'\)直接得到\(Q'\):

  • 即:$$ \begin{align} A'(x) = \frac{P'(x) + Q'(x)}{2} \ B'(x) = \frac{P'(x) - Q'(x)}{2i}\end{align} g $$

  • 结果也存在实部和虚部里直接\(idft\)回去就好了;

  • 注意精度;

    #include<bits/stdc++.h>
    #define ld long double
    #define ll long long
    using namespace std;
    const int N=400010,mod=1e9+7;
    const ld pi=acos(-1);
    int n,l,r,rev[N],fac[N],inv[N],p[N],a[N],b[N];
    char gc(){
    static char*p1,*p2,s[1000000];
    if(p1==p2)p2=(p1=s)+fread(s,1,1000000,stdin);
    return(p1==p2)?EOF:*p1++;
    }
    int rd(){
    int x=0;char c=gc();
    while(c<'0'||c>'9')c=gc();
    while(c>='0'&&c<='9')x=(x<<1)+(x<<3)+c-'0',c=gc();
    return x;
    }
    int pw(int x,int y){
    int re=1;
    while(y){
    if(y&1)re=(ll)re*x%mod;
    y>>=1,x=(ll)x*x%mod;
    }
    return re;
    }
    struct C{
    ld x,y;
    C(ld _x=0,ld _y=0):x(_x),y(_y){};
    C operator +(const C&A)const{return C(x+A.x,y+A.y);}
    C operator -(const C&A)const{return C(x-A.x,y-A.y);}
    C operator *(const C&A)const{return C(x*A.x-y*A.y,x*A.y+y*A.x);}
    C operator /(const ld&A)const{return C(x/A,y/A);}
    C operator !()const{return C(x,-y);}
    };
    void fft(C*A,int len,int f){
    for(int i=0;i<len;++i)if(i<rev[i])swap(A[i],A[rev[i]]);
    for(int i=1;i<len;i<<=1){
    C wn=C(cos(pi/i),f*sin(pi/i));
    for(int j=0;j<len;j+=i<<1){
    C w=C(1,0);
    for(int k=0;k<i;++k,w=w*wn){
    C x=A[j+k],y=w*A[j+k+i];
    A[j+k]=x+y,A[j+k+i]=x-y;
    }
    }
    }
    if(!~f){for(int i=0;i<len;++i){A[i]=A[i]/len;}}
    }
    void mul(int*A,int*B,int n){
    static C t1[N],t2[N],t3[N],t4[N];
    int L,len;
    for(L=0,len=1;len<(n+1)<<1;++L,len<<=1);
    for(int i=0;i<len;++i){rev[i]=(rev[i>>1]>>1)|((i&1)<<(L-1));}
    for(int i=0;i<=n;++i)t1[i]=C(A[i]>>15,A[i]&0x7fff),t2[i]=C(B[i]>>15,B[i]&0x7fff);
    for(int i=n+1;i<len;++i)t1[i]=t2[i]=C(0,0);
    fft(t1,len,1),fft(t2,len,1);
    for(int i=0;i<len;++i){
    C x1=t1[i],y1=!t1[(len-i)&(len-1)];
    C x2=t2[i],y2=!t2[(len-i)&(len-1)];
    t3[i]=(x1+y1)*x2*C(0.5,0);
    t4[i]=(x1-y1)*x2*C(0,-0.5);
    // t3[i]=((x1+y1)*(x2+y2)+(x1+y1)*(x2-y2))*C(0.25,0);
    // t4[i]=((x1-y1)*(x2+y2)+(x1-y1)*(x2-y2))*C(0,-0.25);
    }
    fft(t3,len,-1),fft(t4,len,-1);
    for(int i=0;i<=n<<1;++i){
    A[i]=((((ll)(t3[i].x+0.1)%mod)<<30)+(((ll)(t3[i].y+t4[i].x+0.1)%mod)<<15)+(ll)(t4[i].y+0.1))%mod;
    }
    }
    void solve(int*A,int n){
    if(n==0){A[0]=1;return;}
    if(n==1){A[1]=1;return;}
    if(n&1){
    solve(A,n-1);
    for(int i=n;i;--i)A[i]=(A[i-1]+(ll)A[i]*(n-1))%mod;
    A[0]=0;
    return;
    }
    static int t1[N],t2[N];
    int m=n>>1;
    solve(A,m);
    for(int i=0;i<=m;++i)t1[i]=(ll)A[i]*fac[i]%mod;
    for(int i=0,t=1;i<=m;++i,t=(ll)t*m%mod)t2[m-i]=(ll)t*inv[i]%mod;
    mul(t1,t2,m);
    for(int i=0;i<=m;++i){t1[i]=(ll)inv[i]*t1[m+i]%mod;}
    mul(A,t1,m);
    }
    int main(){
    freopen("gift.in","r",stdin);
    freopen("gift.out","w",stdout);
    n=rd();l=rd();r=rd();
    for(int i=1;i<=n;++i)p[i]=rd();
    for(int i=fac[0]=inv[0]=1;i<=n;++i){
    fac[i]=(ll)fac[i-1]*i%mod;
    inv[i]=pw(fac[i],mod-2);
    }
    solve(a,l-1);
    for(int i=l-1;i<=r;++i){
    for(int j=1;j<=i;++j){
    int t=(ll)a[j]*p[j]%mod;
    if((i-j)&1)b[i]=(b[i]-t+mod)%mod;
    else b[i]=(b[i]+t)%mod;
    }
    for(int j=i+1;j;--j){a[j]=(a[j-1]+(ll)a[j]*i)%mod;}a[0]=0;
    // b[i]=(b[i]-b[i-1]+mod)%mod;
    }
    for(int i=l;i<=r;++i)printf("%d ",(b[i]-b[i-1]+mod)%mod);
    return 0;
    }

【纪中集训2019.3.12】Z的礼物的更多相关文章

  1. 【纪中集训2019.3.12】Mas的仙人掌

    题意: ​ 给出一棵\(n\)个点的树,需要加\(m\)条边,每条边脱落的概率为\(p_{i}\) ,求加入的边在最后形成图中仅在一个简单环上的边数的期望: \(1 \le n \ , m \le 1 ...

  2. 【纪中集训2019.3.27】【集训队互测2018】小A的旅行(白)

    题目 描述 ​ \(0-n-1\)的图,满足\(n\)是\(2\)的整数次幂, $ i \to j $ 有 $ A_{i,j} $ 条路径: ​ 一条路径的愉悦值定义为起点和终点编号的\(and\)值 ...

  3. 【纪中集训2019.3.23】Deadline

    题意 描述 一个二分图\((A,B)\),每个点额外有一个颜色0或者1: 匹配时,只能相同颜色的点匹配: 给出\(A\)中的颜色,问如何分配\(B\)种的颜色使得\((A,B)\)的最大匹配最小: 范 ...

  4. 【纪中集训2019.3.23】IOer

    题目 描述 你要在\(m\)天内,刷\(n\)道题,每天可以刷的题的数目不限: 第\(i\)天可以刷的题目的种类是\(ui+v\): 两种刷题的方案不同当且仅当某天刷题的数量不同或者依次刷题的种类不同 ...

  5. 【纪中集训2019.3.11】Cubelia

    题目: 描述 给出长度为\(n\)的数组\(a\)和\(q\)个询问\(l,r\). 求区间\([l,r]\)的所有子区间的前缀和的最大值之和: 范围: $n \le 2 \times 10^5 , ...

  6. 【纪中集训2019.3.13】fft

    题意: 描述 一共有\(n+m\)道题,其中\(n\)道答案是\(A\),\(m\)道答案是\(B\): 你事先知道\(n和m\),问在最优情况下的期望答错次数,对\(998244353\)取模: 范 ...

  7. 「中山纪中集训省选组D1T1」最大收益 贪心

    题目描述 给出\(N\)件单位时间任务,对于第\(i\)件任务,如果要完成该任务,需要占用\([S_i, T_i]\)间的某个时刻,且完成后会有\(V_i\)的收益.求最大收益. 澄清:一个时刻只能做 ...

  8. 纪中集训 Day 2

    今天(其实是昨天= =)早上起来发现好冷好冷啊= = 吃完饭就准备比赛了,好吧B组难度的题总有一道不知到怎么写QAQ 太弱了啊!!! 蒟蒻没人权啊QAQ 今天第4题不会写,在这里说说吧 题目的意思就是 ...

  9. 纪中集训 Day1

    今天早上起来吃饭,发现纪中伙食真的是太差了!!!什么都不热,早餐的面包还好,然后就迎来了美好的早晨= = 早上做一套题,T1T2果断秒,T3一看就是noi原题,还好看过题解会写,然后就愉快的码+Deb ...

随机推荐

  1. 基于神念TGAM的脑波小车(2)

    将数据处理移植到STM32上,采用串口的DMA接收模式,注意的是DMA_MODE采用Circular,DMA_BufferSize>(8*512+36=4132)(小包8个字节,每秒512个,完 ...

  2. 华为笔试——C++最高分问题

    题目介绍:现在输入一组数据,写入学生的考试分数.已知学生数为N,学生编号为1到N,且0<N<=30000,每个学生都有一个分数:操作数为M且0<M<5000.输入第一行为N M ...

  3. nginx 在ubuntu上使用笔记(绑定域名)

    1. 重启nginx的两个语句: sudo service nginx restart sudo nginx -s reload 2. nginx配置文件路径: etc/nginx/ 尤其是 site ...

  4. route命令详情

    基础命令学习目录首页 原文链接:https://www.cnblogs.com/lpfuture/p/5857738.html 考试题一:linux下如何添加路由(百度面试题) 以上是原题,老男孩老师 ...

  5. $.each()用法

    通过它,你可以遍历对象.数组的属性值并进行处理. 使用说明 each函数根据参数的类型实现的效果不完全一致: 1.遍历对象(有附加参数) $.each(Object, function(p1, p2) ...

  6. 实验二 Java面向对象程序设计 20135321

    课程:Java程序设计   班级:1353    姓名:余佳源  学号:20135321 成绩:             指导教师:娄嘉鹏      实验日期:2015-5-8 实验密级:       ...

  7. mianshi

    https://blog.csdn.net/u012557610/article/details/80350099 https://blog.csdn.net/liuqiyao_01/article/ ...

  8. Beat(2/7)

    目录 摘要 团队部分 个人部分 摘要 队名:小白吃 组长博客:hjj 作业博客:beta冲刺(2/7) 团队部分 后敬甲(组长) 过去两天完成了哪些任务 整理博客 做了点商家数据表格 接下来的计划 做 ...

  9. BNUOJ 52303 Floyd-Warshall Lca+bfs最短路

    题目链接: https://www.bnuoj.com/v3/problem_show.php?pid=52303 Floyd-Warshall Time Limit: 60000msMemory L ...

  10. BETA-2

    前言 我们居然又冲刺了·二 团队代码管理github 站立会议 队名:PMS 530雨勤(组长) 过去两天完成了哪些任务 了解OpenCV下的视频参数及其调用方法 初步编码 接下来的计划 文档工作 速 ...