莫队算法+分块 bzoj2038小z的袜子

题面：

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

Time Limit: 20 Sec Memory Limit: 259 MB
Submit: 13997 Solved: 6334
[Submit][Status][Discuss]
Description

作为一个生活散漫的人，小Z每天早上都要耗费很久从一堆五颜六色的袜子中找出一双来穿。终于有一天，小Z再也无法忍受这恼人的找袜子过程，于是他决定听天由命……
具体来说，小Z把这N只袜子从1到N编号，然后从编号L到R(L 尽管小Z并不在意两只袜子是不是完整的一双，甚至不在意两只袜子是否一左一右，他却很在意袜子的颜色，毕竟穿两只不同色的袜子会很尴尬。
你的任务便是告诉小Z，他有多大的概率抽到两只颜色相同的袜子。当然，小Z希望这个概率尽量高，所以他可能会询问多个(L,R)以方便自己选择。

Input

输入文件第一行包含两个正整数N和M。N为袜子的数量，M为小Z所提的询问的数量。接下来一行包含N个正整数Ci，其中Ci表示第i只袜子的颜色，相同的颜色用相同的数字表示。再接下来M行，每行两个正整数L，R表示一个询问。

Output

包含M行，对于每个询问在一行中输出分数A/B表示从该询问的区间[L,R]中随机抽出两只袜子颜色相同的概率。若该概率为0则输出0/1，否则输出的A/B必须为最简分数。（详见样例）

Sample Input

6 4

1 2 3 3 3 2

2 6

1 3

3 5

1 6

Sample Output

2/5

0/1

1/1

4/15

【样例解释】

询问1：共C(5,2)=10种可能，其中抽出两个2有1种可能，抽出两个3有3种可能，概率为(1+3)/10=4/10=2/5。

询问2：共C(3,2)=3种可能，无法抽到颜色相同的袜子，概率为0/3=0/1。

询问3：共C(3,2)=3种可能，均为抽出两个3，概率为3/3=1/1。

注：上述C(a, b)表示组合数，组合数C(a, b)等价于在a个不同的物品中选取b个的选取方案数。

【数据规模和约定】

30%的数据中 N,M ≤ 5000；

60%的数据中 N,M ≤ 25000；

100%的数据中 N,M ≤ 50000，1 ≤ L < R ≤ N，Ci ≤ N。

嗯刚刚学习了分块
顺便讲了一下莫队
感觉真的很神奇啊%%%

莫队就是处理不带修改的区间查询问题的离线算法
据说是由曼哈顿距离最小生成树而得的
简单来说
就是如果可以将询问按照端点排序
并且可以由此区间端点+-1得到下一个区间
那么结合分块的单点修改和区间查询
可以将左端点排序后分块
再将块内按照右端点进行排序
那么这个复杂度应该是O（n*根号n）
在5*10^4的范围以内都可以

然后小z的袜子这道题
不仅要用莫队
还要用到数学知识···
就是一个公式
取到两只同颜色的袜子的方案数: Σ（i从1到n）cnt[i]∗(cnt[i]−1)(cnt[i]>=2) (cnt[i]表示颜色为i的袜子有几只)

取的总方案数：(r-l+1)*(r-l)
上下一除一约分就好啦

但是代码太复杂不带了自己打了QAQ

所以找了找终于找到一个不鬼畜的比较好懂的代码【虽然有点长】
然后自己加了些注释
如下：

#define ll long long
#include<algorithm>
#include<iostream>
#include<cstdio>
#define N 100010
#include<cmath>
using namespace std;
struct node
{
    ll id,l,r,z,m;//z就是分子，m就是分母 
}ask[N];
ll num[200010],n,block;
ll a[N],q,tot;
bool comp(node aa,node bb)
{
    if(aa.l/block==bb.l/block)
    {
        return aa.r<bb.r;
    }
    return aa.l/block<bb.l/block;
}
bool cmp(node aa,node bb)
{
    return aa.id<bb.id;
}
void add(ll nu)//加1
{
    ll pa=num[nu];
    num[nu]++;
    ll no=num[nu];
    tot-=pa*(pa-1)/2;//修改，减去之前加的那个方案数
    tot+=no*(no-1)/2;//加上加一以后的方案数
}
void pop(ll nu)//减1
{
    ll pa=num[nu];
    num[nu]--;
    ll no=num[nu];
    tot-=pa*(pa-1)/2; 
    tot+=no*(no-1)/2;
}
ll gcd(ll a,ll b)
{
    return b==0?a:gcd(b,a%b);
}
int main()
{
    scanf("%lld",&n);
    scanf("%lld",&q);
    block=sqrt(n);
    for(ll i=1;i<=n;++i) scanf("%lld",&a[i]);
    for(ll i=1;i<=q;++i)
    {
        scanf("%lld%lld",&ask[i].l,&ask[i].r);
        ask[i].id=i;
    }
    sort(ask+1,ask+q+1,comp);
    for(ll i=ask[1].l;i<=ask[1].r;++i) add(a[i]);//先朴素处理第一个区间 
    ask[1].z=tot;
    ll s=(ask[1].r-ask[1].l+1);
    ask[1].m=s*(s-1)/2;
    if(ask[1].z==0) ask[1].m=1;
    else 
    {
        ll it=gcd(ask[1].z,ask[1].m);
        ask[1].z/=it;
        ask[1].m/=it;
    }
    ll L=ask[1].l,R=ask[1].r;
    for(ll i=2;i<=q;++i)
    {
        while(R<ask[i].r) add(a[++R]);
        while(R>ask[i].r) pop(a[R--]);
        while(L>ask[i].l) add(a[--L]);
        while(L<ask[i].l) pop(a[L++]);
        ask[i].z=tot;
        ll s=(R-L+1);
        ask[i].m=s*(s-1)/2;
        if(ask[i].z==0) ask[i].m=1;
        else 
        {
            ll it=gcd(ask[i].z,ask[i].m);
            ask[i].z/=it;
            ask[i].m/=it;//约分 
        }
    }
    sort(ask+1,ask+q+1,cmp);
    for(ll i=1;i<=q;++i)
    printf("%lld/%lld\n",ask[i].z,ask[i].m);
    return 0;
}