数据结构小总结（成都磨子桥技工学校数据结构前12题）

[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]))%