【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=;
 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=;i<=n;i++) scanf("%d",&c[i]);
     int block=int(sqrt(n));
     for (int i=;i<=n;i++) bel[i]=(i/block)+;
     for (int i=;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==) return b;
     return gcd(b,a%b);
 }

 int main(){
     init();
     sort(a+,a+m+,cmp);
     int l=,r=;
     for (int i=;i<=m;i++){
         for (;r<a[i].r;r++) updata(r+,);
         for (;r>a[i].r;r--) updata(r,-);
         for (;l<a[i].l;l++) updata(l,-);
         for (;l>a[i].l;l--) updata(l-,);
         if (a[i].l==a[i].r) {a[i].a=,a[i].b=;continue;}
         a[i].a=ans-(a[i].r-a[i].l+);
         a[i].b=(ll)(a[i].r-a[i].l+)*(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+,a+m+,cmp1);
     for (int i=;i<=m;i++) printf("%lld/%lld\n",a[i].a,a[i].b);
 }