[CF799F]Beautiful fountains rows

3 篇文章 0 订阅
3 篇文章 0 订阅

Description

在一个长度为m的数轴上,有n种球,每种球会出现在区间[l,r]中。
一个合法的区间满足:这个区间里有球,并且每种出现过的球都出现了奇数次
求所有合法的区间的长度之和。
n,m<=2*1e5

Solution

讲课时选了这道题然而似乎讲错题意了
直接数据结构可能很麻烦而且常数大,我们考虑一种取巧的做法:随机化。
给每个求随机一个权值key,那么一个区间合法的条件就是:这个区间里面每个位置的异或和的异或和=所有出现过的球的权值的异或和。
前一个东西可以用前缀和解决。
后一个东西可以写成所有右端点在l右边的球^所有左端点在r左边的球^所有球。
这样也可以用两个前缀和解决。
于是枚举右端点直接求所有的合法左端点即可。用hash/map实现。
注意空区间是会被我们计算的,要减去。
复杂度O(n)/O(n log n)

Ps:上网找了一下发现基本单哈希是过不去的

Code

#include <map>
#include <cstdio>
#include <cstring>
#include <algorithm>
#define fo(i,a,b) for(int i=a;i<=b;i++)
#define fd(i,a,b) for(int i=a;i>=b;i--)
using namespace std;

typedef long long ll;

int read() {
    char ch;
    for(ch=getchar();ch<'0'||ch>'9';ch=getchar());
    int x=ch-'0';
    for(ch=getchar();ch>='0'&&ch<='9';ch=getchar()) x=x*10+ch-'0';
    return x;
}

const int N=2*1e5+5;

int n,m,l,r,num[N];
ll pre_1[N],sl_1[N],sr_1[N],all_1,pre_2[N],sl_2[N],sr_2[N],all_2;
map<ll,ll> cnt,sum;
ll ans;

ll calc(int x) {
    ll res=(ll)x*(x+1)/2*(x+1);
    res-=(ll)x*(x+1)*(2*x+1)/6;
    return res;
}

int main() {
    n=read();m=read();srand(20180316);
    fo(i,1,n) {
        l=read();r=read();
        ll key_1=(ll)rand()*rand()*rand()*rand();
        pre_1[l]^=key_1;pre_1[r+1]^=key_1;
        sl_1[r]^=key_1;sr_1[l]^=key_1;all_1^=key_1;
        ll key_2=(ll)rand()*rand()*rand()*rand();
        pre_2[l]^=key_2;pre_2[r+1]^=key_2;
        sl_2[r]^=key_2;sr_2[l]^=key_2;all_2^=key_2;
        num[l]++;num[r+1]--;
    }
    fo(i,1,m) num[i]+=num[i-1];
    fo(i,1,m) pre_1[i]^=pre_1[i-1],pre_2[i]^=pre_2[i-1];
    fo(i,1,m) pre_1[i]^=pre_1[i-1],pre_2[i]^=pre_2[i-1];
    fo(i,1,m) sr_1[i]^=sr_1[i-1],sr_2[i]^=sr_2[i-1];
    fd(i,m,1) sl_1[i]^=sl_1[i+1],sl_2[i]^=sl_2[i+1];
    cnt[sl_1[1]^sl_2[1]]++;
    fo(i,1,m) {
        ll now=all_1^pre_1[i]^sr_1[i]^all_2^pre_2[i]^sr_2[i];
        ans+=(ll)i*cnt[now];ans-=sum[now];
        cnt[pre_1[i]^sl_1[i+1]^pre_2[i]^sl_2[i+1]]++;
        sum[pre_1[i]^sl_1[i+1]^pre_2[i]^sl_2[i+1]]+=i;
    }
    for(int l=1,r=1;l<=m;l=r+1) {
        r=l;
        while (r<=m&&!num[r]) r++;
        ans-=calc(r-l);
    }
    printf("%I64d\n",ans);
    return 0;
}
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值