CF848,第一次参加div1

开场3分钟杜教秒了A题,未成年的内心收到了震撼。
感觉A题同种字符只要数量一定,对答案的贡献就一定,而且不同种字符对答案的贡献互不影响。
然后打表找了下规律,好像是f[x]=x*(x-1)>>1
接着一波01背包就过了A题,不过感觉 26nn 可能比较悬?幸好这是CF。
看了下榜,大概300名左右,我还是太菜了。
开B题,似乎两个人i,j会相交,就一定满足ti-pi==tj-pj,开个map存一下每个权值ti-pi对应的人就行了,如果有多个,就取最后一个。因为是二维平面,所以要做两遍
写完后发现WA on test5了,好像同一个权值可以不只有一个人,而且除了最后一个,别的也会受到影响,轨迹发生变化。
map<int,int>换成map<int,deque<int> >,发现WA on test4,好像第二遍时map忘清空了,改完就A了。
此时(1:23)还是300名左右。
UOJ群里说C是个动态二维数点,还剩30分钟时看懂题并想清如何动态二维数点,码码码,还剩5分钟时写完,还剩2分钟时发现bug,然后,就没有然后了,Orz lych_cys 24分钟一遍A掉 C题。
事后花了20分钟调完C题,发现居然卡内存(我写的是树套树),然后了20分钟卡内存,将sum为0的节点给垃圾回收,然后才AC
这样算下来或许赛场上不做B,而是去做C可能更划算?
B 1000分,耗时65min WA2发
C 1750分,耗时75min左右,RE2发
突然有些后悔。
下面是蒟蒻的代码。

#include<cstdio>
#include<cctype>
#include<set>
using namespace std;
const int N=100010,M=12000000;
typedef long long ll;
struct memorypool{
    int a[M],xb;
    inline void Init(){
        for(register int i=0;i<M;++i)a[i]=i;
        xb=M-1;
    }
    inline int New(){
        return a[xb--];
    }
    inline void Del(int x){
        a[++xb]=x;
    }
}X;
struct node{
    int l,r;
    ll s;
}a[M];
int xb;
void add(int&i,int rl,int rr,int l,int v){
    if(!i)i=X.New();
    a[i].s+=v;
    if(rl==rr)return;
    int m=(rl+rr)>>1;
    if(l>m)add(a[i].r,m+1,rr,l,v);
        else add(a[i].l,rl,m,l,v);
    if(!a[i].s)X.Del(i),i=0;
}
ll query(int i,int rl,int rr,int l,int r){
    if(!i)return 0;
    if(rl==l && rr==r)return a[i].s;
    int m=(rl+rr)>>1;
    if(l>m)return query(a[i].r,m+1,rr,l,r);
        else if(r<=m)return query(a[i].l,rl,m,l,r);
                else return query(a[i].l,rl,m,l,m)+query(a[i].r,m+1,rr,m+1,r);
}
int n,i,aa[N],lst[N],nxt[N],m,o,u,v;
struct BIT{
    int rt[N];
    inline void addd(int x,int y){//lst[x] x
        for(int i=x+1;i<=n+2;i+=i&-i)add(rt[i],1,n,y,y);
    }
    inline void del(int x,int y){
        for(int i=x+1;i<=n+2;i+=i&-i)add(rt[i],1,n,y,-y); 
    }
    inline ll que(int x,int y,int z){//<=x y->z
        ll ans=0;
        for(int i=x+1;i;i-=i&-i)ans+=query(rt[i],1,n,y,z);  
        return ans;
    }
}tl,tn;
set<int> s[N];
set<int>::iterator it1,it2;
inline char read() {
    static const int IN_LEN = 1000000;
    static char buf[IN_LEN], *s, *t;
    if (s == t) {
        t = (s = buf) + fread(buf, 1, IN_LEN, stdin);
        if (s == t) return -1;
    }
    return *s++;
}
template<class T>
inline void read(T &x) {
    static bool iosig;
    static char c;
    for (iosig = false, c = read(); !isdigit(c); c = read()) {
        if (c == '-') iosig = true;
        if (c == -1) return;
    }
    /* 这里这么写的原因见下 */
    for (x = 0; isdigit(c); c = read())
        x = (x + (x << 2) << 1) + (c ^ '0');
    if (iosig) x = -x;
}
const int OUT_LEN = 2000000;
char obuf[OUT_LEN], *oh = obuf;
inline void print(char c) {
    if (oh == obuf + OUT_LEN) fwrite(obuf, 1, OUT_LEN, stdout), oh = obuf;
    *oh++ = c;
}
template<class T>
inline void print(T x) {
    static int buf[30], cnt;
    if (x == 0) {
        print('0');
    } else {
        if (x < 0) print('-'), x = -x;
        for (cnt = 0; x; x /= 10) buf[++cnt] = x % 10 + 48;
        while (cnt) print((char)buf[cnt--]);
    }
}
inline void flush() {
    fwrite(obuf, 1, oh - obuf, stdout);
}
int main(){
    //freopen("c.txt","r",stdin);
    //freopen("cc.txt","w",stdout);
    X.Init();
    read(n);
    read(m);
    for(i=1;i<=n;++i)s[i].insert(0),s[i].insert(n+1);
    for(i=1;i<=n;++i){
        read(aa[i]);
        lst[i]=*--s[aa[i]].lower_bound(i);
        s[aa[i]].insert(i);
    }
    for(i=1;i<=n;++i)nxt[i]=*s[aa[i]].upper_bound(i);
    for(i=1;i<=n;++i){
        tl.addd(lst[i],i);
        tn.addd(nxt[i],i);
    }
    while(m--){
        read(o),read(u),read(v);
        if(o&2)print(tn.que(n+1,u,v)-tn.que(v,u,v)-tl.que(u-1,u,v)),print('\n');
            else{
                it1=s[aa[u]].upper_bound(u);
                it2=--s[aa[u]].lower_bound(u);
                if(*it1<=n)tl.del(lst[*it1],*it1),lst[*it1]=*it2,tl.addd(lst[*it1],*it1);
                if(*it2>=1)tn.del(nxt[*it2],*it2),nxt[*it2]=*it1,tn.addd(nxt[*it2],*it2);
                s[aa[u]].erase(u);
                it1=s[v].upper_bound(u);
                it2=--s[v].lower_bound(u);
                if(*it1<=n)tl.del(lst[*it1],*it1),lst[*it1]=u,tl.addd(lst[*it1],*it1);
                if(*it2>=1)tn.del(nxt[*it2],*it2),nxt[*it2]=u,tn.addd(nxt[*it2],*it2);          
                s[v].insert(u);
                tl.del(lst[u],u);
                tn.del(nxt[u],u);
                lst[u]=*it2;
                nxt[u]=*it1;
                tl.addd(lst[u],u);
                tn.addd(nxt[u],u);  
                aa[u]=v;
            }
    }
    flush();
    return 0;
}
  • 0
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值