BZOJ4085: [Sdoi2015]音质检测

最新推荐文章于 2020-04-22 15:48:02 发布
最新推荐文章于 2020-04-22 15:48:02 发布
阅读量200 收藏
点赞数
文章标签： 数据结构与算法
原文链接：http://www.cnblogs.com/Yangrui-Blog/p/9623294.html
版权

BZOJ4085: [Sdoi2015]音质检测

由于这题太毒了，导致可能会被某些人卡评测，于是成了一道权限题。。。

本蒟蒻表示没钱氪金。。。

但是可以去洛谷/Vijos搞搞事。。。

但是洛谷上只能评测，题面暂时不全。。。

然而Vijos上数据范围不对，应该是洛谷上的数据范围。。。

所以没钱氪金就凑合着看吧。。。

这里附上Vijos的题面。

1958 音质检测

描述

万老板希望在新的智能音乐播放设备IPOOD中，实现对波文件音质性能的评定。

离散的波文件被考虑为长度为N的整数序列：$A_1,A_2,\cdots,A_N$​。所谓的音质性能检测，可以评定任何的一个区间范围$[L,R]$，其音质性能取决于下述评分：

$\Big\{ \sum_{L<i<R} F[A_{i-1}+1]F[A_{i+1}-1] \Big\} \mod (10^9+7)$

其中F是可归纳定义的数列，满足$F[1]=1$，$F[2]=2$且$F[k+2]=F[k+1]+aF[k]+b$对于任何$k\ge 1$成立。其中$a $和$b$为正整系数。

为了可以为用户提供更好的服务体验，并希望对给定的波文件进行修正优化。这一款设备中，还应该支持对波文件的修改。
对于给定的区间范围$[L,R]$，允许用户将$A[L]$至$A[R]$同时增加一，或同时减少一。

格式

输入格式

输入的第一行有两个正整数，波文件的总长度NN，和总的修改与询问次数QQ。

第二行有两个整数，分别表示系数aa和bb。

之后若干行，一共给出NN个正整数A_1A1​到A_NAN​，满足$1 \le A[i] \le 2*10^9$。

之后$Q$行，每行是下述三种形式之一：

$plus~L~R$：将波文件数列中下标在区间$[L,R]$内的元素每一个都加一。

$minus~L~R$：将波文件数列中下标在区间$[L,R]$内的元素每一个都减一。

$query~L~R$：询问区间$[L,R]$的音质性能评分。

修改和询问中，均保证$L\le R$，且保证$A[i]$严格大于总的修改次数加一（修改操作包括plus和minus两种）。

输出格式

输出若干行，每一行对应一次询问，输出一个整数。

样例1

样例输入1

7 7
1 0
3 4 5 6 7 8 9
query 2 4
query 3 7
plus 3 5
query 2 4
plus 4 7
query 3 7
query 1 7

样例输出1

64
1766
104
7479
7687

样例2

样例输入2

7 12
123456789 987654321
1111111111 1222222222 1333333333
1444444444 1555555555 1666666666 1777777777
query 2 4
query 3 7
plus 3 5
minus 4 6
plus 5 7
query 2 4
query 1 6
plus 4 7
minus 1 4
query 3 7
query 2 6
query 1 7

样例输出2

528103209
239947280
528103209
229970829
524160263
336413855
113033289

限制

对于15%的数据，$\(N\le 8000，Q\le 8000\)$。

对于100%的数据，$\(N\le 100000，Q\le 100000，0\le a,b\le 10^9\)$。

此外：

存在15%的数据，每一次修改的区间都是[1,N]。

还存在30%的数据，a=0，b=1。

来源

SDOI 2015 round2 day1

题解Here!

这个题简直就是道神题。。。

首先，线段树+奇奇怪怪的矩阵快速幂应该都能想到。

关键是怎么做对吧。。。

首先考虑一个性质：我们如果有数列的相邻两项$f[i],f[i+1]$。

那么用这两项向后推$k$项其线性表示系数一定。

即表示为$f[i+k]=A\times f[i]+B\times f[i+1]+C$的形式。

那么这样我们预处理这些系数就好了嘛。

