CodeForces 266E More Queries to Array...(线段树+式子展开)

原创 2016年08月30日 21:43:14

开始觉得是规律题的,自以为是的推了一个规律,结果测试数据都没过。。。。看了love神的博客才发现只是把式子展开就找到规律了。不过挺6的是我虽然想错了,但是维护的的东西没有错,只是改改(改了进两个小时好吗????)就过了
题意是说给你一串数字,然后两种操作:
“= l r x” 是把数组第l位置到r位置的值都变为x
“? l r k”就是求一个公式的和 ∑(l <= i <= r) ai*(i-l+1)^k

因为k比较小嘛,可以展开式子看一下规律,展开后注意要往 构造常数或者构造递归 迭代 方面想。
当k=0时为 ai
当k=1时为 ai*(1-l) + ai*i
当k=2时为 ai*(1-l)^2+ai*2*i*(1-l)+ai*i^2
这样就可以得到:当k>0时,每个k都可以使用之前的(0 , k-1)迭代求出来
其中的系数就是二项式的系数(话说我更喜欢用杨辉三角来算),这样就直接线段树每个节点只维护ai* i^j(0 <= j <= 5)就好了。注意这儿是区间更新,所以要预处理一个对于每个k次方的前缀和,这样lazy操作是就可以节约时间了,还有就是注意取模之前要看看是否是负数,因为负数的话要+mod

#include<set>
#include<map>
#include<queue>
#include<stack>
#include<cmath>
#include<vector>
#include<string>
#include<cstdio>
#include<cstring>
#include<stdlib.h>
#include<iostream>
#include<algorithm>
using namespace std;
#define eps 1E-8
/*注意可能会有输出-0.000*/
#define Sgn(x) (x<-eps? -1 :x<eps? 0:1)//x为两个浮点数差的比较,注意返回整型
#define Cvs(x) (x > 0.0 ? x+eps : x-eps)//浮点数转化
#define zero(x) (((x)>0?(x):-(x))<eps)//判断是否等于0
#define mul(a,b) (a<<b)
#define dir(a,b) (a>>b)
typedef long long ll;
typedef unsigned long long ull;
const int Inf=1<<30;
const double Pi=acos(-1.0);
const ll Mod=1000000007ll;
const int Max=100010<<2;
int dir[6][6]={{0},{1,1},{1,2,1},{1,3,3,1},{1,4,6,4,1},{1,5,10,10,5,1}};//符合二项式
ll ans;//注意线段树维护的是权值与下标的k次方的乘积
ll segtr[6][Max],mark[Max];//存6个值,代表k从0到5 区间标记
char str[5];
ll ini[100005][6];//预处理
void Init()//预处理出一个次方前缀和
{
    for(int j=0;j<6;j++)
        ini[0][j]=0ll;
    for(int i=1;i<100001;i++)//底数为i
    {
        ini[i][0]=1ll;
        for(int j=1;j<6;j++)//次方为j
        {
            ini[i][j]=(ini[i][j-1]*i)%Mod;
        }
    }
    for(int i=1;i<100001;i++)
        for(int j=0;j<6;j++)
        ini[i][j]=(ini[i][j]+ini[i-1][j])%Mod;//前缀和,注意分开求
    return;
}
void Upnow(int now,int next)
{
    for(int i=0;i<6;i++)
        segtr[i][now]=(segtr[i][next]+segtr[i][next|1])%Mod;
    return ;
}
void Create(int sta,int enn,int now)
{
    mark[now]=-1ll;
    if(sta==enn)
    {
        scanf("%I64d",&segtr[0][now]);
        for(int i=1;i<6;i++)
        segtr[i][now]=(segtr[i-1][now]*sta)%Mod;//存下每个数的位置的k次方
        return;
    }
    int mid=dir(sta+enn,1);
    int next=mul(now,1);
    Create(sta,mid,next);
    Create(mid+1,enn,next|1);
    Upnow(now,next);
    return;
}
void Downnow(int now,int mid,int next,int sta,int enn)//lazy操作
{
    if(mark[now]>-1ll)
    {
        mark[next]=mark[now];
        mark[next|1]=mark[now];
        for(int i=0;i<6;i++)
        {
            segtr[i][next]=((mark[now]*(ini[mid][i]-ini[sta-1][i]))%Mod+Mod)%Mod;
            segtr[i][next|1]=((mark[now]*(ini[enn][i]-ini[mid][i]))%Mod+Mod)%Mod;
        }
        mark[now]=-1ll;
    }
    return;
}
void Update(int sta,int enn,int now,int x,int y,int z,int kk1)
{
    if(sta>=x&&enn<=y)
    {
        if(z<0ll)//查询
        {
            ll ans2=segtr[0][now];
            for(int i=1;i<kk1+1;i++)//关键
                ans2=((ans2*((ll)z+1ll)+(ll)dir[kk1][i]*segtr[i][now])%Mod+Mod)%Mod;
            ans=(ans+ans2)%Mod;
        }
        else//修改
        {
        mark[now]=z;
        for(int i=0;i<6;i++)
            segtr[i][now]=((ll)z*(ini[enn][i]-ini[sta-1][i])%Mod+Mod)%Mod;//更新的关键
        }
        return;
    }
    int mid=dir(sta+enn,1);
    int next=mul(now,1);
    Downnow(now,mid,next,sta,enn);
    if(mid>=x)
        Update(sta,mid,next,x,y,z,kk1);
    if(mid<y)
        Update(mid+1,enn,next|1,x,y,z,kk1);
    Upnow(now,next);
    return;
}
int main()
{
    Init();
int n,m,lef,rig,val;
while(~scanf("%d %d",&n,&m))
{
    Create(1,n,1);
    for(int i=0;i<m;i++)
    {
        scanf("%s %d %d %d",str,&lef,&rig,&val);
        if(str[0]=='=')
            Update(1,n,1,lef,rig,val,0);
        else
        {
            ans=0ll;
            Update(1,n,1,lef,rig,-lef,val);
            printf("%I64d\n",ans);
        }
    }
}
    return 0;
}
版权声明:本文为博主原创文章,未经博主允许不得转载。

