莫隊算法：區間衆數出現次數

2020.5.19
昨天學完莫隊之後重寫了一遍mcpc19的那一道區間衆數的，發現了幾個可以改進的點。一個是離散化，港真我不知道原題是怎麼回事兒，上了離散化，後來無論是d-query還是路邊的蒲公英好像完全沒必要上離散化，畢竟是1e5個數字，開個數組記錄下就行了，所以這部分刪掉，然後模擬一下我們查詢區間衆數出現次數的方法，首先考慮加入座標位置爲x，數字爲a[x]，對於當前區間維護的答案的影響，首先要消除對於上一次a[x]出現的影響，那就是sum[cnt[a[x]]] -1， 然後將出現次數cnt[a[x]]加一，比較其與目前出現次數的的大小，然後將sum[cnt[a[x]]] 加一，就完成了一次維護，縮進的del函數同理，就是看sum[cnt[a[x]]] 是不是爲0，如果cnt[a[x]]是目前最大次數，（因爲sum[cnt[a[x]]] 是0所以是唯一最大次數），那麼就對答案進行-1，修正後的區間就是我們所需要的答案。
代碼如下：

#include <bits/stdc++.h>
using namespace std;
#define limit (10000 + 5)//防止溢出
#define INF 0x3f3f3f3f
#define inf 0x3f3f3f3f3f
#define lowbit(i) i&(-i)//一步兩步
#define EPS 1e-6
#define FASTIO  ios::sync_with_stdio(false);cin.tie(0);
#define ff(a) printf("%lld\n",a );
#define pi(a,b) pair<a,b>
#define rep(i, a, b) for(int i = a ; i <= b ; ++i)
#define per(i, a, b) for(int i = b ; i >= a ; --i)
#define mint(a,b,c) min(min(a,b), c)
#define MOD 998244353
#define FOPEN freopen("C:\\Users\\administrator01\\CLionProjects\\untitled24\\data.txt", "rt", stdin)
typedef long long ll;
typedef unsigned long long ull;
ll read(){
    ll sign = 1, x = 0;char s = getchar();
    while(s > '9' || s < '0' ){if(s == '-')sign = -1;s = getchar();}
    while(s >= '0' && s <= '9'){x = x * 10 + s - '0';s = getchar();}
    return x * sign;
}//快讀
void write(ll x){
    if(x / 10) write(x / 10);
    putchar(x % 10 + '0');
}
int n, m;
struct node{
    int l, r, qid, blo;
    int len(){
        return r - l + 1;
    }
    bool operator<(const node &rhs)const{
        if(blo ^ rhs.blo)return l < rhs.l;
        return blo & 1 ? r < rhs.r : rhs.r < r;
    }
}query[limit];
int cnt[limit],a[limit],sum[limit];
int res;
void add(int x){
    --sum[cnt[a[x]]];//減去多餘部分
    ++cnt[a[x]];//增加一個count
    res = max(res, cnt[a[x]]);
    ++sum[cnt[a[x]]];
}
void del(int x){
    --sum[cnt[a[x]]];//撤回一個count
    if(res == cnt[a[x]] && !sum[cnt[a[x]]])--res;//如果沒有了那就增加
    --cnt[a[x]];
    ++sum[cnt[a[x]]];
}
map<double, int>mp;
int tot, vis[limit];
int main(){
#ifdef LOCAL
    FOPEN;
    //freopen("C:\\Users\\administrator01\\CLionProjects\\untitled24\\out.txt", "w", stdout);
#endif
    n = read(), m = read(),
    tot = 0;
    rep(i ,1, n){
        double input;
        scanf("%lf" , &input);//掃一下
        if(!mp[input]){
            a[i] = mp[input] = ++tot;
        }else{
            a[i] = mp[input];
        }
    }
    int block = int(sqrt(n >= 3 ? n * (2.0 / 3) : n));//分塊
    rep(i, 1, m){
        query[i].l = read(), query[i].r = read(), query[i].qid = i;
        query[i].blo = query[i].l / block;//分成塊
    }
    sort(query + 1, query + 1 + m);
    int l = 0 , r = 0;
    res = 0;
    rep(i ,1, m){
        while(l < query[i].l)del(l++);//縮進
        while(l > query[i].l)add(--l);
        while(r < query[i].r)add(++r);
        while(r > query[i].r)del(r--);
        vis[query[i].qid] = res >= query[i].len() / 2 + 1;
    }
    rep(i ,1, m){
        puts(vis[i] ? "usable" : "unusable");
    }
    return 0;
}