题解-洛谷P4229 某位歌姬的故事
题面
\(T\) 组测试数据。有 \(n\) 个音节,每个音节 \(h_i\in[1,A]\),还有 \(m\) 个限制 \((l_i,r_i,g_i)\) 表示 \(\max_{k=l_i}^{r_i}h_k=g_i\)。求满足条件的 \(h_i\) 的方案数膜 \(998244353\)。
数据范围:\(1\le T\le 20\),\(1\le l_i\le r_i\le n\le 9\cdot 10^8\),\(1\le g_i\le A\le 9\cdot 10^8\),\(1\le m\le 500\)。
蒟蒻语
这是一眼题,蒟蒻 \(7\) 个小时就秒掉了。
题解做法时间复杂度 \(\Theta(Tm\log m)\)。
蒟蒻解
转化条件:\(\forall k\in[l_i,r_i],h_k\le g_i\) 并且 \(\exists k\in [l_i,r_i],h_k\ge g_i\)。
可以离散化,所以下文中可以把 \(n\) 看成与 \(m\) 相等的范围。
很明显对于 \(g_i<g_j\) 如果一个区间满足了 \(i\) 条件就无法满足 \(j\) 条件了,所以可以按 \(g\) 从小到大排序条件,每次处理完以后把区间删除;但是又会发现对于 \(g\) 相等的区间情况很复杂,所以可以排序后对每个 \(g\) 分开处理覆盖部分答案求积。删除区间可以用并查集(参考),然后最后对于每个没删除的音节答案乘以 \(A\) 即可。
然后只需假设 \(g\) 一定,并且区间覆盖全部音节。这时有多个棘手的条件,根据条件的形式容易想到子集反演。
设 \(p=(\forall k\in[l_i,r_i],h_k\le g_i)\),\(q=(\exists k\in [l_i,r_i],h_k\ge g_i)\),所以:
\]
这个子集反演是可以用 \(\rm dp\) 维护的,\(f(i,j)\) 表示 \(\forall k\in[i,j),h_k<g_i\) 前 \(i\) 个音节的方案数。
\]
\]
假设有个以 \(i\) 为左端点的条件 \([i,t)\),可以用 \(\rm dp\) 维护反演老套路打个负标记并用类似 \(\rm min-max\) 卷积的方法:
\]
操作比看起来简单:把 \(j\ge t\) 的 \(f(i,j)\) 直接置为 \(0\),然后 \(f(i,t)\leftarrow\sum_{j<t} f(i,j)=\sum f(i,j)\)。
然后是一波骚操作:把 \(i\) 这维滚掉(用离散化和并查集),然后线段树维护 \(j\) 这一维(啊!\(\rm dp\) 沦陷了!)。
最后还有个小细节:要把 \(g\) 相等但是 \([l_i,r_i]\in[l_j,r_j]\) 的 \(j\) 条件不放到 \(\rm dp\) 里,最后删除区间时统计。
时间复杂度 \(\Theta(Tm\log m)\)。 蒟蒻的题解都写成这样了,有人看得懂就怪了。
代码
#include <bits/stdc++.h>
using namespace std;
//Start
typedef long long ll;
typedef double db;
#define mp(a,b) make_pair((a),(b))
#define x first
#define y second
#define be(a) (a).begin()
#define en(a) (a).end()
#define sz(a) int((a).size())
#define pb(a) push_back(a)
#define R(i,a,b) for(int i=(a),I=(b);i<I;i++)
#define L(i,a,b) for(int i=(b)-1,I=(a)-1;i>I;i--)
const int iinf=0x3f3f3f3f;
const ll linf=0x3f3f3f3f3f3f3f3f;
//Data
const int Q=501,N=Q<<1,mod=998244353;
int Pow(int a,int x){
if(!a) return 0; int res=1;
for(;x;a=1ll*a*a%mod,x>>=1)if(x&1) res=1ll*a*res%mod;
return res;
}
int n,m,h,l[Q],r[Q],g[Q],dn,d[N],te[N+1];
int find(int u){return u==te[u]?u:(te[u]=find(te[u]));}
map<int,vector<int>> mp;
//SegmentTree
const int T=N<<2; int val[T],ml[T];
#define mid ((l+r)>>1)
void clear(int k=0,int l=0,int r=dn){
val[k]=0,ml[k]=1; if(r-l==1) return;
clear(k*2+1,l,mid),clear(k*2+2,mid,r);
}
void pushup(int k){(val[k]=val[k*2+1]+val[k*2+2])%=mod;}
void pushmul(int k,int v){val[k]=1ll*val[k]*v%mod,ml[k]=1ll*ml[k]*v%mod;}
void pushdown(int k){if(ml[k]^1) pushmul(k*2+1,ml[k]),pushmul(k*2+2,ml[k]),ml[k]=1;}
void add(int x,int v,int k=0,int l=0,int r=dn){
if(r<=x||x+1<=l) return; if(r-l==1) return void((val[k]+=v)%=mod);
pushdown(k),add(x,v,k*2+1,l,mid),add(x,v,k*2+2,mid,r),pushup(k);
}
void mul(int x,int y,int v,int k=0,int l=0,int r=dn){
if(r<=x||y<=l) return; if(x<=l&&r<=y) return pushmul(k,v);
pushdown(k),mul(x,y,v,k*2+1,l,mid),mul(x,y,v,k*2+2,mid,r),pushup(k);
}
//DP
int dp(int x,vector<int> a){
mul(0,dn,0),add(0,1); vector<int> b;
for(int i:a) l[i]=find(l[i]),r[i]=find(r[i]);
sort(be(a),en(a),[&](int p,int q){return l[p]!=l[q]?l[p]<l[q]:r[p]<r[q];});
for(int i:a){while(sz(b)&&r[b.back()]>=r[i]) b.pop_back();b.pb(i);}
R(t,0,sz(b)){
int i=find(l[b[t]]),j=find(r[b[t]]); if(i>=j) return 0;
int ni=j;if(t+1<sz(b)) ni=min(ni,find(l[b[t+1]]));
mul(j,dn,0),add(j,(mod-val[0])%mod);
for(int k=i;k<ni;k=te[k]=find(k+1))
mul(0,k+1,Pow(x,d[k+1]-d[k])),mul(k+1,dn,Pow(x-1,d[k+1]-d[k]));
}
int res=val[0];
for(int i:a) for(int k=find(l[i]);k<find(r[i]);k=te[k]=find(k+1)) res=1ll*res*Pow(x,d[k+1]-d[k])%mod;
return res;
}
//Main
int Main(){
cin>>n>>m>>h,dn=0,d[dn++]=0,d[dn++]=n,mp.clear();
R(i,0,m) cin>>l[i]>>r[i]>>g[i],--l[i],d[dn++]=l[i],d[dn++]=r[i];
sort(d,d+dn),dn=unique(d,d+dn)-d,clear();
R(i,0,m) l[i]=lower_bound(d,d+dn,l[i])-d,r[i]=lower_bound(d,d+dn,r[i])-d,mp[g[i]].pb(i);
iota(te,te+dn+1,0); int res=1;
for(auto u:mp) res=1ll*res*dp(u.x,u.y)%mod;
for(int k=find(0);k<dn-1;k=te[k]=find(k+1)) res=1ll*res*Pow(h,d[k+1]-d[k])%mod;
return res;
}
int main(){
ios::sync_with_stdio(0);
cin.tie(0),cout.tie(0);
int ta; cin>>ta; while(ta--) cout<<Main()<<'\n';
return 0;
}
祝大家学习愉快!
题解-洛谷P4229 某位歌姬的故事的更多相关文章
- Loj #2331. 「清华集训 2017」某位歌姬的故事
Loj #2331. 「清华集训 2017」某位歌姬的故事 IA 是一名会唱歌的女孩子. IOI2018 就要来了,IA 决定给参赛选手们写一首歌,以表达美好的祝愿.这首歌一共有 \(n\) 个音符, ...
- 题解-洛谷P4859 已经没有什么好害怕的了
洛谷P4859 已经没有什么好害怕的了 给定 \(n\) 和 \(k\),\(n\) 个糖果能量 \(a_i\) 和 \(n\) 个药片能量 \(b_i\),每个 \(a_i\) 和 \(b_i\) ...
- 题解 洛谷P5018【对称二叉树】(noip2018T4)
\(noip2018\) \(T4\)题解 其实呢,我是觉得这题比\(T3\)水到不知道哪里去了 毕竟我比较菜,不大会\(dp\) 好了开始讲正事 这题其实考察的其实就是选手对D(大)F(法)S(师) ...
- 题解 洛谷P1562 【还是N皇后】
原题:洛谷P1562 这个题的原理和8皇后的原理是一模一样的,就是必须要用n个皇后把每一个行填满,同时满足每一列,每一行,每一条对角线只有一个棋子.但如果按照原来的方法暴打的话只有60分(优化亲测无效 ...
- 题解 洛谷 P3396 【哈希冲突】(根号分治)
根号分治 前言 本题是一道讲解根号分治思想的论文题(然鹅我并没有找到论文),正 如论文中所说,根号算法--不仅是分块,根号分治利用的思想和分块像 似却又不同,某一篇洛谷日报中说过,分块算法实质上是一种 ...
- 题解-洛谷P5410 【模板】扩展 KMP(Z 函数)
题面 洛谷P5410 [模板]扩展 KMP(Z 函数) 给定两个字符串 \(a,b\),要求出两个数组:\(b\) 的 \(z\) 函数数组 \(z\).\(b\) 与 \(a\) 的每一个后缀的 L ...
- 题解-洛谷P4724 【模板】三维凸包
洛谷P4724 [模板]三维凸包 给出空间中 \(n\) 个点 \(p_i\),求凸包表面积. 数据范围:\(1\le n\le 2000\). 这篇题解因为是世界上最逊的人写的,所以也会有求凸包体积 ...
- 题解-洛谷P5217 贫穷
洛谷P5217 贫穷 给定长度为 \(n\) 的初始文本 \(s\),有 \(m\) 个如下操作: \(\texttt{I x c}\),在第 \(x\) 个字母后面插入一个 \(c\). \(\te ...
- 题解 洛谷 P2010 【回文日期】
By:Soroak 洛谷博客 知识点:模拟+暴力枚举 思路:题目中有提到闰年然后很多人就认为,闰年是需要判断的其实,含有2月29号的回文串,前四位是一个闰年那么我们就可以直接进行暴力枚举 一些小细节: ...
随机推荐
- 聊一聊sockmap 以及ebpf 实例演示
eBPF实质上是一个内核注入技术 用户态可以用C来写运行的代码,再通过一个Clang&LLVM的编译器将C代码编译成BPF目标码: 用户态通过系统调用bpf()将BPF目标码注入到内核当中,并 ...
- 加解密 C语言实现
1.加密的基本原理 加密分为对称加密和非对称加密,对称加密就是加密方和解密放用同一个密钥. 加密是分组加密,即将明文数据分成多个密钥大小的块,依次和密钥运算,输出密文. padding,由于加密需要分 ...
- git命令学习之clone指定分支代码
今天要拉取一个项目,但是是一个指定分支,本来我以为直接git clone就行,但是发现好像不能,报错: Cloning into 'lecture'...fatal: unable to update ...
- Kernel RBD的QOS配置方案
前言 关于qos的讨论有很多,ceph内部也正在实现着一整套的基于dmclock的qos的方案,这个不是本篇的内容,之前在社区的邮件列表看过有研发在聊qos的相关的实现的,当时一个研发就提出了在使用k ...
- 部署sftp服务
部署sftp服务有风险,可能造成ssh无法连接到服务器,因此写个脚本定时覆盖一下,保证ssh可以正常使用. 创建数据目录并赋权,创建账号密码,修改ssh文件. * mkdir /sftp groupa ...
- GreenDao增删改查
3.GreenDao增删改查 (1)插入 常用API //这是最简单的插入语句,新增一行数据,返回值为行号 public long insert(T entity) //传递一个数组,新增多行数据 p ...
- 循序渐进VUE+Element 前端应用开发(27)--- 数据表的动态表单设计和数据存储
在我们一些系统里面,有时候会需要一些让用户自定义的数据信息,一般这些可以使用扩展JSON进行存储,不过每个业务表的显示项目可能不一样,因此需要根据不同的表单进行设计,然后进行对应的数据存储.本篇随笔结 ...
- BUUCTF 不一样的flag writeup
感谢BUUCTF提供的学习平台 https://buuoj.cn 题目:不一样的flag 工具:x64dbg 这是一道内存的迷宫题,迷宫是402000处的字符串 根据经验,这应该(a行*b列)的字符, ...
- webug第三关:你看到了什么?
第三关:你看到了什么? 右键源码 扫描到test目录
- Linux中配置环境变量
Linux中环境变量的搭建(推荐用法) 第一步:进入到/etc/profile.d文件夹下 cd /etc/profile.d 第二步:创建并编辑一个my_env.sh文件 vim my_env.sh ...