【CodeForces】266E More Queries to Array... 线段树

E. More Queries to Array... time limit per test 5 seconds memory limit per test 256 megabyte...
  • u013368721
  • u013368721
  • 2014年06月26日 14:23
  • 787

CodeForces 266E More Queries to Array...

因为k最大只有5,所以我建六个线段树,分别用于保存ai*i^k的值,这样就能在logn内算出值,所以此题总复杂度nlogn。 #include #include #include #inclu...
  • sssogs
  • sssogs
  • 2013年03月27日 00:28
  • 911

CodeForces 266E More Queries to Array... 线段树

题意:两种操作: 1 区间【l,r】变为一个值 2 查询   ,k 解题:线段树每个点记录a[l]*l^k,统计区间和,直接查询即可。然后把要查询的值多项式展开转化一下,当k=5时,就要同时查询...
  • translata
  • translata
  • 2016年02月27日 16:21
  • 110

CF 266E More Queries to Array...(线段树)

转载请注明出处,谢谢http://blog.csdn.net/ACM_cxlove?viewmode=contents    by---cxlove 题目:给出一个数列,修改操作是将区间的数改为...
  • ACM_cxlove
  • ACM_cxlove
  • 2013年03月20日 20:42
  • 1803

codeforces 266E More Queries to Array 线段树

题目大意:给你一串数列,a[1]~a[n]。有两种操作 操作1:", (1 ≤ l ≤ r ≤ n; 0 ≤ x ≤ 109). 表示将区间[l,r]中的数全赋值为x。 操作2:"", (1 ≤ l...
  • dyx404514
  • dyx404514
  • 2013年03月14日 15:55
  • 763

Codeforces Round #163 (Div. 2) E. More Queries to Array...

题意:有N个数,M个操作。(1)"= l r k",表示把区间[l,r]的数全部变成k。(2)"? l r k",查询区间[l,r]范围里。 k比较小,将式子拆开成多项式,一项一项加。数学的东西感觉...
  • shiqi_614
  • shiqi_614
  • 2013年01月23日 16:19
  • 1393

codeforces 482B. Interesting Array【线段树区间更新】

题目:codeforces 482B. Interesting Array 题意:给你一个值n和m中操作,每种操作就是三个数 l ,r,val。就是区间l---r上的与的值为val,最后问你...
  • y990041769
  • y990041769
  • 2014年10月29日 19:34
  • 2043

Codeforces #371(Div.2)C.Sonya and Queries【map应用】【思维】

传送门 :C. Sonya and Queries 描述: C. Sonya and Queries time limit per test 1 second memory limit p...
  • guhaiteng
  • guhaiteng
  • 2016年09月14日 14:13
  • 490

codeforces 266E 多颗线段树维护

http://codeforces.com/contest/266/problem/E   题意明显,略掉。把表达式展开,其中i一项,(1-l)一项,根据二项式定理展开,发现只需要维护a[...
  • liverpippta
  • liverpippta
  • 2013年08月04日 15:52
  • 833

Codeforces 650D. Zip-line (动态LIS) (可持久化线段树 或 离线+树状数组)

题意: 给定一个长度为n的数列,和m个询问,每个询问的格式是:将原数组的第a个数改成b之后,数组的最长上升子序列(LIS)的长度。 做法:可持久化线段树 令 h[ ] 为原数组,LIS_L[ i ...
  • u012891242
  • u012891242
  • 2016年03月11日 17:38
  • 1346
内容举报
返回顶部
收藏助手
不良信息举报
您举报文章:CodeForces 266E More Queries to Array...(线段树+式子展开)
举报原因:
原因补充:

(最多只允许输入30个字)