zoukankan      html  css  js  c++  java
  • 【题解】洛谷P1975排序

    分块,注意重复的值之间的处理。跟普通分块的操作一样的啦,具体可以参见‘不勤劳的图书管理员’。

    #include <bits/stdc++.h>
    using namespace std;
    #define maxn 500000
    #define lowbit(i) i & (-i) 
    #define int long long
    int n, m, cnt, ans, B, c[200][maxn];
    struct node
    {
        int num, id, rank;
    }a[30000];
    
    int read()
    {
        int x = 0;
        char c;
        c = getchar();
        while(c < '0' || c > '9') c = getchar();
        while(c >= '0' && c <= '9') x = x * 10 + c - '0', c = getchar();
        return x;
    }
    
    int Get_ans(int x, int y)
    {
        if(a[x].rank < a[y].rank) return 1;
        else if(a[x].rank > a[y].rank) return -1;
        return 0;
    }
    
    bool cmp1(node a, node b)
    {
        if(a.num != b.num) return a.num < b.num;
        return a.id < b.id;
    }
    
    bool cmp2(node a, node b)
    {
        return a.id < b.id;
    }
    
    void update(int opt, int x, int sum)
    {
        if(!x) return;
        for(int i = x; i <= cnt; i += lowbit(i))
            c[opt][i] += sum;
    }
    
    int query(int opt, int x)
    {
        if(x < 0) return 0;
        int ans = 0;
        for(int i = x; i; i -= lowbit(i))
            ans += c[opt][i];
        return ans;
    }
    
    signed main()
    {
        n = read(); 
        B = sqrt(n);
        for(int i = 1; i <= n; i ++)
            a[i].num = read(), a[i].id = i;
        sort(a + 1, a + 1 + n, cmp1);
        a[0].rank = 1;
        for(int i = 1; i <= n; i ++)
            if(a[i].num == a[i - 1].num) a[i].rank = a[i - 1].rank;
            else a[i].rank = ++ cnt;
        sort(a + 1, a + 1 + n, cmp2);
        for(int i = 1; i <= n; i ++)
            update(i / B, a[i].rank, 1);
        for(int i = n; i >= 1; i --)
        {
            ans += query(n / B + 1, a[i].rank - 1);
            update(n / B + 1, a[i].rank, 1);
        }
        m = read();
        cout << ans << endl;
        for(int i = 1; i <= m; i ++)
        {
            int x = read(), y = read();
            if(x > y) swap(y, x);
            int p = x / B, q = y / B;
            if(a[x].rank < a[y].rank) ans ++;
            else if(a[x].rank > a[y].rank) ans --;
            if(p == q)
            {
                for(int i = x + 1; i < y; i ++)
                    ans += Get_ans(x, i) + Get_ans(i, y);
                swap(a[x], a[y]);
                printf("%lld
    ", ans);
                continue;
            }
            for(int i = x + 1; i < (p + 1) * B; i ++)
                if(i <= n) ans += Get_ans(x, i) + Get_ans(i, y);
            for(int i = q * B; i < y; i ++)
                ans += Get_ans(i, y) + Get_ans(x, i);
            for(int i = p + 1; i < q; i ++)
            {
                ans = (ans - query(i, a[x].rank - 1) + query(i, a[y].rank - 1));
                ans = (ans - query(i, a[x].rank)  + query(i, a[y].rank));
            }
            update(p, a[x].rank, -1), update(q, a[y].rank, -1);
            update(q, a[x].rank, 1), update(p, a[y].rank, 1);
            swap(a[x], a[y]);
            printf("%lld
    ", ans);
        }
        return 0;
    }
  • 相关阅读:
    夜游遂宁滨江路
    易中天讲座十句人生感悟(发人深省,耐人寻味)
    遥望死海
    一直被忽略的成功之道:勤快并非优点,成功需要懒惰
    给别人以宽容,给自己以信心
    合理支配“财富”:经理人运用时间的12种典型模式
    三月的清晨
    学习的三重境界(想成功的人都不可错过)
    上班
    持续开发你的事业智慧:企业家保持冲锋势头的路径
  • 原文地址:https://www.cnblogs.com/twilight-sx/p/8490705.html
Copyright © 2011-2022 走看看