CodeChef June Challenge 2017

好气啊，本来以为比赛时间还有很多，结果回家养病两天回到学校怎么比赛就结束了（雾），大约是小高考弄错了时间？

挑3道有意思的写写题解吧。

Cloning

题目大意：给一个序列，每次询问两个等长区间，问区间内的数是否排序后至多只有一个对应位不同。

题解：主席树维护一下hash，求出最大的k，使得两个区间中的前k大在排序后相同，然后判一下后缀即可。

#include<cstdio>
#include<algorithm>
#define MN 110000
#define ull unsigned long long
using namespace std;

int read_p,read_ca;
inline int read(){
    read_p=;read_ca=getchar();
    while(read_ca<''||read_ca>'') read_ca=getchar();
    while(read_ca>=''&&read_ca<='') read_p=read_p*+read_ca-,read_ca=getchar();
    return read_p;
}
const ull base=;
int T,n,m,ro[MN],num=,a,b,c,d;
ull H[MN];
struct na{int l,r,s;ull h;na(){s=h=;}}t[MN*];
inline void add(int &p,int x,int l,int r,int k){
    p=++num;t[p]=t[x];
    t[p].h+=H[k];t[p].s++;
    if (l==r) return;
    int mid=l+r>>;
    if (k<=mid) add(t[p].l,t[x].l,l,mid,k);else add(t[p].r,t[x].r,mid+,r,k);
}
inline int aski(int a,int b,int c,int d,int l,int r){
    if (t[b].h-t[a].h==t[d].h-t[c].h&&t[b].s-t[a].s==t[d].s-t[c].s) return t[b].s-t[a].s;
    if (l==r) return min(t[b].s-t[a].s,t[d].s-t[c].s);
    int mid=l+r>>;
    int mmh=aski(t[a].l,t[b].l,t[c].l,t[d].l,l,mid);
    if (mmh==t[t[b].l].s-t[t[a].l].s&&mmh==t[t[d].l].s-t[t[c].l].s) mmh+=aski(t[a].r,t[b].r,t[c].r,t[d].r,mid+,r);
    return mmh;
}
inline int aska(int a,int b,int c,int d,int l,int r){
    if (t[b].h-t[a].h==t[d].h-t[c].h&&t[b].s-t[a].s==t[d].s-t[c].s) return t[b].s-t[a].s;
    if (l==r) return min(t[b].s-t[a].s,t[d].s-t[c].s);
    int mid=l+r>>;
    int mmh=aska(t[a].r,t[b].r,t[c].r,t[d].r,mid+,r);
    if (mmh==t[t[b].r].s-t[t[a].r].s&&mmh==t[t[d].r].s-t[t[c].r].s) mmh+=aska(t[a].l,t[b].l,t[c].l,t[d].l,l,mid);
    return mmh;
}
int main(){
    int i;
    T=read();
    H[]=;
    for (i=;i<=1e5;i++) H[i]=H[i-]*base;
    while(T--){
        n=read();m=read();num=;
        for (i=;i<=n;i++) add(ro[i],ro[i-],,1e5,read());
        while(m--){
            a=read();b=read();c=read();d=read();
            puts(aska(ro[a-],ro[b],ro[c-],ro[d],,1e5)+aski(ro[a-],ro[b],ro[c-],ro[d],,1e5)>=b-a?"YES":"NO");
        }
    }
} 

Euler Sum

题目大意：求$\sum_{1}^{n}\left\lfloor i*e \right\rfloor$

题解：连分数求出e的近似分数以后类欧即可

z,m=1,0

def work(x):
    global z,m
    z,m=m,z
    z+=m*x
