SCOI2013 day1 多项式的运算(polynomial)

题目描述：

首先注意！！这道题的修改是将x^L到x^R项进行的！！！不是下标！！也就是说mulx操作会让一些项整体右移！！

好吧，不得不说电子科大的题是世界上最坑的题，电子科大的样例是世界上最坑的样例。这道题本不难，但要是看错题了…………………………就只有20分了…………………………更可悲的是，我的线段树都写错了，于是又悲惨地挂零了……………………

好吧，既然要整体右移，那么就很容易想到用splay来维护了。每次mulx只需要把第r项与第r+1项合并，然后在l前添加一项系数为零的项即可。其余操作就非常容易了，记两个标记add和mul，表示这个项被乘上了mul后加上了add，如何打标记和下放标记什么的就不用说了吧………………而query因为只有10次，直接暴力即可。

注意：1.如果你看到这道题以为是线段树，又看到4s的时限，那么请仔细想想题目到底是不是这个意思。

　　　2.致以后的SCOIer：不要以为过样例了就基本没错了，要知道电子科大的样例不管你错成什么样子都能过的。

代码（求众犇轻虐）：

#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<algorithm>
using namespace std;
typedef long long ll;
#define inc(a,b) a+=b,a=a>=MOD?a-MOD:a
const int MOD=20130426;
const int MAXP=200010;

struct _Node
{
    int v,add,mul,size;
    _Node *s[2],*fa;
}pool[MAXP],*stk[MAXP],root,null;

int idx,stat[MAXP];

_Node* new_Node(_Node *fa)
{
    _Node *ret=pool+(idx++);
    ret->v=ret->add=0;
    ret->size=ret->mul=1;
    ret->s[0]=ret->s[1]=&null;
    ret->fa=fa;
    return ret;
}
void _update(_Node *now)
{
    now->size=now->s[0]->size+now->s[1]->size+1;
}
void init()
{
    root.s[0]=new_Node(&root);
    root.s[0]->s[1]=new_Node(root.s[0]);
    _Node *lst=root.s[0]->s[1];
    for (int i=100001;i>=1;--i)
    {
        lst->s[0]=new_Node(lst);
        lst->s[0]->size=i;
        lst=lst->s[0];
    }
    _update(root.s[0]->s[1]);
    _update(root.s[0]);
}
void _push_down(_Node *now)
{
    if (now->mul!=1)
    {
        now->s[0]->mul=(ll)now->s[0]->mul*now->mul%MOD;
        now->s[1]->mul=(ll)now->s[1]->mul*now->mul%MOD;
        now->s[0]->add=(ll)now->s[0]->add*now->mul%MOD;
        now->s[1]->add=(ll)now->s[1]->add*now->mul%MOD;
        now->s[0]->v=(ll)now->s[0]->v*now->mul%MOD;
        now->s[1]->v=(ll)now->s[1]->v*now->mul%MOD;
        now->mul=1;
    }
    if (now->add)
    {
        inc(now->s[0]->add,now->add);
        inc(now->s[0]->v,now->add);
        inc(now->s[1]->add,now->add);
        inc(now->s[1]->v,now->add);
        now->add=0;
    }
}
void _rot(_Node *now,int d)
{
    _Node *s=now->s[d],*fa=now->fa;
    (now->s[d]=s->s[!d])->fa=now;
    (s->s[!d]=now)->fa=s;
    (fa->s[fa->s[1]==now]=s)->fa=fa;
    _update(now);
}
void _splay(_Node *now,_Node *&aim)
{
    _Node *&p=now->fa,*g;
    int dp,dg;
    while (p!=aim && now!=aim)
    {
        g=p->fa;
        if ((dp=p->s[1]==now) == (dg=g->s[1]==p))
            _rot(g,dg),_rot(p,dp);
        else _rot(p,dp),_rot(g,dg);
    }
    if (now!=aim) _rot(aim,aim->s[1]==now);
    _update(now);
}
_Node* select(int k)
{
    _Node *now=root.s[0];
    while (1)
    {
        _push_down(now);
        if (now->s[0]->size==k) return now;
        if (now->s[0]->size<k) k-=now->s[0]->size+1,now=now->s[1];
        else now=now->s[0];
    }
}
void mulx(int a,int b)
{
    _splay(select(b-1),root.s[0]);
    _splay(select(b+1),root.s[0]->s[1]);
    _Node *now=root.s[0]->s[1];
    _push_down(now->s[0]);
    _Node tmp=*(now->s[0]);
    now->s[0]=&null;
    inc(now->v,tmp.v);
    _splay(select(a-1),root.s[0]);
    _splay(select(a),root.s[0]->s[1]);
    _Node *&cur=root.s[0]->s[1]->s[0];
    cur=new_Node(root.s[0]->s[1]);
    _update(root.s[0]->s[1]);
    _update(root.s[0]);
}
void mul(int a,int b,int v)
{
    _splay(select(a-1),root.s[0]);
    _splay(select(b+1),root.s[0]->s[1]);
    _Node *now=root.s[0]->s[1]->s[0];
    now->mul=(ll)now->mul*v%MOD;
    now->add=(ll)now->add*v%MOD;
    now->v=(ll)now->v*v%MOD;
}
void add(int a,int b,int v)
{
    _splay(select(a-1),root.s[0]);
    _splay(select(b+1),root.s[0]->s[1]);
    _Node *now=root.s[0]->s[1]->s[0];
    inc(now->add,v);
    inc(now->v,v);
}
int get_v(int a,int b,int v)
{
    _splay(select(a-1),root.s[0]);
    _splay(select(b+1),root.s[0]->s[1]);
    int now_v=1,ans=0,top;
    stk[top=1]=root.s[0]->s[1]->s[0];
    memset(stat,0,sizeof(stat));
    while (top)
    {
        _Node *now=stk[top];
        if (!stat[now-pool])
        {
            _push_down(now);
            stat[now-pool]=1;
            if (now->s[0]!=&null) 
                stk[++top]=now->s[0];
            continue;
        }
        if (stat[now-pool]==1)
        {
            inc(ans,((ll)now->v*now_v)%MOD);
            stat[now-pool]=2;
            now_v=((ll)now_v*v)%MOD;
            continue;
        }
        if (stat[now-pool]==2)
        {
            if (now->s[1]!=&null)
                stk[++top]=now->s[1];
            stat[now-pool]=3;
            continue;
        }
        --top;
    }
    return ans;
}
int main()
{
    freopen("polynomial.in","r",stdin);
    freopen("polynomial.out","w",stdout);
    char opt[10];
    int x,y,z,n;
    scanf("%d",&n);
    init();
    while (n--)
    {
        scanf("%s",opt);
        switch (opt[0])
        {
            case 'a':
                scanf("%d%d%d",&x,&y,&z);
                add(x+1,y+1,z%MOD);
                break;
            case 'm':
                if (opt[3]=='x')
                {
                    scanf("%d%d",&x,&y);
                    mulx(x+1,y+1);
                }
                else
                {
                    scanf("%d%d%d",&x,&y,&z);
                    mul(x+1,y+1,z%MOD);
                }
                break;
            case 'q':
                scanf("%d",&x);
                printf("%d\n",get_v(1,100001,x));
                break;
        }
    }
    return 0;
}