Subsequence Count (线段树)
Time Limit: 1000 ms Memory Limit: 256 MB
Description
给定一个01串 $S_{1 \cdots n}$ 和 $Q$ 个操作。
操作有两种类型:
1、将 $[l, r]$ 区间的数取反(将其中的0变成1,1变成0)。
2、询问字符串 $S$ 的子串 $S_{l \cdots r}$ 有多少个不同的子序列。由于答案可能很大,请将答案对 $10^9 + 7$ 取模。
在数学中,某个序列的子序列是从最初序列通过去除某些元素但不破坏余下元素的相对位置(在前或在后)而形成的新序列。
Input
第一行包含两个整数 $N$ 和 $Q$ ,分别表示字符串长度和操作次数。
第二行包含一个字符串 $S$ 。
接下来 $Q$ 行,每行3个整数 $type, l, r$ ,其中 $type$ 表示操作类型, $l, r$ 表示操作区间为 $[l, r]$ 。
Output
对于每一个 $type = 2$ 的询问,输出一个整数表示答案。
由于答案可能很大,请将答案对 $10^9 + 7$ 取模。
Sample Input |
Sample Output |
4 4 |
11 |
HINT
数据范围与约定
对于5%的数据, $N \leq 20, Q = 1$
对于10%的数据, $N \leq 1000, Q = 1$
对于20%的数据, $N \leq 10^5, Q \leq 10$
对于另外30%的数据, $1 \leq N \leq 10^5, 1 \leq Q \leq 10^5, type = 2$
对于100%的数据, $1 \leq N \leq 10^5, 1 \leq Q \leq 10^5$
题解
这道题很有意思。
首先考虑一下不带修改的解法。
设$f_{i,0}$表示$s_1...s_i$中,以$0$结尾的子序列数量;$f_{i,1}$表示$s_1...s_i$中,以$1$结尾的子序列数量。
则有方程:
若$s_i$为0:$\begin{aligned}f_{i,0}&=f_{i-1,0}+f_{i-1,1}+1\\f_{i,1}&=f_{i-1,1}\end{aligned}$
若$s_i$为1:$\begin{aligned}f_{i,0}&=f_{i-1,0}\\f_{i,1}&=f_{i-1,0}+f_{i-1,1}+1\end{aligned}$
发现这是一类线性递推,如果用一个1x3的矩阵表示原来的$f_{i,0}$与$f_{i,1}$:$\begin{pmatrix} f_0&f_1&1 \end{pmatrix}\\$(最后的1仅作为辅助计算),乘上一个3x3的转移矩阵来得到下一位的状态呢?
如果序列中这一位$s_i$为0,则在后面乘上这样一个转移矩阵$G_0$:
$$\begin{pmatrix} f_0&f_1&1\end{pmatrix}*\begin{pmatrix} 1&0&0\\1&1&0\\1&0&1 \end{pmatrix}=\begin{pmatrix} f_0+f_1+1&f_1&1\end{pmatrix}$$
如果这一位$s_i$为1,则在后面乘上另一个转移矩阵$G_1$:
$$\begin{pmatrix} f_0&f_1&1\end{pmatrix}*\begin{pmatrix} 1&1&0\\0&1&0\\0&1&1 \end{pmatrix}=\begin{pmatrix} f_0&f_0+f_1+1&1\end{pmatrix}$$
那么我们用线段树存储每一位的转移矩阵,查询时直接查询$[l,r]$的矩阵乘积,乘上初始矩阵(其实初始矩阵为$(0,0,1)$乘了相当于没乘),所以直接输出查询矩阵的$[3][1]+[3][2]$即可
处理区间数值翻转操作
最基础的想法就是,将线段树$[l,r]$叶子节点对应的转换矩阵换成另一个转换矩阵。
观察$G_0=\begin{pmatrix} 1&0&0\\1&1&0\\1&0&1 \end{pmatrix}$ 与$G_1=\begin{pmatrix} 1&1&0\\0&1&0\\0&1&1 \end{pmatrix}$
本质上,只需把第一第二行交换一下,再将第一第二列交换一下,它们都能变成对方。
第一第二行交换,相当于在$G$前乘上一个矩阵$A$。第一第二列交换,相当于在$G$后乘上这个矩阵$A$。
$$A=\begin{pmatrix} 0&1&0\\1&0&0\\0&0&1\end{pmatrix}$$
那么:
$A*G1*A=G2$ $A*G2*A=G1$
我们暂且看回原来的模型:计算一个矩阵序列。
如果要对$[l,r]$的数字翻转,假设矩阵序列是$a*b*c*d*e$,考虑如何变换:
按照我们的预想处理方式,那应该变成$(A*a*A)*(A*b*A)*(A*c*A)*(A*d*A)*(A*e*A)$。
此时我们发现,$A*A$居然是单位矩阵...
于是就变成了$A*(a*b*c*d*e)*A$。
相当于对$a*b*c*d*e$直接手动1、2行交换,1、2列交换。
回到线段树,如果要翻转,直接在对应区间维护的矩阵进行 行交换列交换,维护并下传标记即可。
神题啊!
#include <cstdio>
#include <iostream>
using namespace std;
typedef long long ll;
const int N=,Mod=1e9+;
int n,q;
char in[N];
struct Mat{
ll a[][];
void flip(){
for(int i=;i<;i++) swap(a[i][],a[i][]);
swap(a[][],a[][]);
swap(a[][],a[][]);
}
friend Mat operator * (Mat x,Mat y){
Mat ret;
for(int i=;i<;i++)
for(int j=;j<;j++){
ret.a[i][j]=;
for(int k=;k<;k++)
ret.a[i][j]=(ret.a[i][j]+x.a[i][k]*y.a[k][j])%Mod;
}
return ret;
}
};
const Mat stand[]={{,,,,,,,,},{,,,,,,,,}};
struct SegmentTree{
int root,cnt,ch[N*][],rev[N*];
Mat info[N*];
void build(int &u,int l,int r){
if(!u) u=++cnt;
if(l==r){
info[u]=stand[in[l]==''];
return;
}
int mid=(l+r)>>;
build(ch[u][],l,mid);
build(ch[u][],mid+,r);
pushup(u);
}
void flip(int u,int l,int r,int L,int R){
if(L<=l&&r<=R){
rev[u]^=;
info[u].flip();
return;
}
pushdown(u);
int mid=(l+r)>>;
if(L<=mid) flip(ch[u][],l,mid,L,R);
if(mid<R) flip(ch[u][],mid+,r,L,R);
pushup(u);
}
Mat query(int u,int l,int r,int L,int R){
if(L<=l&&r<=R) return info[u];
pushdown(u);
int mid=(l+r)>>;
if(R<=mid) return query(ch[u][],l,mid,L,R);
if(mid<L) return query(ch[u][],mid+,r,L,R);
return query(ch[u][],l,mid,L,R)*query(ch[u][],mid+,r,L,R);
}
inline void pushup(int u){
info[u]=info[ch[u][]]*info[ch[u][]];
}
inline void pushdown(int u){
if(!rev[u]) return;
rev[ch[u][]]^=; rev[ch[u][]]^=;
info[ch[u][]].flip(); info[ch[u][]].flip();
rev[u]=;
}
}seg;
int main(){
scanf("%d%d%s",&n,&q,in+);
seg.build(seg.root,,n);
int t,l,r;
while(q--){
scanf("%d%d%d",&t,&l,&r);
if(t==)
seg.flip(seg.root,,n,l,r);
else{
Mat ans=seg.query(seg.root,,n,l,r);
printf("%lld\n",(ans.a[][]+ans.a[][])%Mod);
}
}
return ;
}
奇妙代码
Subsequence Count (线段树)的更多相关文章
- HDU 6155 Subsequence Count 线段树维护矩阵
Subsequence Count Time Limit: 10000/5000 MS (Java/Others) Memory Limit: 256000/256000 K (Java/Oth ...
- [HDU6155]Subsequence Count(线段树+矩阵)
DP式很容易得到,发现是线性递推形式,于是可以矩阵加速.又由于是区间形式,所以用线段树维护. https://www.cnblogs.com/Miracevin/p/9124511.html 关键在于 ...
- HDU.6155.Subsequence Count(线段树 矩阵)
题目链接 首先考虑询问[1,n]怎么做 设 f[i][0/1]表示[1,i]以0/1结尾的不同子序列个数 则 \(if(A[i]) f[i][1] = f[i-1][0] + f[i-1][1] + ...
- FZU 2105 Digits Count(线段树)
Problem 2105 Digits Count Accept: 302 Submit: 1477 Time Limit: 10000 mSec Memory Limit : 262144 KB P ...
- 【BZOJ3638】Cf172 k-Maximum Subsequence Sum 线段树区间合并(模拟费用流)
[BZOJ3638]Cf172 k-Maximum Subsequence Sum Description 给一列数,要求支持操作: 1.修改某个数的值 2.读入l,r,k,询问在[l,r]内选不相交 ...
- CSU - 1551 Longest Increasing Subsequence Again —— 线段树/树状数组 + 前缀和&后缀和
题目链接:http://acm.csu.edu.cn/csuoj/problemset/problem?pid=1551 题意: 给出一段序列, 删除其中一段连续的子序列(或者不删), 使得剩下的序列 ...
- Codeforces 750E - New Year and Old Subsequence(线段树维护矩阵乘法,板子题)
Codeforces 题目传送门 & 洛谷题目传送门 u1s1 我做这道 *2600 的动力是 wjz 出了道这个套路的题,而我连起码的思路都没有,wtcl/kk 首先考虑怎样对某个固定的串计 ...
- CF280D k-Maximum Subsequence Sum(线段树)
在做这题时我一开始把\(tag\)写入了结构体 #include <iostream> #include <cstdio> #include <cstring> # ...
- FZU 2105Digits Count(线段树 + 成段更新)
Description Given N integers A={A[0],A[1],...,A[N-1]}. Here we have some operations: Operation 1: AN ...
- fzu 2105 Digits Count ( 线段树 ) from 第三届福建省大学生程序设计竞赛
http://acm.fzu.edu.cn/problem.php?pid=2105 Problem Description Given N integers A={A[0],A[1],...,A[N ...
随机推荐
- CSDN博客测试目录
经常写博客,但是一般没怎么注意些目录,最近看别人写的博客都有目录,所以我也想在以后写好目录,这样子也方便阅读. 这里就写一个实验: 这里一级目录 这里是一级目录下的文本.林肯公园 这里是1.1目录 这 ...
- SCOPE_IDENTITY()
@@IDENTYITY,SCOPE_IDENTITY的主要区别:在有触发器中而且触发器的内容里面含有插入标识符的操作的时候,@@IDENTITY则返回的是触发器里面新插入标识符的值而SCOPE_IDE ...
- 如何书写优雅、漂亮的SQL脚本?
本篇来聊聊如何书写漂亮.整洁.优雅的SQL脚本,下面这些是我个人总结.整理出来的.姑且做个抛砖引玉吧,呵呵,欢迎大家一起来讨论. 我们首先来看看一段创建数据表的脚本(如下所示),你是否觉得有什么不 ...
- jsp页面取值
一般就用el表达式 ${recordList[4].baseRate8.split("/")[0] } <s:date name="recordList[#id]. ...
- 重定向stdin stdout stderr |
在Linux下,当一个用户进程被创建的时候,系统会自动为该进程创建三个数据 流,也就是题目中所提到的这三个.那么什么是数据流呢(stream)? 我们知道,一个程序要运行,需要有输入.输出,如果出错, ...
- R 调用 python
上一篇说了python使用 rpy2 调用 R,这里介绍R如何调用python.R的强项在于统计方面,尤其是专业的统计分析,统计检验以及作图功能十分强大,但是在通用性方面,就远不如Python了,比如 ...
- oracle的分组查询和连接查询
分组函数: 六个常用的分组函数: AVG,SUM,MIN,MAX,COUNT,WM_CONCAT: 行转列 PS:分组函数默认会自动过滤控制,可以使用NVL函数使分组函数无法忽略空值: 未使用NVL函 ...
- 编译安装 python 2.7
下载python2.7 Python-2.7.6.tgz 下载链接:http://pan.baidu.com/s/1c0AJDDI 配置./configure 编译make 安装 make insta ...
- Laravel5.5核心架构理解
1.依赖注入 方法传入组件名,框架会自动实例化,方法内可直接使用 例如最常用的requert对象 2.服务容器 其实,Laravel 的核心就是一个 IoC 容器,Laravel 的核心本身十分轻量, ...
- CocosCreator游戏开发---菜鸟学习之路(三)如何在CocosCreator中使用Pomelo
PS(废话): 这段时间都在研究网易的Pomelo框架,作为新手小白,自然遇到了不少坑爹的事情.(当然也有可能是因为自己技术不过关的原因所以导致在很多基础的问题上纠结了很久.)网上也搜索了好久,但是基 ...