【学习笔记】CF930E Coins Exhibition

news/2024/11/18 3:27:07/

感觉像是之前做过的题的加强版😅

考虑容斥哪些区间不合法。直接处理比较困难,考虑将所有区间按右端点排序,并将端点离散化(将右端点 + 1 +1 +1,转化为左闭右开区间),设 d p i , j , k dp_{i,j,k} dpi,j,k表示只考虑前 i i i个区间,以及 [ 1 , j ) [1,j) [1,j)这段前缀,上一个选择的区间类型是 k ∈ [ 0 , 1 ] k\in [0,1] k[0,1]时的答案。转移如下:

  • d p i , j , k ← d p i − 1 , j , k dp_{i,j,k}\gets dp_{i-1,j,k} dpi,j,kdpi1,j,k
  • d p i , j , k ′ ← − d p i − 1 , l i , k × 1 2 r i − l i dp_{i,j,k'}\gets -dp_{i-1,l_i,k}\times \frac{1}{2^{r_i-l_i}} dpi,j,kdpi1,li,k×2rili1,条件: r i ≤ j r_i\le j rij,可以是相同类型的区间也可以是不同类型的区间
  • d p i , j , k ← − d p i − 1 , j , k dp_{i,j,k}\gets -dp_{i-1,j,k} dpi,j,kdpi1,j,k,条件: l i ≥ j l_i\ge j lij,且必须是相同类型区间
  • d p i , j , k ← − d p i − 1 , l i , k × 1 2 j − l i dp_{i,j,k}\gets -dp_{i-1,l_i,k}\times \frac{1}{2^{j-l_i}} dpi,j,kdpi1,li,k×2jli1,条件: l i < j < r i l_i<j< r_i li<j<ri,且必须是相同类型区间

最后答案要乘上 2 K 2^K 2K

显然,这些操作都可以用线段树去维护。

有没有更好的方法?

注意到,第三种转移加上第一种转移是将 ≤ l i \le l_i li D P DP DP值推平成 0 0 0,那么我们维护一个指针 p p p表示 [ 1 , p ] [1,p] [1,p]这段前缀的 D P DP DP值都是 0 0 0,如果 l i ≤ p l_i\le p lip那么什么都不做;否则我们暴力将指针移动到 l i l_i li,然后根据转移的范围在差分数组上打标记即可。

复杂度 O ( n log ⁡ n ) O(n\log n) O(nlogn)

remark \text{remark} remark 我低估了这道题的思维难度。。。(主要是后半部分)

#include<bits/stdc++.h>
#define ll long long
#define pb push_back
#define inf 0x3f3f3f3f
#define fi first
#define se second
using namespace std;
const int N=4e5+5;
const int mod=1e9+7;
int n,m,K,lsh[N<<1],cnt;
struct node{int l,r,t;bool operator <(const node &a)const{return r<a.r;}
}a[N];
ll fpow(ll x,ll y=mod-2){ll z(1);for(;y;y>>=1){if(y&1)z=z*x%mod;x=x*x%mod;}return z;
}
int get(int x){return lower_bound(lsh+1,lsh+1+cnt,x)-lsh;
}
int p[2];
ll c[2][N][2];
void add(ll &x,ll y){x=(x+y)%mod;
}
ll calc(int f,int x){return (c[f][x][0]+c[f][x][1]*fpow(mod+1>>1,lsh[x]))%mod;
}
int main(){ios::sync_with_stdio(false);cin.tie(0),cout.tie(0);cin>>K>>n>>m;for(int i=1;i<=n+m;i++){cin>>a[i].l>>a[i].r,a[i].r++,a[i].t=(i<=n);lsh[++cnt]=a[i].l,lsh[++cnt]=a[i].r;}sort(lsh+1,lsh+1+cnt),cnt=unique(lsh+1,lsh+1+cnt)-lsh-1;n+=m;for(int i=1;i<=n;i++)a[i].l=get(a[i].l),a[i].r=get(a[i].r);sort(a+1,a+1+n);c[0][1][0]=c[1][1][0]=1;ll res=1;for(int i=1;i<=n;i++){int l=a[i].l,r=a[i].r,f=a[i].t;if(p[f]>=l)continue;while(p[f]<l){add(c[f][p[f]+1][0],c[f][p[f]][0]);add(c[f][p[f]+1][1],c[f][p[f]][1]);p[f]++;}ll x=calc(f,l),y=-x*fpow(mod+1>>1,lsh[r]-lsh[l])%mod;add(res,y);add(c[f^1][r][0],y);add(c[f][r][0],y);add(c[f][l+1][1],-x*fpow(2,lsh[l]));add(c[f][r][1],x*fpow(2,lsh[l]));}res=res*fpow(2,K)%mod;cout<<(res+mod)%mod;
}

http://www.ppmy.cn/news/1188869.html

相关文章

redis的分布式锁

分布式锁是一种用在分布式系统中实现同步和互斥访问的机制。 1、分布式锁概念 满足分布式系统或者集群模式下&#xff0c;多进程可见并且互斥的锁。 分布式锁的核心思想就是让分布式集群中的节点都适用同一把锁&#xff0c;只要大家使用的是同一把锁&#xff0c;就能锁住线程…

dbeaver查看表,解决证书报错current license is non-compliant for [jdbc]

http://localhost:9200/_license { “license” : { “status” : “active”, “uid” : “b91ae0e0-b04d-4e20-8730-cf0bca7b2035”, “type” : “basic”, “issue_date” : “2023-02-22T14:33:27.648Z”, “issue_date_in_millis” : 1677076407648, “max_nodes” : 10…

随记-动态切换数据库链接

动态切换数据库链接 如果你有多个数据库链接&#xff0c;并且在启动服务之前没有写入相关的配置文件&#xff0c;这种情况下该如何使用动态的增加数据库链接并使用呢&#xff1f; 前置条件 这里假设你的数据库配置信息已经写入默认的数据库中&#xff08;这里暂且称为应用程序…

write、fwrite、fstream二进制I/O的效率对比

测试数据 写文件测试数据 注&#xff1a;每项数据为3次运行的结果取平局值&#xff0c;运行时长为10ms。 单次传入字节数(B)写入字节数&#xff08;MB&#xff09;writefwritefstream.write10.005 0.124 0.127 40.019 0.453 0.450 160.064 1.840 1.860 640.363 4.503 4.277 …

什么是 DevOps

DevOps是一套融合软件开发&#xff08;Dev&#xff09;和 IT 运营&#xff08;Ops&#xff09;的实践&#xff0c;旨在缩短应用程序开发周期并确保以高软件质量持续交付&#xff0c;通过采用 DevOps 实践&#xff0c;您可以帮助组织更可靠、更快速、更高效地交付软件。 什么是…

高效改名,文件夹名称替换:一键批量替换文件夹名中间部分内容

在我们的日常生活和工作中&#xff0c;经常需要处理大量的文件夹&#xff0c;其中有些文件夹名称可能包含我们需要替换的内容。但如果我们一个一个地手动修改文件夹名称&#xff0c;不仅耗时而且容易出错。为了解决这个问题&#xff0c;我们可以使用云炫文件管理器高效的文件夹…

罢工效应,韩电池业在美建厂面临挑战 | 百能云芯

美国汽车行业的罢工事件仍在持续&#xff0c;对于远在数千公里之外的韩国电池制造商来说&#xff0c;这引发了不小的担忧&#xff0c;他们担心生产成本会因此大幅上升。 据彭博资讯报道&#xff0c;韩国电池制造商LG新能源、SK On和三星SDI&#xff0c;已与美国的三家汽车制造巨…

Pod资源限制与探针

目录 资源限制 资源限制概述 资源单位 示例 示例1 示例2 探针 探针的三种规则 Probe支持三种检查方法 示例 示例1&#xff1a;exec方式 示例2&#xff1a;httpGet方式 示例3&#xff1a;tcpSocket方式 示例4&#xff1a;就绪检测 示例5&#xff1a;就绪检测2 …