[NOI2017] 整数，蚯蚓排队，泳池

已于 2022-08-19 20:09:57 修改
阅读量522 收藏
点赞数 2
分类专栏： 信息竞赛解题 文章标签： 算法 数据结构 哈希算法 动态规划 多项式
于 2022-08-17 17:12:47 首次发布
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/weixin_43960287/article/details/126388858
版权

前言

本来只是想写一下[NOI2017] 泳池的题解，不过前两题题解太短了干脆挼一起了。

[NOI2017] 整数

考虑大力数据结构，直接开一个长度为 30 n + 233 30n+233 30n+233（或者+114514） 的线段树，用最朴素的方法维护每个bit的值。

加上一个正整数 a a a 可以直接把序列上对应位置长为 log ⁡ a \log a loga 的那一段数拿下来做加法，然后把前 log ⁡ a \log a loga 位再贴回去。最后考虑和的第 log ⁡ a + 1 \log a+1 loga+1 位，最多进位为1，我们找到序列上前面第一个0的位置然后做一个区间取反即可。

减法类似，若需要补位则找前面第一个1的位置即可。

复杂度 O ( n ( log ⁡ n + log ⁡ a ) ) O(n(\log n+\log a)) O(n(logn+loga))，理论上常数较大，但是用zkw线段树可以轻松进900ms。

#include<bits/stdc++.h>//JZM yyds!!
#define ll long long
#define lll __int128
#define uns unsigned
#define fi first
#define se second
#define IF (it->fi)
#define IS (it->se)
#define END putchar('\n')
#define lowbit(x) ((x)&-(x))
#define inline jzmyyds
using namespace std;
const int MAXN=1e6+5;
const ll INF=1e18;
ll read(){
	ll x=0;bool f=1;char s=getchar();
	while((s<'0'||s>'9')&&s>0){if(s=='-')f^=1;s=getchar();}
	while(s>='0'&&s<='9')x=(x<<1)+(x<<3)+(s^48),s=getchar();
	return f?x:-x;
}
int ptf[50],lpt;
void print(ll x,char c='\n'){
	if(x<0)putchar('-'),x=-x;
	ptf[lpt=1]=x%10;
	while(x>9)x/=10,ptf[++lpt]=x%10;
	while(lpt>0)putchar(ptf[lpt--]^48);
	if(c>0)putchar(c);
}
const int MAXM=1<<25|(MAXN<<5);
int n,m;
struct zkw{
	bool f[MAXM],h[MAXM],lz[MAXM];
	int p,dep;
	void init(int n){for(p=1,dep=0;p<n+2;p<<=1,dep++);}
	void cover(int x,bool e){f[x]=h[x]=e,lz[x]=1;}
	void upd(int x){f[x]=f[x<<1]&f[x<<1|1],h[x]=h[x<<1]|h[x<<1|1];}
	void pushd(int x){
		if(lz[x]){
			f[x<<1]=f[x<<1|1]=h[x<<1]=h[x<<1|1]=f[x];
			lz[x<<1]=lz[x<<1|1]=1,lz[x]=0;
		}
	}
	void add(int l,int r,bool e){
		l=p^(l-1),r=p^(r+1);
		for(int i=dep;i>0;i--)pushd(l>>i),pushd(r>>i);
		while(l^1^r){
			if(~l&1)cover(l^1,e);
			if(r&1)cover(r^1,e);
			l>>=1,r>>=1,upd(l),upd(r);
		}for(l>>=1;l;l>>=1)upd(l);
	}
	uns ll query(int l,int r){
		l^=p,r^=p;
		for(int i=dep;i>0;i--)
			for(int j=l>>i,lm=r>>i;j<=lm;j++)pushd(j);
		uns ll res=0;
		for(int i=r;i>=l;i--)res=res<<1|f[i];
		return res;
	}
	void fz(int l,int r,uns ll d){
		l^=p,r^=p;
		for(int i=dep;i>0;i--)
			for(int j=l>>i,lm=r>>i;j<=lm;j++)pushd(j);
		for(int i=l;i<=r;i++,d>>=1)f[i]=h[i]=d&1;
		for(l>>=1,r>>=1;l;l>>=1,r>>=1)
			for(int j=l;j<=r;j++)upd(j);
	}
	int sch0(int l){
		for(int i=(l^=p,dep);i>0;i--)pushd(l>>i);
		for(;l>1;l>>=1)if((~l&1)&&!f[l^1]){l^=1;break;}
		for(;l<p;pushd(l),l<<=1,l^=f[l]);
		return l^p;
	}
	int sch1(int l){
		for(int i=(l^=p,dep);i>0;i--)pushd(l>>i);
		for(;l>1;l>>=1)if((~l&1)&&h[l^1]){l^=1;break;}
		for(;l<p;pushd(l),l<<=1,l^=!h[l]);
		return l^p;
	}
	void chg(int x,bool e){
		x^=p;
		for(int i=dep;i>0;i--)pushd(x>>i);
		f[x]=h[x]=e;
		for(int siz=(x>>=1,1);x;x>>=1,siz<<=1)upd(x);
	}
	bool sch(int x){
		x^=p;
		for(int i=dep;i>0;i--)pushd(x>>i);
		return f[x];
	}
}T;
int main()
{
	m=read(),read(),read(),read();
	n=m*30+233,T.init(n);
	for(int id=1;id<=m;id++){
		int op=read();
		if(op==1){
			ll a=read();
			int l=read()+1,r=l+29;
			if(a==0)continue;
			ll b=T.query(l,r),mx=1ll<<30;
			T.fz(l,r,(a+b+mx)&(mx-1));
			if(a+b>=mx){
				int c=T.sch0(r);
				if(c>r+1)T.add(r+1,c-1,0);
				T.chg(c,1);
			}else if(a+b<0){
				int c=T.sch1(r);
				if(c>r+1)T.add(r+1,c-1,1);
				T.chg(c,0);
			}
		}else print(T.sch(read()+1));
	}
	return 0;
}

