【BZOJ 2038】[2009国家集训队]小Z的袜子(hose)

Description

HH有一串由各种漂亮的贝壳组成的项链。HH相信不同的贝壳会带来好运，所以每次散步 完后，他都会随意取出一段贝壳，思考它们所表达的含义。HH不断地收集新的贝壳，因此， 他的项链变得越来越长。有一天，他突然提出了一个问题：某一段贝壳中，包含了多少种不同 的贝壳？这个问题很难回答。。。因为项链实在是太长了。于是，他只好求助睿智的你，来解 决这个问题。

Input

第一行：一个整数N，表示项链的长度。 第二行：N个整数，表示依次表示项链中贝壳的编号（编号为0到1000000之间的整数）。 第三行：一个整数M，表示HH询问的个数。 接下来M行：每行两个整数，L和R（1 ≤ L ≤ R ≤ N），表示询问的区间。

Output

M行，每行一个整数，依次表示询问对应的答案。

Sample Input

6
1 2 3 4 3 5
3
1 2
3 5
2 6

Sample Output

2
2
4

HINT

对于20%的数据，N ≤ 100，M ≤ 1000；
对于40%的数据，N ≤ 3000，M ≤ 200000；
对于100%的数据，N ≤ 50000，M ≤ 200000。

 

莫队算法。。。感觉是一种有规律的暴力

他们写的比我好，不想多说什么了

#include<cstdio>
#include<cmath>
#include<algorithm>
#define ll long long
using namespace std;
const int N=50010;
struct node{int l,r,id;
            ll a,b;}a[N];
int n,m;
ll ans,sum[N];
int c[N],bel[N];
bool cmp(node a,node b){
    if (bel[a.l]==bel[b.l]) return (a.r<b.r);
    return a.l<b.l;
}
 
bool cmp1(node a,node b){return a.id<b.id;}
ll sqr(ll x){return x*x;}
void init(){
    scanf("%d%d",&n,&m);
    for (int i=1;i<=n;i++) scanf("%d",&c[i]);
    int block=int(sqrt(n));
    for (int i=1;i<=n;i++) bel[i]=(i/block)+1;
    for (int i=1;i<=m;i++){
        scanf("%d%d",&a[i].l,&a[i].r);
        a[i].id=i;
    }
}
 
void updata(int wz,int ch){
    ans-=sqr(sum[c[wz]]);
    sum[c[wz]]+=ch;
    ans+=sqr(sum[c[wz]]);
}
 
ll gcd(ll a,ll b){
    if (a%b==0) return b;
    return gcd(b,a%b);
}
 
int main(){
    init();
    sort(a+1,a+m+1,cmp);
    int l=1,r=0;
    for (int i=1;i<=m;i++){
        for (;r<a[i].r;r++) updata(r+1,1);
        for (;r>a[i].r;r--) updata(r,-1);
        for (;l<a[i].l;l++) updata(l,-1);
        for (;l>a[i].l;l--) updata(l-1,1);
        if (a[i].l==a[i].r) {a[i].a=0,a[i].b=1;continue;}
        a[i].a=ans-(a[i].r-a[i].l+1);
        a[i].b=(ll)(a[i].r-a[i].l+1)*(a[i].r-a[i].l);
        ll tmp=gcd(a[i].a,a[i].b);
        a[i].a/=tmp;a[i].b/=tmp;
    }
    sort(a+1,a+m+1,cmp1);
    for (int i=1;i<=m;i++) printf("%lld/%lld\n",a[i].a,a[i].b);
}