bzoj2038: [2009国家集训队]小Z的袜子(hose)（莫队+分块模板题）

转载于 2017-08-16 11:21:57 发布 · 228 阅读

ACM 专栏收录该内容

101 篇文章

订阅专栏

本文介绍了一种使用莫队算法和分块技术解决特定问题的方法。该问题无法通过线段树直接处理，因为需要对访问顺序进行预处理以减少不必要的遍历次数。文中详细解释了莫队算法的时间复杂度分析，并提供了完整的实现代码。

题目链接：http://www.lydsy.com/JudgeOnline/problem.php?id=2038

思路：

转自：http://blog.csdn.net/bossup/article/details/39236275

这题用线段树貌似处理不了，因为左区间和右区间合并后还要重算一遍。如果直接暴力遍历每个访问的话，时间复杂度最坏可达m*n，这时就需要莫队+分块来预处理访问顺序，来减小重复的遍历。

这题我还没推出求概率的公式，好弱，，但是精髓在于莫队时间复杂度的分析：

一、i与i+1在同一块内，r单调递增，所以r是O(n)的。由于有n^0.5块,所以这一部分时间复杂度是n^1.5。
二、i与i+1跨越一块，r最多变化n，由于有n^0.5块，所以这一部分时间复杂度是n^1.5
三、i与i+1在同一块内时变化不超过n^0.5，跨越一块也不会超过2*n^0.5，不妨看作是n^0.5。由于有n个数，所以时间复杂度是n^1.5.

代码:

#include <iostream>
#include <cstdio>
#include <algorithm>
#include <cmath>
#include <cstring>


using namespace std;
typedef long long ll;
const int maxn = 50000+5;
ll block[maxn],val[maxn],ans,c[maxn];
struct node
{
    ll l,r,fz,fm;
    int id;
}q[maxn];
bool cmp1(node a,node b)
{
    if(block[a.l]==block[b.l])
        return a.r<b.r;
    return a.l<b.l;
}
bool cmp2(node a,node b)
{
    return a.id<b.id;
}

void update(int x,int v)
{
    ans -= val[c[x]]*val[c[x]];
    val[c[x]] += v;
    ans += val[c[x]]*val[c[x]];
}

ll gcd(ll x,ll y)
{
    if(!y)
        return x;
    return gcd(y,x%y);
}

int main()
{
    int n,m;
    while(~scanf("%d%d",&n,&m))
    {
        memset(val,0,sizeof(val));
        for(int i = 1;i<=n;i++)
            scanf("%lld",&c[i]);
        int len = sqrt(1.0*n);
        for(int i = 1;i<=n;i++)
            block[i] = (i-1)/len+1;
        for(int i = 0;i<m;i++)
        {
            scanf("%lld%lld",&q[i].l,&q[i].r);
            q[i].id = i;
        }
        sort(q,q+m,cmp1);
        ll ql = 1,qr = 0;
        ans = 0;
        for(int i = 0;i<m;i++)
        {
            while(qr<q[i].r)
            {
                ++qr;
                update(qr,1);
            }
            while(qr>q[i].r)
            {
                update(qr,-1);
                --qr;
            }
            while(ql<q[i].l)
            {
                update(ql,-1);
                ++ql;
            }
            while(ql>q[i].l)
            {
                --ql;
                update(ql,1);
            }
            if(q[i].l==q[i].r)
            {
                q[i].fz = 0,q[i].fm  = 1;
            }
            else
            {
                q[i].fz = ans-(qr-ql+1);
                q[i].fm = (qr-ql+1)*(qr-ql);
                ll g = gcd(q[i].fm,q[i].fz);
                q[i].fm/=g;
                q[i].fz/=g;
            }
        }
        sort(q,q+m,cmp2);
        for(int i = 0;i<m;i++)
            printf("%lld/%lld\n",q[i].fz,q[i].fm);
    }
    return 0;
}