注意到维护的是一个乘积的形式，那么我们要维护这个必须得维护$8$个量。

对！你没有看错！$8$个！

将其写成$3\times 3$矩阵的形式转移会比较科学。

注意$a==0$时的特判。

以上是正解。

然后某位巨佬用了$BSGS$预处理$f$函数，简直无敌了：链接。

然而我表示并不能码出来，于是自己$YY$了一个方案。

思路大致差不多，但是维护的东西多了点。

我们可以用线段树维护$f[a_{i-1}],f[a_{i-1}+1],f[a_{i-1}-1],f[a_{i+1}],f[a_{i+1}+1],f[a_{i+1}-1]$以及两两乘积和。

加一的时候可以直接用递推式用保存的$f$值求出新的值。

减一的时候也可以解方程求值。

然后开始码码码。。。

好长啊。。。

记得尽量降低时间、空间复杂度，能$return$的地方尽量$return$，以免不必要的搜索导致$TLE/MLE$。

还有一件事，注意开$long\ long$！

还有一件事，做减法的时候如果减了两次，一定要先取模，再$+MOD$，再取模！不然会求出来负数。。。

还有一件事，注意线段树的边界问题，为这个我还浪费了一次提交，药丸。。。

附上炒鸡长的代码（码风较丑不要打我。。。）：