def mmh(n,z,m):
    if n==0 or z==0 or m==0:
        return 0
    p=int(z//m)
    z-=int(p*m)
    N=int(n*z//m)
    return n*(n+1)//2*p+int(n*N)-int(mmh(N,m,z))

i=3000
while (i>0):
    work(i*2)
    work(1)
    work(1)
    i-=1
z+=m
n=int(input())
mmh=0
f=1
while(n>0):
    p=z//m
    z-=p*m
    N=n*z//m
    mmh+=f*(n*(n+1)//2*p+n*N)
    f*=-1
    n=N
    z,m=m,z
print(mmh) 

Persistent oak

题意：不想写啊

题解：链剖后可持久化线段树即可，注意打标记，下放标记的时候都要新建节点，不然会影响到历史版本。

#include<cassert>
#include<cstdio>
#include<cstring>
#include<algorithm>
#define MN 210000
#define lp t[p].l
#define rp t[p].r
using namespace std;

int read_p,read_ca;
inline int read(){
    read_p=;read_ca=getchar();
    while(read_ca<''||read_ca>'') read_ca=getchar();
    while(read_ca>=''&&read_ca<='') read_p=read_p*+read_ca-,read_ca=getchar();
    return read_p;
}
struct na{int y,ne;}b[MN<<];
struct tree{int l,r,M,m,d;bool re;}t[MN*];
int T,n,m,ty,l[MN],fa[MN],si[MN],de[MN],son[MN],to[MN],x,y,w[MN],df[MN],pdf[MN],lo[MN],ro[MN],nm,NUM,num;
char s[];
inline void in(int x,int y){b[++num].y=y;b[num].ne=l[x];l[x]=num;}
inline int min(int a,int b){return a<b?a:b;}
void dfs(int x){
    si[x]=;son[x]=;
    for (int i=l[x];i;i=b[i].ne)
    if (fa[b[i].y]=x,de[b[i].y]=de[x]+,dfs(b[i].y),si[x]+=si[b[i].y],si[son[x]]<si[b[i].y]) son[x]=b[i].y;
}
void DFS(int x,int p){
    to[x]=p;df[x]=++nm;pdf[nm]=x;
    if (son[x]) DFS(son[x],p);
    for (int i=l[x];i;i=b[i].ne)
    if (b[i].y!=son[x]) DFS(b[i].y,b[i].y);
    lo[x]=nm;
}
inline void hb(int &p,int d,bool bo=){
    if (!p) return;
    t[++NUM]=t[p];p=NUM;
    if (bo) t[p].re=,t[p].m=t[p].M,t[p].d=;else t[p].m+=d,t[p].d+=d,assert(t[p].m<=t[p].M);
}
inline void pd(int &p){
    t[++NUM]=t[p];p=NUM;
    if (t[p].re) hb(lp,,),hb(rp,,),t[p].re=;
    if (t[p].d) hb(lp,t[p].d),hb(rp,t[p].d),t[p].d=;
}
inline void gx(int &p){
    t[p].m=min(t[lp].m,t[rp].m);
    t[p].M=min(t[lp].M,t[rp].M);
}
void build(int &p,int l,int r){
    p=++NUM;t[p].l=t[p].r=t[p].d=t[p].re=;
    if (l==r){
        t[p].M=t[p].m=w[pdf[l]];
        return;
    }
    int mid=l+r>>;
    build(t[p].l,l,mid);build(t[p].r,mid+,r);
    gx(p);
}
void add(int &p,int x,int l,int r,int L,int R,int d){
    p=++NUM;t[p]=t[x];
    pd(p);x=p;
    if (L==l&&R==r) hb(p,d);else{
        int mid=l+r>>;
        if (R<=mid) add(lp,t[x].l,l,mid,L,R,d);else
        if (L>mid) add(rp,t[x].r,mid+,r,L,R,d);else
        add(lp,t[x].l,l,mid,L,mid,d),add(rp,t[x].r,mid+,r,mid+,R,d);
        gx(p);
    }
}
int ask(int &p,int l,int r,int k){
    pd(p);
    if (l==r) return t[p].m;
    int mid=l+r>>;
    if (k<=mid) return ask(lp,l,mid,k);else return ask(rp,mid+,r,k);
}
int mask(int &p,int l,int r,int L,int R){
    pd(p);
    if (t[p].m>=) return -;
    if (l==r) return pdf[l];
    int mid=l+r>>,s;
    if (R<=mid) return mask(lp,l,mid,L,R);else
    if (L>mid) return mask(rp,mid+,r,L,R);else
    if (s=mask(rp,mid+,r,mid+,R),s!=-) return s;else return mask(lp,l,mid,L,mid);
}
inline void reset(int &p,int x,int l,int r,int L,int R){
    p=++NUM;t[p]=t[x];
    pd(p);x=p;
    if (L==l&&R==r) hb(p,,);else{
        int mid=l+r>>;
        if (R<=mid) reset(lp,t[x].l,l,mid,L,R);else
        if (L>mid) reset(rp,t[x].r,mid+,r,L,R);else
        reset(lp,t[x].l,l,mid,L,mid),reset(rp,t[x].r,mid+,r,mid+,R);
        gx(p);
    }
}
inline void dec(int &p,int x,int u,int d){
    while (u){
        add(p,x,,n,df[to[u]],df[u],d);x=p;
        u=fa[to[u]];
    }
}
inline void work(int &p,int x,int y,int d){
    int u=y,s=-;
    while (u){
        add(p,x,,n,df[to[u]],df[u],-d);x=p;
        if (s==-) s=mask(p,,n,df[to[u]],df[u]);
        u=fa[to[u]];
    }
    if (s!=-){
        printf("%d\n",s-);
        y=w[s]-ask(p,,n,df[s]);
        reset(p,p,,n,df[s],lo[s]);
        dec(p,p,fa[s],y);
    }else puts("");
}
int main(){
    int i;
    T=read();
    while (T--){
        n=read();m=read();n++;w[]=1e9;NUM=num=nm=;
        memset(ro,,sizeof(ro));memset(l,,sizeof(l));
        for (i=;i<=n;i++) in(read()+,i),w[i]=read();
        dfs();DFS(,);build(ro[],,n);
        for (i=;i<=m;i++){
            x=read();ty=read();
            if (ty==){
                y=read()+;
                work(ro[i],ro[x],y,read());
            }else{
                y=read()+;
                printf("%d\n",ty=(w[y]-ask(ro[x],,n,df[y])));
                reset(ro[i],ro[x],,n,df[y],lo[y]);
                dec(ro[i],ro[i],fa[y],ty);
            }
        }
    }
}