洛谷noip 模拟赛 day1 T3

最新推荐文章于 2021-09-25 14:41:22 发布

weixin_33953249

最新推荐文章于 2021-09-25 14:41:22 发布

阅读量123

点赞数

原文链接：http://www.cnblogs.com/lyzuikeai/p/7263132.html

版权

T7983 大芳的逆行板载

题目背景

大芳有一个不太好的习惯：在车里养青蛙。青蛙在一个n厘米（11n毫米s）的Van♂杆子上跳来跳去。她时常盯着青蛙看，以至于突然逆行不得不开始躲交叉弹。有一天他突发奇想，在杆子上每1厘米为一个单位，瞎涂上了墨水，并且使用mOgic，使青蛙跳过之处墨水浓度增加x。当然，他还会闲着无聊滴几滴墨水再涂♂抹均匀。

他现在无时无刻都想知道，第l厘米到第r厘米墨水的浓度是多少？

哦不!等等，他现在找到了一个计算器，可以输入几个数字与x，计算他们的x次幂和，所以。。。他想知道的是第l厘米到第r厘米墨水的浓度的x次幂和是多少？

题目描述

大芳有3种~~舰长技能~~骚操作

续：把青蛙放到第l厘米处，戳青蛙使其跳至r。效果：第l厘米至第r厘米墨水浓度增加x
抚♂摸：擦干杆子某一部分，重新滴加墨水并抹匀。效果：使第l厘米至第r厘米墨水浓度都变成x

最后一种是：

压线逆行，将车流看做⑨弹幕找安定点，掏出计算器，大喊板载后计算：

第l厘米至第r厘米墨水浓度的x次幂和是几何？记得答案要

模

输入输出格式

输入格式：

第一行

第二行

接下来每行4个数字，依次为：操作编号（1、2或3），

输出格式：

每次进行3操作，输出一行表示答案

记得膜模1000000007

输入输出样例

输入样例#1：

5 5
19844 14611 26475 4488 6967 
2 1 3 15627
2 1 2 30113
2 3 5 14686
2 5 5 32623
3 1 2 8

输出样例#1：

466266421

说明

对于20%的数据，满足 $n,m≤1000n,m\leq 1000n,m≤1000$

对于另外20%的数据，满足 $k≤1k\leq 1k≤1$

对于另外20%的数据，满足 $k≤2k\leq 2k≤2$

对于另外20%的数据，满足 $n,m≤50000n,m\leq 50000n,m≤50000$

对于100%的数据，满足 $n,m≤100000,0≤k≤10n,m\leq 100000,0\leq k \leq 10n,m≤100000,0≤k≤10$

操作1，2对应的 $x≤109+7x\le 10^9+7x≤109+7$

$(x+c)^n=sigma i=0..n (c^i)*C(n,i)*(x^(n-i))$

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#define LL long long
using namespace std;
const int M=1<<18,mod=1000000007;
LL read(){
    LL ans=0,f=1,c=getchar();
    while(c<'0'||c>'9'){if(c=='-') f=-1; c=getchar();}
    while(c>='0'&&c<='9'){ans=ans*10+(c-'0'); c=getchar();}
    return ans*f;
}
LL C[15][15],c[15],q[15],num[M],w; 
int h,n,m,L,R;
struct node{
    LL f,tag;
    LL m[11];
}tr[M];
void calc_f(int x,int l,int r,LL v){
    LL sz=r-l+1;
    tr[x].m[1]=sz*v%mod;
    for(int i=2;i<=10;i++) tr[x].m[i]=(tr[x].m[i-1]*v)%mod;
}
void push(int x,int k,LL v){
    LL sum=1;
    for(int i=1;i<=k;i++) sum=sum*v%mod,tr[x].m[k]=(tr[x].m[k]+sum*C[k][i]%mod*tr[x].m[k-i]%mod)%mod;
}
void calc_tag(int x,LL v){
    for(int i=10;i>=1;i--) 
    push(x,i,v);
}
void up(int x){
    int l=x<<1,r=x<<1^1;
    for(int i=0;i<=10;i++) tr[x].m[i]=(tr[l].m[i]+tr[r].m[i])%mod;
}
void down(int x,int l,int r){
    int ll=x<<1,rr=x<<1^1;
    LL mid=(l+r)>>1;
    if(tr[x].f!=-1){
        tr[ll].f=tr[x].f; tr[rr].f=tr[x].f;
        calc_f(ll,l,mid,tr[x].f);
        calc_f(rr,mid+1,r,tr[x].f);
        tr[x].f=-1;    tr[ll].tag=tr[rr].tag=0;
    }
    if(tr[x].tag){
        if(l==r) return ;
        calc_tag(ll,tr[x].tag);
        calc_tag(rr,tr[x].tag);
        tr[x<<1].tag+=tr[x].tag; tr[x<<1^1].tag+=tr[x].tag; tr[x].tag=0;
    }
}
void modify_tag(int x,int l,int r){
    if(L<=l&&r<=R){
        tr[x].tag+=w;
        if(l==r){
            tr[x].m[1]+=w;
            for(int i=2;i<=10;i++) tr[x].m[i]=tr[x].m[i-1]*tr[x].m[1]%mod;
        }
        else calc_tag(x,w);
        return ;
    }
    down(x,l,r);
    int mid=(l+r)>>1;
    if(L<=mid) modify_tag(x<<1,l,mid);
    if(R>mid) modify_tag(x<<1^1,mid+1,r);
    up(x);
}
void modify_f(int x,int l,int r){
    if(L<=l&&r<=R){
        tr[x].f=w;
        tr[x].tag=0;
        calc_f(x,l,r,w);
        return ;
    }
    down(x,l,r);
    int mid=(l+r)>>1;
    if(L<=mid) modify_f(x<<1,l,mid);
    if(R>mid) modify_f(x<<1^1,mid+1,r);
    up(x);
}
LL query(int x,int l,int r,int k){
    if(L<=l&&r<=R) return tr[x].m[k];
    down(x,l,r);
    LL mid=(l+r)>>1;
    LL sum=0;
    if(L<=mid) sum=sum+query(x<<1,l,mid,k);
    if(R>mid) sum=sum+query(x<<1^1,mid+1,r,k);
    return sum%mod;
}
void build(int x,int l,int r){
    tr[x].f=-1;
    if(l==r){
        tr[x].m[0]=1;
        tr[x].m[1]=num[l]; 
        for(int i=2;i<=10;i++) tr[x].m[i]=(tr[x].m[i-1]*tr[x].m[1])%mod; 
        return ;
    }
    int mid=(l+r)>>1;
    build(x<<1,l,mid);
    build(x<<1^1,mid+1,r);
    up(x);
}
void prepare(){
    q[0]=1; q[1]=1;
    for(int i=2;i<=10;i++) q[i]=q[i-1]*i;
    for(int i=1;i<=10;i++)
     for(int j=1;j<=i;j++)
      C[i][j]=q[i]/(q[j]*q[i-j]);
}
int main()
{
    prepare();
    n=read(); m=read();
    for(int i=1;i<=n;i++) num[i]=read();
    build(1,1,n); 
    for(int i=1;i<=m;i++){
        h=read();
        if(h==1){
            L=read(); R=read(); w=read();
            modify_tag(1,1,n);
        }
        else if(h==2){
            L=read(); R=read(); w=read();
            modify_f(1,1,n);
        }
        else{
            L=read(); R=read(); 
            int k=read();
            printf("%lld\n",query(1,1,n,k));
        }//printf("[%lld]\n",tr[1].m[1]);
    }
    return 0;
}