[pixiv] https://www.pixiv.net/member_illust.php?mode=medium&illust_id=34352147
暑假的作业,颓颓的我总算是写完了
线段树
线段树是一个高级玩意,不仅可以求区间和,区间最大等等的简单问题,灵活运用还有好多变种。自从学了主席树,知道了null自环这种东西后,用在线段树上也是得心应手
c3
给一个长为N的数列,有M次操作,每次操作是以下两种之一:
(1)修改数列中的一个数
(2)求数列中某连续一段的和
赤裸裸的线段树
c4
给一个长为N的数列,有M次操作,每次操作时以下三种之一:
(1)修改数列中的一个数
(2)求数列中某连续一段所有数的两两乘积的和 mod 1000000007
(3)求数列中某连续一段所有相邻两数乘积的和 mod 1000000007
数据剧毒无比,有负数,取模就出问题了。对于区间维护答案,主要就是如何合并区间。操作3好合并,只要记录每个区间的头、尾的数,把左右儿子区间的和加起来,再加上中间两个数的乘积。关键是操作2,要是没有见识过这个脑筋急转弯,我可能一辈子都不会:
给出N个数, 每次可以合并两个数, 合并的代价是两个数的乘积, 合并得到的数是两个数的和。
问最后把所有数合并成一个数的最小代价。 求这个最小代价对10^9+7取模的结果。
N <= 5000000。
题解是:
显然无论怎么合并答案都是一样的, 任意两个数的乘积恰好会对答案贡献一次。
直接搞就好了
于是这道题的操作2就迎刃而解了
由于这道题坑很多,我就不放我wa掉的代码了
大神的AC代码
#include<iostream>
#include<cstdio>
#include<algorithm>
using namespace std;
const int maxn=100005;
const int mod=1000000007;
struct Node{
int s1,s2,s3;
}node[maxn<<2];
int a[maxn];
int n,m;
int slow_mult(int a,int p){
if(a<0)a=(a+(mod<<1))%mod;
if(p<0)p =(p+(mod<<1))%mod;
if(a<p)swap(a,p);
if(p==0)return 0;
if(p==1)return a%mod;
int tmp=slow_mult(a,p>>1);
if(p&1)return ((tmp<<1)%mod+a%mod)%mod;
else return (tmp<<1)%mod;
}
inline void update(int root,int l,int r){
node[root].s1=(node[root<<1].s1+node[root<<1|1].s1)%mod;
node[root].s2=(node[root<<1].s2+node[root<<1|1].s2)%mod;
int m=(l+r)>>1;
node[root].s3=((node[root<<1].s3+node[root<<1|1].s3)%mod+slow_mult(a[m],a[m+1]))%mod;
}
void build(int root,int l,int r){
if(l==r){
node[root].s1=a[l];
node[root].s2=slow_mult(a[l],a[r]);
node[root].s3=0;
return ;
}
int m=(l+r)>>1;
build(root<<1,l,m),build(root<<1|1,m+1,r);
update(root,l,r);
}
void modify(int root,int l,int r,int x,int val){
if(l==r){
node[root].s1=val;
node[root].s2=slow_mult(val , val);
node[root].s3=0;
a[x]=val;
return;
}
int m=(l+r)>>1;
if(x<=m)modify(root<<1,l,m,x,val);
else modify(root<<1|1,m+1,r,x,val);
update(root,l,r);
}
int query1(int root,int l,int r,int x,int y){
if(x<=l&&r<=y)return node[root].s1;
int m=(l+r)>>1,ret=0;
if(x<=m&&l<=y)ret+=query1(root<<1,l,m,x,y);
if(y>=m+1&&r>=x)ret+=query1(root<<1|1,m+1,r,x,y);
return ret % mod;
}
int query2(int root,int l,int r,int x,int y){
if(x<=l&&r<=y)return node[root].s2;
int m=(l+r)>>1,ret=0;
if(x<=m&&l<=y)ret+=query2(root<<1,l,m,x,y);
if(y>=m+1&&r>=x)ret+=query2(root<<1|1,m+1,r,x,y);
return ret % mod;
}
int query3(int root,int l,int r,int x,int y){
if(x<=l&&r<=y)return node[root].s3;
int m=(l+r)>>1,ret=0,flag=1;
if(x<=m&&l<=y)ret+=query3(root<<1,l,m,x,y),flag*=-1;
if(y>=m+1&&r>=x)ret+=query3(root<<1|1,m+1,r,x,y),flag*=-1;
if(flag^1)return ret%mod;
else return(ret%mod+slow_mult(a[m],a[m+1]))%