题面

洛谷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)\),所以:

\[p(S)\&q(S)=\sum_{T\in S}(-1)^{|T|}p(S-T)\neg q(T),\neg q(T)=(\forall k\in[l_i,r_i],h_k<g_i)
\]

这个子集反演是可以用 \(\rm dp\) 维护的,\(f(i,j)\) 表示 \(\forall k\in[i,j),h_k<g_i\) 前 \(i\) 个音节的方案数。

\[\forall j\le i,f(i,j)\leftarrow g\cdot f(i-1,j)
\]
\[\forall j> i,f(i,j)\leftarrow (g-1)\cdot f(i-1,j)
\]

假设有个以 \(i\) 为左端点的条件 \([i,t)\),可以用 \(\rm dp\) 维护反演老套路打个负标记并用类似 \(\rm min-max\) 卷积的方法:

\[\forall j,f(i,\max(j,t))\leftarrow f(i,\max(j,t))-f(i,j)
\]

操作比看起来简单:把 \(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 某位歌姬的故事的更多相关文章

  1. Loj #2331. 「清华集训 2017」某位歌姬的故事

    Loj #2331. 「清华集训 2017」某位歌姬的故事 IA 是一名会唱歌的女孩子. IOI2018 就要来了,IA 决定给参赛选手们写一首歌,以表达美好的祝愿.这首歌一共有 \(n\) 个音符, ...

  2. 题解-洛谷P4859 已经没有什么好害怕的了

    洛谷P4859 已经没有什么好害怕的了 给定 \(n\) 和 \(k\),\(n\) 个糖果能量 \(a_i\) 和 \(n\) 个药片能量 \(b_i\),每个 \(a_i\) 和 \(b_i\) ...

  3. 题解 洛谷P5018【对称二叉树】(noip2018T4)

    \(noip2018\) \(T4\)题解 其实呢,我是觉得这题比\(T3\)水到不知道哪里去了 毕竟我比较菜,不大会\(dp\) 好了开始讲正事 这题其实考察的其实就是选手对D(大)F(法)S(师) ...

  4. 题解 洛谷P1562 【还是N皇后】

    原题:洛谷P1562 这个题的原理和8皇后的原理是一模一样的,就是必须要用n个皇后把每一个行填满,同时满足每一列,每一行,每一条对角线只有一个棋子.但如果按照原来的方法暴打的话只有60分(优化亲测无效 ...

  5. 题解 洛谷 P3396 【哈希冲突】(根号分治)

    根号分治 前言 本题是一道讲解根号分治思想的论文题(然鹅我并没有找到论文),正 如论文中所说,根号算法--不仅是分块,根号分治利用的思想和分块像 似却又不同,某一篇洛谷日报中说过,分块算法实质上是一种 ...

  6. 题解-洛谷P5410 【模板】扩展 KMP(Z 函数)

    题面 洛谷P5410 [模板]扩展 KMP(Z 函数) 给定两个字符串 \(a,b\),要求出两个数组:\(b\) 的 \(z\) 函数数组 \(z\).\(b\) 与 \(a\) 的每一个后缀的 L ...

  7. 题解-洛谷P4724 【模板】三维凸包

    洛谷P4724 [模板]三维凸包 给出空间中 \(n\) 个点 \(p_i\),求凸包表面积. 数据范围:\(1\le n\le 2000\). 这篇题解因为是世界上最逊的人写的,所以也会有求凸包体积 ...

  8. 题解-洛谷P5217 贫穷

    洛谷P5217 贫穷 给定长度为 \(n\) 的初始文本 \(s\),有 \(m\) 个如下操作: \(\texttt{I x c}\),在第 \(x\) 个字母后面插入一个 \(c\). \(\te ...

  9. 题解 洛谷 P2010 【回文日期】

    By:Soroak 洛谷博客 知识点:模拟+暴力枚举 思路:题目中有提到闰年然后很多人就认为,闰年是需要判断的其实,含有2月29号的回文串,前四位是一个闰年那么我们就可以直接进行暴力枚举 一些小细节: ...

随机推荐

  1. Scrum转型(一) 为什么敏捷和Scrum

    1.1 为什么敏捷 由于传统的瀑布模型管理方法无法满足现代某些软件产品开发过程的特点,我们需要使用敏捷的方法(例如,Scrum是一个让我们关注于在短时间里交付高质量商业价值的敏捷框架). 需求频繁变动 ...

  2. 储存与RAID--独立磁盘阵列

    存储:  专门用来插硬盘的机器,作用是增加插口,可以多插硬盘. 这种有策略保证硬盘坏了,数据不丢.而本地磁盘坏了,会导致数据丢失,故一般操作系统放在本地磁盘.而数据放在存储盘. 存储里依然有:cpu( ...

  3. SQL SERVER数据库使用过程中系统提示死锁处理办法

    马上双节(国庆节.中秋节)了,这篇文章是双节前的最后一篇,祈祷过节期间,数据库稳定运行,服务器正常发挥.祝大家假期愉快!!!! 任何的数据库都会出现死锁的情况,特别是一些大型的复杂业务,数据库架构的设 ...

  4. tp5获取当前域名

    //头部引入 use think\Request; //获取当前域名 $request = Request::instance(); $domain = $request->domain();

  5. Android ALSPS驱动分析

    一.alsps的初始化函数和重要结构体 epl2182_init // Epl2182.c (kernel-3.10\drivers\misc\mediatek\alsps\epl2182-new) ...

  6. EasyRecovery——信息时代的“后悔药”

    前几日,小编在豆瓣潜水的时候看到这么一个帖子,说是一对小情侣吵架,女方一气之下把男方的博士论文和资料全删了,求一个办法让男友消气. 站在吃瓜的角度,小编和广大群众看法一致,希望两人直接分手,放男方一条 ...

  7. 巧妙利用Camtasia制作网课

    随着互联网的快速发展,网络学习变得非常流行.这种躺在床上就可以获取知识的方法让大家渐渐地都喜欢上了学习,那么我们是否想要了解一下网课的幕后制作呢. 今天我给大家带来的便是巧妙利用Camtasia进行网 ...

  8. Improving Commonsense Question Answering by Graph-based Iterative Retrieval over Multiple Knowledge Sources —— 基于多知识库迭代检索的常识问答系统

    基于多知识库迭代检索的问答系统 论文地址 背景 常识问答任务需要引入外部知识来帮助模型更好地理解自然语言问题,现有的解决方案大都采用两阶段框架: 第一阶段 -- 从广泛的知识来源中找到与给定问题相关的 ...

  9. Session 与 sql 会话,mysql 权限设置,mybatis 逆向工程

    Session 与 Sql 会话注意点: 通过 sqlSessionFactoty 工厂建立的与sql的会话,在进行相应的插入操作后,需要进行 commit 操作才会让数据库执行插入更新操作.如何主键 ...

  10. 2017年第八届蓝桥杯【C++省赛B组】B、C、D、H 题解

    可能因为我使用暴力思维比较少,这场感觉难度不低. B. 等差素数列 #暴力 #枚举 题意 类似:\(7,37,67,97,127,157\) 这样完全由素数组成的等差数列,叫等差素数数列. 上边的数列 ...