[NOI2017] 蚯蚓排队

这题没啥好说的，直接维护链表+Hash即可。

代码我没写，看粪兔的。

[NOI2017] 泳池

看到这个“面积最大值恰好为 k k k”，可以很自然地转化为“面积都不超过 k k k”减去“面积都不超过 k − 1 k-1 k−1”。

考虑怎么求矩形面积不超过 k k k 的概率。我们知道一种求最大矩形面积的方法是把每个位置从岸边出发的最长安全距离弄成整数序列，然后用笛卡尔树，所以可以类似地定义一个区间DP来求。由于每个位置等价，所以直接设状态为区间长度：

令 d p [ i ] [ j ] dp[i][j] dp[i][j] 表示长度为 i i i 的区间，区间最小值为 j j j 的概率，我们只需要满足 i × j ≤ k i\times j\le k i×j≤k，转移就直接枚举第一个取到最小值的位置，用一下后缀和优化即可。这里我们只需要考虑最小值 > 0 >0 >0 的区间，所以区间长度不超过 k k k。

注意到状态总数是 ∑ i = 1 k k i = k log ⁡ k \sum_{i=1}^k\frac{k}{i}=k\log k ∑i=1k​ik​=klogk，且转移总复杂度为 O ( ∑ i = 1 k k i ⋅ i ) = O ( k 2 ) O(\sum_{i=1}^k\frac{k}{i}\cdot i)=O(k^2) O(∑i=1k​ik​⋅i)=O(k2) 可以通过。

显然整个序列被值为0的位置划分成了若干段，我们需要对每种划分组合的概率求和。

遇到划分计数的问题又不得不用到生成函数了。设 f ( x ) = ∑ i = 1 k + 1 ( ∑ j > 0 d p [ i − 1 ] [ j ] ) ⋅ ( 1 − q ) x i f(x)=\sum_{i=1}^{k+1}(\sum_{j>0}dp[i-1][j])\cdot (1-q)x^i f(x)=∑i=1k+1​(∑j>0​dp[i−1][j])⋅(1−q)xi 表示每种长度的非0区间后跟一个0的概率的生成函数，那么
A n s = 1 1 − q [ x n + 1 ] ( 1 + f ( x ) + f 2 ( x ) + . . . ) = 1 1 − q [ x n + 1 ] 1 1 − f ( x ) Ans=\frac{1}{1-q}[x^{n+1}]\left(1+f(x)+f^2(x)+...\right)=\frac{1}{1-q}[x^{n+1}]\frac{1}{1-f(x)} Ans=1−q1​[xn+1](1+f(x)+f2(x)+...)=1−q1​[xn+1]1−f(x)1​还剩下一个求多项式分式 1 1 − f ( x ) \frac{1}{1-f(x)} 1−f(x)1​ 远项的问题，用波斯坦-茉莉算法即可。

