bzoj2038(小z的袜子——莫队)

最新推荐文章于 2018-08-24 12:13:12 发布

原创最新推荐文章于 2018-08-24 12:13:12 发布 · 311 阅读

0 ·

CC 4.0 BY-SA版权

莫队专栏收录该内容

2 篇文章

订阅专栏

本文详细介绍了莫队算法的应用及实现细节，通过一道具体的题目来展示莫队算法的基本写法，包括分块策略、区间更新操作等，并强调了在处理大规模数据时避免整数溢出的重要性。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

莫队的模板题。莫队的注意点记在蓝书上了。

本题关键的是莫队的基本写法需要记清楚

#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
using namespace std;
typedef long long ll;
const int N=60005;
int n,m,c[N];
int pos[N];
ll ans,ml,mr,sum[N];
struct aa//包含询问信息的结构体
{
    int l,r,id;
    ll ansl,ansr;
}qu[N];
ll gcd(ll a,ll b) 
{
    return  b==0?a:gcd(b,a%b);
}
bool cmp1(aa a,aa b)
{
    if (pos[a.l]!=pos[b.l]) return pos[a.l]<pos[b.l];
    return a.r<b.r;
}
bool cmp2(aa a,aa b)
{
    return a.id<b.id;
}
void del(int i)//删除
{
    ans-=sum[c[i]]*(sum[c[i]]-1)/2;
    sum[c[i]]--;
    ans+=sum[c[i]]*(sum[c[i]]-1)/2;
}
void ins(int i)//插入
{
    ans-=sum[c[i]]*(sum[c[i]]-1)/2;
    sum[c[i]]++;
    ans+=sum[c[i]]*(sum[c[i]]-1)/2;
}
void solve()
{
    ml=1;mr=0;//区间初始位置
    for (int i=1;i<=m;i++)
    {
        for (;mr<qu[i].r;mr++) ins(mr+1);
        for (;mr>qu[i].r;mr--) del(mr);
        for (;ml<qu[i].l;ml++) del(ml);
        for (;ml>qu[i].l;ml--) ins(ml-1);//莫队找区间
        if (qu[i].l==qu[i].r) {qu[i].ansl=1;continue;}
        if (ans==0) {qu[i].ansr=1;continue;}
         
        qu[i].ansr=(ll)(qu[i].r-qu[i].l+1)*(qu[i].r-qu[i].l)/2;//这里必须要加ll
        qu[i].ansl=ans;
         
         
        ll d=gcd(qu[i].ansl,qu[i].ansr);
        qu[i].ansl=qu[i].ansl/d;
        qu[i].ansr=qu[i].ansr/d; 
    }
}
int main()
{
    scanf("%d%d",&n,&m);
    int block=int(sqrt(n*1.0));
    for (int i=1;i<=n;i++) pos[i]=(i-1)/block+1//分块
     
    for (int i=1;i<=n;i++) scanf("%d",&c[i]);
    for (int i=1;i<=m;i++) scanf("%d%d",&qu[i].l,&qu[i].r),qu[i].id=i;
     
    sort(qu+1,qu+m+1,cmp1);
    solve();
    sort(qu+1,qu+m+1,cmp2);
     
    for (int i=1;i<=m;i++) printf("%lld/%lld\n",qu[i].ansl,qu[i].ansr);
    return 0; 
}

总结

关于这道题，有一个地方必须要开ll，注意这里的n可能达到50000，而50000^2>maxint，会爆int的。

刚开始想到有可能爆int错了，但是没有考虑清楚，找到爆点。实际上。以后关于这些乘法呀加法呀要更加注意爆int。！！！！