#include<iostream>
#include<algorithm>
#include<cstdio>
#define LSON rt<<1
#define RSON rt<<1|1
#define DATA(x,i,j) a[x].data[i-1][j-1]
#define SUM(x,i,j) a[x].sum[i-1][j-1]//为了省空间只能这么干了。。。
#define ADD(x) a[x].c
#define DEL(x) a[x].d
#define LSIDE(x) a[x].l
#define RSIDE(x) a[x].r
#define WIDTH(x) (RSIDE(x)-LSIDE(x)+1)
#define MAXN 300010
#define MOD 1000000007LL
using namespace std;
int n,m;
long long A,B,inv,f[MAXN][4];
struct node{
	long long val[4][4];
	inline void clean(){
		for(int i=1;i<=3;i++)
		for(int j=1;j<=3;j++)
		val[i][j]=(i==j);
	}
	friend node operator *(const node x,const node y){
		node ret;
		for(int i=1;i<=3;i++)
		for(int j=1;j<=3;j++){
			ret.val[i][j]=0;
			for(int k=1;k<=3;k++)ret.val[i][j]=(ret.val[i][j]+x.val[i][k]*y.val[k][j]%MOD)%MOD;
		}
		ret.val[3][3]=1;
		return ret;
	}
}one,two,power[35];
struct Segment_Tree{
	long long sum[3][3],data[3][3];
	int c,d,l,r;
}a[MAXN<<2];
inline long long read(){
	long long date=0,w=1;char c=0;
	while(c<'0'||c>'9'){if(c=='-')w=-1;c=getchar();}
	while(c>='0'&&c<='9'){date=date*10+c-'0';c=getchar();}
	return date*w;
}
long long mexp(long long a,long long b,long long c){
	long long s=1;
	while(b){
		if(b&1)s=s*a%c;
		a=a*a%c;
		b>>=1;
	}
	return s;
}
node pow(long long k){
	node s;
	s.clean();
	int i=1;
	while(k){
		if(k&1)s=s*power[i];
		k>>=1;
		i++;
	}
	return s;
}
void function(int i){
	node s=one*pow(f[i][0]-2);
	f[i][1]=s.val[1][2];f[i][2]=s.val[1][1];
}
void build(){
	one.val[1][1]=2;one.val[1][2]=one.val[1][3]=1;
	two.val[1][1]=two.val[1][2]=two.val[3][3]=1;
	two.val[2][1]=A;two.val[3][1]=B;
	for(int i=1;i<=32;i++){
		power[i]=two;
		two=two*two;
	}
}
inline void pushup(int rt){
	for(int i=1;i<=3;i++)
	for(int j=1;j<=3;j++){
		DATA(rt,i,j)=(DATA(LSON,i,j)+DATA(RSON,i,j))%MOD;
		SUM(rt,i,j)=(SUM(LSON,i,j)+SUM(RSON,i,j))%MOD;
	}
}
inline void add(int rt,int l){
	for(int i=1;i<=2;i++)SUM(rt,l,i)=SUM(rt,l,i+1);
	SUM(rt,l,3)=(SUM(rt,l,2)+A*SUM(rt,l,1)%MOD+B*WIDTH(rt)%MOD)%MOD;
	if(l==1){
		for(int i=1;i<=2;i++)
		for(int j=1;j<=3;j++)
		DATA(rt,i,j)=DATA(rt,i+1,j);
		for(int i=1;i<=3;i++)DATA(rt,3,i)=(DATA(rt,2,i)+A*DATA(rt,1,i)%MOD+B*SUM(rt,2,i)%MOD)%MOD;
	}
	else{
		for(int i=1;i<=3;i++)
		for(int j=1;j<=2;j++)
		DATA(rt,i,j)=DATA(rt,i,j+1);
		for(int i=1;i<=3;i++)DATA(rt,i,3)=(DATA(rt,i,2)+A*DATA(rt,i,1)%MOD+B*SUM(rt,1,i)%MOD)%MOD;
	}
}
inline void del(int rt,int l){
	if(A==0){
		for(int i=2;i>=1;i--)SUM(rt,l,i+1)=SUM(rt,l,i);
		SUM(rt,l,1)=(SUM(rt,l,2)-B*WIDTH(rt)%MOD+MOD)%MOD;
		if(l==1){
			for(int i=2;i>=1;i--)
			for(int j=1;j<=3;j++)
			DATA(rt,i+1,j)=DATA(rt,i,j);
			for(int i=1;i<=3;i++)DATA(rt,1,i)=(DATA(rt,2,i)-B*SUM(rt,2,i)%MOD+MOD)%MOD;
		}
		else{
			for(int i=1;i<=3;i++)
			for(int j=2;j>=1;j--)
			DATA(rt,i,j+1)=DATA(rt,i,j);
			for(int i=1;i<=3;i++)DATA(rt,i,1)=(DATA(rt,i,2)-B*SUM(rt,1,i)%MOD+MOD)%MOD;
		}
		return;
	}
	for(int i=2;i>=1;i--)SUM(rt,l,i+1)=SUM(rt,l,i);
	SUM(rt,l,1)=((SUM(rt,l,3)-SUM(rt,l,2)-B*WIDTH(rt)%MOD)*inv%MOD+MOD)%MOD;
	if(l==1){
		for(int i=2;i>=1;i--)
		for(int j=1;j<=3;j++)
		DATA(rt,i+1,j)=DATA(rt,i,j);
		for(int i=1;i<=3;i++)DATA(rt,1,i)=((DATA(rt,3,i)-DATA(rt,2,i)-B*SUM(rt,2,i)%MOD)*inv%MOD+MOD)%MOD;
	}
	else{
		for(int i=1;i<=3;i++)
		for(int j=2;j>=1;j--)
		DATA(rt,i,j+1)=DATA(rt,i,j);
		for(int i=1;i<=3;i++)DATA(rt,i,1)=((DATA(rt,i,3)-DATA(rt,i,2)-B*SUM(rt,1,i)%MOD)*inv%MOD+MOD)%MOD;
	}
}
inline void pushdown_sign(int rt,int l,int c){
	if(!c)return;
	if(c>0)for(int i=1;i<=c;i++)add(rt,l);
	else for(int i=-1;i>=c;i--)del(rt,l);
}
inline void pushdown(int rt){
	if(LSIDE(rt)==RSIDE(rt))return;
	ADD(LSON)+=ADD(rt);DEL(LSON)+=DEL(rt);
	pushdown_sign(LSON,1,ADD(rt));pushdown_sign(LSON,2,DEL(rt));
	ADD(RSON)+=ADD(rt);DEL(RSON)+=DEL(rt);
	pushdown_sign(RSON,1,ADD(rt));pushdown_sign(RSON,2,DEL(rt));
	ADD(rt)=DEL(rt)=0;
}
void buildtree(int l,int r,int rt){
	LSIDE(rt)=l;RSIDE(rt)=r;ADD(rt)=DEL(rt)=0;
	if(l==r){
		for(int i=1;i<=3;i++){
			SUM(rt,1,i)=f[l-1][i];
			SUM(rt,2,i)=f[l+1][i];
		}
		for(int i=1;i<=3;i++)
		for(int j=1;j<=3;j++)
		DATA(rt,i,j)=SUM(rt,1,i)*SUM(rt,2,j)%MOD;
		return;
	}
	int mid=LSIDE(rt)+RSIDE(rt)>>1;
	buildtree(l,mid,LSON);
	buildtree(mid+1,r,RSON);
	pushup(rt);
}
void update_add(int l,int r,int c,int rt){
	if(l<=LSIDE(rt)&&RSIDE(rt)<=r){
		ADD(rt)+=c;
		pushdown_sign(rt,1,c);
		return;
	}
	pushdown(rt);
	int mid=LSIDE(rt)+RSIDE(rt)>>1;
	if(l<=mid)update_add(l,r,c,LSON);
	if(mid<r)update_add(l,r,c,RSON);
	pushup(rt);
}
void update_del(int l,int r,int c,int rt){
	if(l<=LSIDE(rt)&&RSIDE(rt)<=r){
		DEL(rt)+=c;
		pushdown_sign(rt,2,c);
		return;
	}
	pushdown(rt);
	int mid=LSIDE(rt)+RSIDE(rt)>>1;
	if(l<=mid)update_del(l,r,c,LSON);
	if(mid<r)update_del(l,r,c,RSON);
	pushup(rt);
}
long long query(int l,int r,int rt){
	if(l>r)return 0;
	long long ans=0;
	if(l<=LSIDE(rt)&&RSIDE(rt)<=r)return DATA(rt,3,1);
	pushdown(rt);
	int mid=LSIDE(rt)+RSIDE(rt)>>1;
	if(l<=mid)ans=(ans+query(l,r,LSON))%MOD;
	if(mid<r)ans=(ans+query(l,r,RSON))%MOD;
	return ans;
}
void work(){
	char ch[7];
	int x,y,l,r;
	while(m--){
		scanf("%s",ch);x=read();y=read();
		if(ch[0]=='p'){
			l=x+1;r=y+1;
			if(r>n)r=n;
			if(l<=r)update_add(l,r,1,1);
			l=x-1;r=y-1;
			if(l<1)l=1;
			if(l<=r)update_del(l,r,1,1);
		}
		else if(ch[0]=='m'){
			l=x+1;r=y+1;
			if(r>n)r=n;
			if(l<=r)update_add(l,r,-1,1);
			l=x-1;r=y-1;
			if(l<1)l=1;
			if(l<=r)update_del(l,r,-1,1);
		}
		else printf("%lld\n",query(x+1,y-1,1));
	}
}
void init(){
	n=read();m=read();A=read();B=read();
	inv=mexp(A,MOD-2,MOD);
	build();
	for(int i=1;i<=n;i++){
		f[i][0]=read();
		function(i);
		f[i][3]=(f[i][2]+f[i][1]*A%MOD+B)%MOD;
	}
	buildtree(1,n,1);
}
int main(){
	init();
	work();
    return 0;
}

 