多项式长度 O ( k ) O(k) O(k)，所以可以暴力卷积。

#include<bits/stdc++.h>//JZM yyds!!
#define ll long long
#define lll __int128
#define uns unsigned
#define fi first
#define se second
#define IF (it->fi)
#define IS (it->se)
#define END putchar('\n')
#define lowbit(x) ((x)&-(x))
#define inline jzmyyds
using namespace std;
const int MAXN=1e3+10;
const ll INF=1e18;
ll read(){
	ll x=0;bool f=1;char s=getchar();
	while((s<'0'||s>'9')&&s>0){if(s=='-')f^=1;s=getchar();}
	while(s>='0'&&s<='9')x=(x<<1)+(x<<3)+(s^48),s=getchar();
	return f?x:-x;
}
int ptf[50],lpt;
void print(ll x,char c='\n'){
	if(x<0)putchar('-'),x=-x;
	ptf[lpt=1]=x%10;
	while(x>9)x/=10,ptf[++lpt]=x%10;
	while(lpt>0)putchar(ptf[lpt--]^48);
	if(c>0)putchar(c);
}
const ll MOD=998244353;
ll ksm(ll a,ll b,const ll&mo=MOD){
	ll res=1;
	for(;b;b>>=1,a=a*a%mo)if(b&1)res=res*a%mo;
	return res;
}
ll n,p[MAXN],dp[MAXN][MAXN],P[MAXN],Q[MAXN],f[MAXN<<1];
int k;
ll solve(int m){
	if(m<0)return 0;
	if(!m)return ksm(p[0],n);
	memset(dp,0,sizeof(dp));
	for(int i=0;i<=m+1;i++)dp[0][i]=1;
	for(int i=1;i<=m;i++)
		for(int j=m/i;j>=0;j--){
			ll&f=dp[i][j];
			for(int k=1;k<=i;k++)
				(f+=dp[k-1][j+1]*p[j]%MOD*dp[i-k][j])%=MOD;
			(f+=dp[i][j+1])%=MOD;
		}
	memset(P,0,sizeof(P));
	memset(Q,0,sizeof(Q));
	P[0]=Q[0]=1;
	for(int i=1;i<=m+1;i++)Q[i]=dp[i-1][1]*(MOD-p[0])%MOD;
	for(int N=n+1;N;N>>=1){
		for(int i=0;i<=(m+2)<<1;i++)f[i]=0;
		for(int i=0;i<=m+1;i++)
			for(int j=0;j<=m+1;j++)
				(f[i+j]+=P[i]*((j&1)?MOD-Q[j]:Q[j]))%=MOD;
		for(int i=0;i<=m+1;i++)P[i]=f[i<<1|(N&1)];
		for(int i=0;i<=(m+2)<<1;i++)f[i]=0;
		for(int i=0;i<=m+1;i++)
			for(int j=0;j<=m+1;j++)
				(f[i+j]+=Q[i]*((j&1)?MOD-Q[j]:Q[j]))%=MOD;
		for(int i=0;i<=m+1;i++)Q[i]=f[i<<1];
	}return P[0]*ksm(Q[0]*p[0]%MOD,MOD-2)%MOD;
}
int main()
{
	n=read(),k=read(),p[0]=read();
	p[0]=(MOD+1-ksm(read(),MOD-2)*p[0]%MOD)%MOD;
	for(int i=1;i<=k;i++)p[i]=ksm(MOD+1-p[0],i)*p[0]%MOD;
	print((solve(k)-solve(k-1)+MOD)%MOD);
	return 0;
}
偶耶XJX
关注
  • 2
    点赞
  • 踩
  • 0
    收藏
    觉得还不错? 一键收藏
  • 1
    评论
