专栏首页数据结构与算法ZR#330. 【18 提高 3】矿石(容斥)

ZR#330. 【18 提高 3】矿石(容斥)

题意

题目链接

Sol

挺显然的,首先对每个矿排序

那么答案就是$2^x - 2^y$

$x$表示能覆盖到它的区间,$y$表示的是能覆盖到它且覆盖到上一个的区间

第一个可以差分维护

第二个直接vector暴力插入扫就行,

时间复杂度:$O(nlogn)$

#include<cstdio>
#include<algorithm>
#include<bitset>
#include<vector>
#define Pair pair<int, ull>
#define MP(x, y) make_pair(x, y)
#define fi first
#define se second
#define ull unsigned long long 
#define LL long long
#define int long long  
using namespace std;
const int MAXN = 3 * 1e6 + 10, INF = 1e9 + 7, mod = 998244353;
inline int read() {
    char c = getchar(); int x = 0, f = 1;
    while(c < '0' || c > '9') {if(c == '-') f = -1; c = getchar();}
    while(c >= '0' && c <= '9') x = x * 10 + c - '0', c = getchar();
    return x * f;
}
int N, M;
int l[MAXN], r[MAXN], a[MAXN], date[MAXN], cnt = 0, sum[MAXN], num[MAXN];
vector<int> v[MAXN];
int fp(int a, int p) {
    int base = 1;
    while(p) {
        if(p & 1) base = (base * a) % mod;
        a = (a * a) % mod; p >>= 1;
    }
    return base % mod;
}
main() {
    N = read(); M = read();
    for(int i = 1; i <= N; i++) l[i] = read(), r[i] = read(), date[++cnt] = l[i], date[++cnt] = r[i];
    for(int i = 1; i <= M; i++) a[i] = read(), date[++cnt] = a[i];
    sort(a + 1, a + M + 1);
    sort(date + 1, date + cnt + 1);
    cnt = unique(date + 1, date + cnt + 1) - date - 1;
    for(int i = 1; i <= N; i++) {
        l[i] = lower_bound(date + 1, date + cnt + 1, l[i]) - date;
        r[i] = lower_bound(date + 1, date + cnt + 1, r[i]) - date;
        sum[l[i]]++; sum[r[i] + 1]--; 
        v[l[i]].push_back(r[i]);
    }
    for(int i = 1; i <= M; i++)  a[i] = lower_bound(date + 1, date + cnt + 1, a[i]) - date;
    for(int i = 1; i <= cnt; i++) 
        sum[i] += sum[i - 1];
    
   // for(int i = 1; i <= cnt; i++)
    //    printf("%d ", num[i]); puts("");
    int ans = 0;
    for(int i = 1; i <= M; i++) {
           int base = 0;
        for(int j = a[i - 1] + 1; j <= a[i]; j++) {
            for(int k = 0; k < v[j].size(); k++) {
                if(v[j][k] >= a[i]) base++;
            }
        }
        ans = (ans + fp(2, sum[a[i]]) - fp(2, sum[a[i]] - base) + mod) % mod;
    }
    printf("%lld\n", (ans + mod) % mod);
    return 0;
}
/*
3 2
7 11
1 5
3 8
4
7
*/

本文参与腾讯云自媒体分享计划,欢迎正在阅读的你也加入,一起分享。

我来说两句

0 条评论
登录 后参与评论

相关文章

  • cf1027F. Session in BSU(并查集 匈牙利)

    $n$个人,每个人可以在第$a_i$天或第$b_i$,一天最多考一场试,问在最优的情况下,最晚什么时候结束

    attack
  • BZOJ3693: 圆桌会议(Hall定理 线段树)

    我的思路:对于区间分两种情况讨论,一种是完全包含,另一种是部分包含。 第一种情况非常好判断,至于计算对于一个区间[l, r]的$\sum a[i]$就可以了,但...

    attack
  • agc016D - XOR Replace(图论 智商)

    不难看出,我们把所有数$xor$起来的数替换掉之后再次$xor$,得到的一定是被替换掉的数。

    attack
  • cf1027F. Session in BSU(并查集 匈牙利)

    $n$个人,每个人可以在第$a_i$天或第$b_i$,一天最多考一场试,问在最优的情况下,最晚什么时候结束

    attack
  • BZOJ3693: 圆桌会议(Hall定理 线段树)

    我的思路:对于区间分两种情况讨论,一种是完全包含,另一种是部分包含。 第一种情况非常好判断,至于计算对于一个区间[l, r]的$\sum a[i]$就可以了,但...

    attack
  • Typecho文章代码高亮功能

    Typecho是一款由国人开发的博客程序,它的特点是简洁小巧,扩展性强,并且内置支持Markdown语法写作,因而很受技术博客作者的欢迎。但是默认的Typec...

    HCG_Sky
  • LOJ#6281. 数列分块入门 5

    内存限制:256 MiB时间限制:500 ms标准输入输出 题目类型:传统评测方式:文本比较 上传者: hzwer 提交提交记录统计讨论 1 测试数据 题...

    attack
  • 1066 图像过滤 (15 分)

    版权声明:本文为博主原创文章,遵循 CC 4.0 BY-SA 版权协议,转载请附上原文出处链接和本声明。 ...

    韩旭051
  • 图解面试题:如何找到喜欢的电影?

    某电影平台(类似豆瓣、猫眼电影)用3个表来记录电影信息。“电影表”中是电影编号、电影名称、电影描述信息。

    猴子聊数据分析
  • 调整数组顺序使奇数位于偶数前面

    输入一个整数数组,实现一个函数来调整该数组中数字的顺序,使得所有奇数位于数组的前半部分,所有偶数位于数组的后半部分。

    编程珠玑

扫码关注云+社区

领取腾讯云代金券