bzoj 2038 小Z的袜子

题目链接

分析：

这道题就是暴力找，但是有一种方法，通过巧妙地重新整理查询的次序，来降低复杂度。这种算法叫莫队算法，适用于离线的区间查询问题。

如果我们已经知道了对于查询[L, R]的答案，那么对于这道题，[L-1,R],[L+1,R]等等的问题查询其实只需要$O(1)$的转移，这种情况的计算就可以用莫队算法巧妙地解决。

方法是这样的，对所有查询离线，然后重新分块排序。对于$N$个查询，我们分出$\sqrt{N}$ 个块，按照L为第一关键字分块，R为第二关键字在块内排序，将$N$个查询排序。

然后就是用currentL和currentR两个指针暴力地查找了。这样的复杂度就不是$O(N^2)$而是$O(N^{1.5})$了。这样的，对于每个查询，左指针最多挪动$\sqrt{N}$个长度，右指针在块内递增，在块之间最多挪动$N$个长度，而一共最多有$\sqrt{N}$个块。这样总的复杂度就是$O(N\sqrt{N})$。很好理解。

---

很好理解啊＝ ＝
个人感觉在区间的端点需要好好想清楚，嗯，这样。
关于起始位置curl＝1，curr＝0，这是非常讲究的，当初看别人这么写不大理解，自己写了一遍就明白了原因。这样正好能够避免update函数在初始位置可能会发生的错误。

代码：

/*****************************************************/
//#pragma comment(linker, "/STACK:1024000000,1024000000")
#include <map>
#include <set>
#include <stack>
#include <queue>
#include <cmath>
#include <string>
#include <vector>
#include <cstdio>
#include <cstring>
#include <sstream>
#include <cstdlib>
#include <iostream>
#include <algorithm>

using namespace std;

#define   offcin        ios::sync_with_stdio(false)
#define   DEBUG         freopen("#.txt", "r", stdin)
#define   sigma_size    26
#define   lson          l,m,v<<1
#define   rson          m+1,r,v<<1|1
#define   slch          v<<1
#define   srch          v<<1|1
#define   ll            long long
#define   ull           unsigned long long
#define   lowbit(x)     (x&-x)

const int    INF    = 0x3f3f3f3f;
const ll     INFF   = 1e18;
const double pi     = acos(-1.0);
const double inf    = 1e18;
const double eps    = 1e-9;
const ll     mod    = 1e9+7;
const int    maxmat = 10;
const ull    BASE   = 133333331;

/*****************************************************/
inline void RI(int &x) {
      char c;
      while((c=getchar())<'0' || c>'9');
      x=c-'0';
      while((c=getchar())>='0' && c<='9') x=(x<<3)+(x<<1)+c-'0';
}
inline ll bits(ll x) {
      return !x ? x : bits(x - lowbit(x)) + 1;
}
/*****************************************************/

const int maxn = 5e4 + 5;

int c[maxn];
int N, M;
int b[maxn];
long long cnt[maxn];
long long res, ans;
int cur_l, cur_r;
pair<long long, long long> a[maxn];

struct Node {
    int l, r, id;
    bool operator <(const Node &rhs) const {
        if (b[l] == b[rhs.l])
            return r < rhs.r;
        return b[l] < b[rhs.l];
    }
}node[maxn];

void update(long long &col, int x) {
    if (x == 1) {
        res += col;
        col ++;
    }
    else {
        col --;
        res -= col;
    }
}

void solve(int L, int R, int id) {
    while (cur_r < R) update(cnt[c[++ cur_r]], 1);
    while (cur_r > R) update(cnt[c[cur_r --]], -1);
    while (cur_l > L) update(cnt[c[-- cur_l]], 1);
    while (cur_l < L) update(cnt[c[cur_l ++]], -1);
    ans = R - L + 1;
    ans = ans * (ans - 1) / 2;
    int k = __gcd(res, ans);
    if (res == 0) ans = 1, k = 1;
    a[id] = make_pair(res / k, ans / k);
}

int main(int argc, char const *argv[]) {
    scanf("%d%d", &N, &M);

    int block = ceil(sqrt((double)N));
    for (int i = 1; i <= N; i ++) {
        scanf("%d", c + i);
        b[i] = (i - 1) / block;
    }

    for (int i = 1; i <= M; i ++) {
        scanf("%d%d", &node[i].l, &node[i].r);
        node[i].id = i;
    }

    sort(node + 1, node + M + 1);
    cur_l = 1, cur_r = 0;
    res = ans = 0;
    for (int i = 1; i <= M; i ++)
        solve(node[i].l, node[i].r, node[i].id);

    for (int i = 1; i <= M; i ++)
        printf("%lld/%lld\n", a[i].first, a[i].second);

    return 0;
}