专栏目录
noi 2017 简要题解
执子之手,与子偕老
01-10 908
回顾noi 2017 DAY1 整数 压位维护序列。用线段树维护一段0后第一个1，一段A-1后第一个&amp;lt;(A - 1)的数 推推转移式子即可 蚯蚓排队 维护hash值，map到1023333的hash表里。 每次分裂只会有k2个串被删掉，合并又至多多出k2个。 O(nk+k^2 * c) 泳池 直接用大佬的吧 感觉dp还有一种推法，从左往右dp，但是方程没有这么清晰 常系数递推可以用特征多项式...
NOI2017题解合集（施工中）
待成熟的葡萄
07-28 4701
Day1queueAnalysis先考虑一种最暴力的算法：使用哈希表储存所有长度不超过K=50K=50的子串，合并和分裂时我们只用修改交界处的K2K^2个子串。时间复杂度O(mK2)O(mK^2)。 考虑精细地实现程序：合并时我们只插入原本没有插入过的，也就是严格跨边界的。虽然单次复杂度可能达到O(K2)O(K^2)，但是均摊意义下是可以保证复杂度的： 设势函数Φ\Phi表示哈希表里面的串的个数
[NOI2017]蚯蚓排队 题解
Hacheylight
01-05 312
题目就是要求我们能够支持三个操作： 1.把两个字符串合并 2.把两个字符串分裂 3.询问有多少一些个串是这么多串中出现了多少次 我们发现数据范围中的K很小，于是我们就考虑k与什么有关，也就是 我们发现第3中操作可以通过哈希表解决，（显而易见的吧） 那怎么维护哈希表？ 我们思考，两个字符串a,ba,ba,b合并会增加哪些哈希值？ 其实就是那些即有aaa的部分又有bbb部分的子串。 这样的串 ...
[NOI 2017]蚯蚓排队
weixin_30421525的博客
07-09 115
Description 题库链接 蚯蚓幼儿园有 \(n\) 只蚯蚓。幼儿园园长神刀手为了管理方便，时常让这些蚯蚓们列队表演。 所有蚯蚓用从 \(1\) 到 \(n\) 的连续正整数编号。每只蚯蚓的长度可以用一个正整数表示，根据入园要求，所有蚯蚓的长度都不超过 \(6\) 。神刀手希望这些蚯蚓排成若干个队伍，初始时，每只蚯蚓各自排成一个仅有一只蚯蚓的队伍，该蚯蚓既在队首，也在队尾。 神刀手将会依次进...
Luogu3823 NOI2017 蚯蚓排队
wennitao的博客
07-07 218
题面：luogu3823 题意：nnn只蚯蚓，每只蚯蚓有一个≤6\le 6≤6的长度，初始时每只蚯蚓一支队伍。给出mmm个操作。 三种操作： 将以iii结尾的队伍和以jjj开头的队伍合并，且iii的队伍在前。 将iii和iii后面一只蚯蚓处断开，分为两支队伍。 定义以第xxx只蚯蚓开始的长度为kkk的数字串为：从xxx开始的连续kkk只蚯蚓，它们的长度数字连成的字符串。定义给定一个字符串sss和...
BZOJ4944: [Noi2017]泳池
weixin_30343157的博客
07-30 112
BZOJ4944: [Noi2017]泳池 题目背景 久莲是个爱玩的女孩子。 暑假终于到了，久莲决定请她的朋友们来游泳，她打算先在她家的私人海滩外圈一块长方形的海域作为游泳场。 然而大海里有着各种各样的危险，有些地方水太深，有些地方有带毒的水母出没。 她想让圈出来的这一块海域都是安全的。 题目描述 经过初步的分析，她把这块海域抽象成了一个底边长为NN米，高为10011001米...
NOI2017天津市选拔赛(TJOI2017)题目及测试数据
07-04
【标题】"NOI2017天津市选拔赛(TJOI2017)题目及测试数据"涉及的是中国信息学奥林匹克竞赛（National Olympiad in Informatics, NOI）的一个地方选拔赛，即2017年天津赛区的比赛。NOI是由中国计算机学会主办的一项...
noi模拟题2
07-22
全国信息学奥林匹克（NOI）是一项面向中学生的编程竞赛，旨在培养青少年的计算机科学技能和创新能力。"noi模拟题2"很可能是一个包含了与NOI竞赛相关的问题集，供参赛者或热衷于编程竞赛的人进行练习和提高。这个...
NOI 国家集训队论文 1999-2009 2013-2017
04-18
这些资料主要涵盖了1999年至2017年间NOI（全国青少年信息学奥林匹克竞赛）国家集训队的论文，以及2014年和2015年的NOI冬令营讲者PPT，是信息学竞赛领域的重要参考资料。NOI是中国最高级别的中学生信息学竞赛，而国家...
IBM智能运维分析NOI-2017(1).pptx
最新发布
01-25
IBM智能运维分析NOI-2017(1)IBM智能运维分析NOI-2017(1)IBM智能运维分析NOI-2017(1)IBM智能运维分析NOI-2017(1)IBM智能运维分析NOI-2017(1)IBM智能运维分析NOI-2017(1)IBM智能运维分析NOI-2017(1)IBM智能运维分析NOI...
noi2017初赛c语言试题,NOIP2017普及组初赛试题及答案
weixin_35713476的博客
05-22 723
原标题：NOIP2017普及组初赛试题及答案清北学堂信息学金牌教研团今天为学生们整理了NOIP2017普及组初赛试题及答案，供学生们参考哦！——此文2100多文字，大概需要60分钟看完一、单项选择题(共 20 题，每题 1.5 分，共计 30 分；每题有且仅有一个正确选项)1.在 8 位二进制补码中，10101011 表示的数是十进制下的( )。A. 43 B. -85 C. -43 D. -84...
NOI2017 游记
lych的博客
07-21 5115
很幸运地搭上了前往绍一的末班车，抱着稳赚不亏的心态踏入了绍一        夕阳余晖洒下，高高的玻璃外墙上，泛起点点金光 Day -1        报到日        到了绍一已经错过了饭店，在外吔饭又辗转不少，回到寝室有点头痛，，qaq。于是睡了一觉，不知道和我们同寝室的另外两位绍一大佬怎么样啊，，        然而事实是他们至始至终都没有出现过，，        晚上开了一道
noi2017解题报告（部分）
wzf_2000的博客
08-07 2418
D1T1: 每30位压成一个int，然后压位线段树就可以过了。#include <bits/stdc++.h> #define gc getchar() #define N 4500009 #define mid (l+r>>1) #define ll long long using namespace std; const int Max=1<<30; int n,m,t1,t2,t3; str
【NOI2017】 整数（二进制计数器）
DD(XYX)的硫氰化氘博客
05-28 204
摘要（选填）：rt
[NOI2017]整数
weixin_33697898的博客
05-28 128
https://www.luogu.org/problemnew/show/P3822 题解 如果只有加法没有减法，那么我们朴素进位的时间复杂度是均摊\(O(1)\)的。 那么减法的话我们维护一个增量的数，一个减量的数，就是一个存所有\(a>0\)的变量，一个存所有\(a<0\)的变量。 那么我们考虑一个位置的数是\(0\)还是\(1\)，我们先对这个位置做个减法，然后唯一对这个值有影...
[字符串HASH][复杂度分析] NOI2017 .Day1 T2 蚯蚓排队
CE玩家
07-24 1542
因为k只有50，所以先把所有询问的字符串拆开，加到hash表里，然后用链表模拟整个过程，暴力把新产生的或消失的字符串在hash表里更新答案，复杂度就是对的 具体做法和时间复杂度分析可以看lzz的知乎回答 https://www.zhihu.com/question/62597216#include <cstdio> #include <iostream> #include <algorithm>
【BZOJ4943】【NOI2017】蚯蚓排队（哈希）
小蒟蒻yyb的博客
07-03 380
题面 BZOJ 洛谷 UOJ 题解 记得去年看网络同步赛的时候是一脸懵逼的。 昨天看到zsyzsyzsy做了，今天就看了看。。 这不是HashHashHash的傻逼题吗。。。 呜。。。 开链表存一下前驱后继，因为询问长度不超过505050 那么，总的串的个数不超过50n50n50n，暴力hashhashhash一下然后拿哈希表存一下就好了。 我开始用的map然后TLE飞了...
[NOI2017]蚯蚓排队
weixin_34269583的博客
05-17 96
题意：有一些蚯蚓，每个蚯蚓都有一个小于等于6的值，有m次操作 　　操作1：将两队蚯蚓连到一起 　　操作2：在蚯蚓x的队列中将x及x之前和x之后分为两队 　　操作3：给一个字符串s和整数k，对于s每个长度为k的字串，看有多少队蚯蚓中有这个子串，将每个子串的值乘起来作为答案 　　观察数据，有一个条件是k<=50，所以直接采用哈希做，每次合并就暴力枚举新产生的长度<=50的串，...
[noi2022] 二次整数规划问题
09-18
二次整数规划问题是指在给定的约束条件下，要求在整数解中最小化二次目标函数的优化问题。它是数学规划中一类重要的问题。 二次整数规划可以形式化表示为： minimize f(x) = ½ x^T Q x + c^T x subject to Ax ≤ ...

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论 1
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值