【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);

 }

秒客网