转载于:https://www.cnblogs.com/Yangrui-Blog/p/9623294.html

weixin_30342827
关注
bzoj 4085: [Sdoi2015]音质检测
reflash的博客
04-17 1650
空？
【SDOI2015】【BZOJ4085】音质检测quality
CreationAugust is 14 years old forever
04-18 1431
4085: [Sdoi2015]qualityTime Limit: 80 Sec Memory Limit: 512 MB Submit: 79 Solved: 27 [Submit][Status][Discuss] DescriptionInputOutputSample InputSample OutputHINTSource题面请去Vijos看.(目前Vijos需要翻墙线段树+矩
洛谷 P3328 【[SDOI2015]音质检测】
12-15 163
这题我做的好麻烦啊。。。 一开始想分块来着，后来发现可以直接线段树 首先考虑一个性质，我们如果有数列的相邻两项f[i]和 f[i+1]那么用这两项向后推k项其线性表示系数一定（表示为f[i+k]=a∗f[i]+b∗f[i+1]+c的形式)，那么这样我们预处理这些系数，注意到维护的是一个乘积的形式，那么我们要维护这个必须得维护8个量，将其写成3 * 3矩阵的形式转移会比较科学,注意a=0的...
音质检测的算法创新与实现
audio_mqiu的专栏
02-26 3194
如需转载请注明出处！ 近来做音质检测相关的研究工作，预期以算法自动筛查的方法过滤掉大规模音频库中的参数劣质音频，并作为后续新上传音频的音质检测门槛，最终在一定程度上保障参数优质音频。 独创了音质检测算法（至今已有的资料还未看到有相同或相似算法），通过本地小规模音频库检测，性能优异（召回率=100%， 正确率=100%） 目前该音质检测系统覆盖范围： 1. 量化精度低音频，
BZOJ 4085:[Sdoi2015]quality(round 2 音质检测)(数据结构)
weixin_30773135的博客
06-10 110
居然在考场上把这道题打出来了觉得自己也是有点吊啊（虽然后面就没时间做其他题了囧而且还被卡常数了。。。） 题解自己写了一份TEX的就直接放上来吧。。。。 好啦，在谈点什么别的 什么？你在bz上TLE了？注意一下你的矩阵乘法，这个程序的大部分时间几乎都是跑矩阵乘法的，我是从900000次到600000次在到300000次最后预处理那些2的次幂才过的。把OJ卡了好久真是对不起啊QAQ COD...
BZOJ 4698: Sdoi2008 Sandy的卡片
LZJ209的博客
12-30 1万+
这道题看到加几个数就能变成另一个数就知道是差分啦，然后显而易见的就是最后的答案为所有串的最长公共串+1，那我们之后怎么办呢，我们先对所有串建一个广义后缀自动机，然后在自动机上搜寻每一个单词，对单词的每一个节点一直跳parent指针，路径上每一个点的次数加一，统计出所有访问次数为字符串总数的长度最大值就好啦，为了不重复访问，我们将每一个点弄一个标记，为了防止繁琐的清空操作，弄一个时间戳就好啦。我这个傻
bzoj 4818: [Sdoi2017]序列计数 动态规划+矩阵乘法
beginend
04-11 1768
题意Alice想要得到一个长度为nn的序列，序列中的数都是不超过mm的正整数，而且这nn个数的和是pp的倍数。 Alice还希望，这nn个数中，至少有一个数是质数。 Alice想知道，有多少个序列满足她的要求。 对100%100\%的数据，1≤n≤109,1≤m≤2×107,1≤p≤1001\leq n \leq 10^9,1\leq m \leq 2×10^7,1\leq p\leq 100
bzoj3943: [Usaco2015 Feb]SuperBull
Cynthia_wjyi
10-25 984
3943: [Usaco2015 Feb]SuperBull Time Limit: 10 Sec  Memory Limit: 128 MB Submit: 126  Solved: 86 [Submit][Status][Discuss] Description Bessie and her friends are playing hoofball in the annual S
音频质量-一种异音分析算法原理.pdf
10-27
本文描述了一种高级用于音频质量评估的杂音、异音检测原理及实现方法
SDOI2015 总结
qq_42835841的博客
04-22 237
文章目录4343【SDOI2015】排序4344【SDOI2015】寻宝游戏4346【SDOI2015】星际战争4347【SDOI2015】约数个数和4348【SDOI2015】道路修建5432【SDOI2015】音质检测5434【SDOI2015】双旋转字符串 4343【SDOI2015】排序 dfsdfsdfs，发现顺序无影响，所以按长度依次dfsdfsdfs1,2,4⋯ ,2n1,2,4\c...
[BZOJ4085][Sdoi2015]quality && 线段树 (TLE)
shiyukun1998的专栏
06-11 970
线段树维护F[ai-1], F[ai-1 + 1], F[ai-1 - 1], F[ai+1], F[ai+1 + 1], F[ai+1 - 1] 以及两两乘积和 加一的时候可以直接用递推式用保存的F值求出新的值 减一的时候也可以解方程求值 然后 写码农题 然而此题常数太卡 在jzoj上都600ms了 BZ上还是过不去 弃疗了 #include #include #include #i
3992: [SDOI2015]序列统计
蒟蒻的Bzoj记录
02-06 252
题目链接题目大意：给定n(n<=10^9)，质数m(3<=m<=8000)，1<=x=m，以及一个[0,m-1]区间内的集合S，求有多少长度为n的数列满足每个元素都属于集合S且所有元素的乘积mod m后=x题解：f[i][j]表示放到前i个数乘积为j的方案数，O(nm2)O(nm^2)转移显然 求出MM的原根gg以及每个数ii的指标ind[i]ind[i] (a×b)modM=x(a\times
【SDOI2015】BZOJ3992 序列统计题解（DP+原根+NTT）
hzk_cpp的博客
03-13 211
题目：BZOJ3992. 题目大意：给定集合SSS，求数列aia_iai​的数量，其中∏i=1nai≡x&amp;amp;ThinSpace;&amp;amp;ThinSpace;(mod&amp;amp;ThinSpace;&amp;amp;ThinSpace;m)\prod_{i=1}^{n}a_i\equiv x\,\,(mod\,\,m)∏i=1n​ai​≡x(modm). 1≤m≤8000,1≤n≤109,1≤x≤m1...
Reverb音效算法研究实现
audio_mqiu的专栏
02-25 9040
如需转载请注明出处！ 1. Reverb物理原型 在一定的空间范围内（如房间中）进行语音/音频播放时(讲话或者Audio播放)。对于语音/音频接收端而言，其接收到的音频信号除了播放端直接传递过来的信号（直达声）外，还包含了经过各种复杂物理反射后的信号（反射声）。反射声信号要延迟于直达声信号到达，并且由于物理反射作用其能量要发生衰减。并且空间模型不同（房间、礼堂），反射声的延迟及能量衰减
【SDOI2015】游记Day0
CreationAugust is 14 years old forever
04-10 1149
下午13:00从老校出发踏上征程，路上除了睡觉就是看风景听音乐和队友互相嘲讽 然后16:30左右抵达目的地 入住了一家宾馆叫海通酒店还是什么 一开始一看在街里面而且外观看起来有点破破烂烂。。。印象不太好 但是入住之后就感觉超级high 有网线有电源有干净的卫浴和床铺 灯光是那种淡淡的昏黄色，感觉不错蛮有气氛的，这种光休息起来也比较舒服吧 说实话做了一下午车感觉挺累的但是到了晚上吃饭的时
音频节奏检测（Onset Detection）
热门推荐
矩阵实验室
12-17 3万+
1. 前言最近市场上出现一些多个视频拼接而成MV，其原理是根据音频的节拍变换切换视频。 我在这里讲述下如何进行音频节拍检测。2. 音频检测一般流程3.3.1 原始音频频谱以1024为窗口（即每次读取1024个采样点）进行量化 WaveDecoder decoder = new WaveDecoder( new FileInputStream( "samples/sample.wav"
[SDOI2015]序列统计 （NTT：快速数论变换 + 快速幂）
总理同学的编程尝试
12-01 601
调试了一个多小时可算是调试出来了，我的NTT真是漏洞百出，还好最后AC了…
利用梅尔倒谱系数（MFCC）及空间聚类算法实现音色识别
zclhit
02-08 9217
写在前面 2016年4月参加了哈尔滨工业大学深圳研究生院举办的创新创业比赛，司职算法组长，切入点定在了音色识别和相似明星音才艺展示推荐算法上，不才，拿到了一等奖，趁佳节未散与大家分享。 项目进度安排 2016年1月~2016年3月：前期工作中了解学习了语音信号处理的基本原理，查阅有关文献了解到了声音音色信息的描述方式，梅尔（Mel）倒谱系数的意义，推导了相关公式并实现了MFCC（Mel Fr
BZOJ2839:集合计数(容斥,组合数学)
最新发布
06-08
这个问题是一个比较经典的容斥原理的应用题目。我们可以考虑使用容斥原理来求解这个问题。具体来说，我们可以考虑求出所有集合的个数，然后减去至少包含一个特定元素的集合的个数，再加上至少包含两个特定元素的集合...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值