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

已于 2022-08-19 20:09:57 修改
阅读量563 收藏
点赞数 2
分类专栏： 信息竞赛解题 文章标签： 算法 数据结构 哈希算法 动态规划 多项式
于 2022-08-17 17:12:47 首次发布
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/weixin_43960287/article/details/126388858
版权
这篇博客详细介绍了如何利用线段树解决ZOJ NOI2017中关于整数操作的问题，包括加减法的高效实现。同时，博主还提及了蚯蚓排队问题的解决方案，通过维护链表和哈希表来求解。此外，文章还探讨了游泳池面积概率问题，转换为求不超过k的矩形面积概率，并利用生成函数和多项式运算求解。
摘要由CSDN通过智能技术生成

前言

本来只是想写一下[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
关注
专栏目录
【BZOJ4944】【NOI2017】泳池 概率DP 常系数线性递推 特征多项式 多项式取模
ez_yww的博客
11-30 2455
题目大意 　　有一个1001×n1001×n1001\times n的的网格，每个格子有qqq的概率是安全的，1−q1−q1-q的概率是危险的。 　　定义一个矩形是合法的当且仅当： 这个矩形中每个格子都是安全的 必须紧贴网格的下边界 　　问你最大的合法子矩形大小为kkk的概率是多少。 　　n≤109,k≤1000n≤109,k≤1000n\leq {10}^9,k...
【BZOJ4943】【UOJ315】【NOI2017】蚯蚓
cz_xuyixuan的博客
05-30 396
【题目链接】BZOJUOJ【思路要点】先说做法：用链表模拟蚯蚓的队列。对于每一个1号操作，计算所有新增的，长度在50以内的子串的哈希值，并加入哈希表中。对于每一个2号操作，计算所有被切断的，长度在50以内的子串的哈希值，并在哈希表中除去。对于每一个3号操作，计算每个子串的哈希值，并在哈希表中找到它们的出现次数，并相乘得到答案。显然，2号操作的总复杂度为\(O(CK^2)\)，3号操作的总复杂度为\...
[NOI2017]蚯蚓排队
weixin_34269583的博客
05-17 116
题意：有一些蚯蚓，每个蚯蚓都有一个小于等于6的值，有m次操作 　　操作1：将两队蚯蚓连到一起 　　操作2：在蚯蚓x的队列中将x及x之前和x之后分为两队 　　操作3：给一个字符串s和整数k，对于s每个长度为k的字串，看有多少队蚯蚓中有这个子串，将每个子串的值乘起来作为答案 　　观察数据，有一个条件是k<=50，所以直接采用哈希做，每次合并就暴力枚举新产生的长度<=50的串，...
NOI2017 蚯蚓排队
Ac_chen_qiu的博客
09-17 834
给出一个正整数 k 和一个长度至少为 k 的数字串 s ，对于 s 的每个长度为 k 的连续子串 t （这样的子串共有 ∣s∣−k+1 个，其中 ∣s∣ 为 s 的长度），定义函数 f(t)，询问所有这些 f(t) 的。给出 i ，令 i 号蚯蚓与紧挨其后的一只蚯蚓分离为两个队伍，具体来说，在分离之后， i 号蚯蚓在其中一个队伍的队尾，原本紧挨其后的那一只蚯蚓在另一个队伍的队首，其余蚯蚓保持队伍的前后关系不变。”的一列为“Yes”，表示该测试点的所有蚯蚓的长度均为 1，并且所有询问串 s 的每一位也均为。
【BZOJ4943】【NOI2017】蚯蚓排队（哈希）
小蒟蒻yyb的博客
07-03 398
题面 BZOJ 洛谷 UOJ 题解 记得去年看网络同步赛的时候是一脸懵逼的。 昨天看到zsyzsyzsy做了，今天就看了看。。 这不是HashHashHash的傻逼题吗。。。 呜。。。 开链表存一下前驱后继，因为询问长度不超过505050 那么，总的串的个数不超过50n50n50n，暴力hashhashhash一下然后拿哈希表存一下就好了。 我开始用的map然后TLE飞了...
BZOJ4944: [Noi2017]泳池
weixin_30343157的博客
07-30 136
BZOJ4944: [Noi2017]泳池 题目背景 久莲是个爱玩的女孩子。 暑假终于到了，久莲决定请她的朋友们来游泳，她打算先在她家的私人海滩外圈一块长方形的海域作为游泳场。 然而大海里有着各种各样的危险，有些地方水太深，有些地方有带毒的水母出没。 她想让圈出来的这一块海域都是安全的。 题目描述 经过初步的分析，她把这块海域抽象成了一个底边长为NN米，高为10011001米...
NOI2017天津市选拔赛(TJOI2017)题目及测试数据
07-04
【标题】"NOI2017天津市选拔赛(TJOI2017)题目及测试数据"涉及的是中国信息学奥林匹克竞赛（National Olympiad in Informatics, NOI）的一个地方选拔赛，即2017年天津赛区的比赛。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题解合集（施工中）
待成熟的葡萄
07-28 4741
Day1queueAnalysis先考虑一种最暴力的算法：使用哈希表储存所有长度不超过K=50K=50的子串，合并和分裂时我们只用修改交界处的K2K^2个子串。时间复杂度O(mK2)O(mK^2)。 考虑精细地实现程序：合并时我们只插入原本没有插入过的，也就是严格跨边界的。虽然单次复杂度可能达到O(K2)O(K^2)，但是均摊意义下是可以保证复杂度的： 设势函数Φ\Phi表示哈希表里面的串的个数
洛谷P3825 [NOI2017]游戏（BZOJ4945）
forezxl的博客
01-19 294
2-SAT DFS 洛谷题目传送门 BZOJ题目传送门 每张地图都只有两种车能使用且只能用一辆车，这样的限制满足2-SAT。 设每种地图能使用的两辆车为id[x]和id[x]^1。 先假设没有x地图存在，那么对于一个限制(x,p,y,q)(x,p,y,q)(x,p,y,q)，我们可以进行如下连边： 1.地图x无法放第p种车。此时略过。 2.地图x可以放第p种车，但是地图y无法放第q种...
[NOI 2017]蚯蚓排队
weixin_30421525的博客
07-09 135
Description 题库链接 蚯蚓幼儿园有 \(n\) 只蚯蚓。幼儿园园长神刀手为了管理方便，时常让这些蚯蚓们列队表演。 所有蚯蚓用从 \(1\) 到 \(n\) 的连续正整数编号。每只蚯蚓的长度可以用一个正整数表示，根据入园要求，所有蚯蚓的长度都不超过 \(6\) 。神刀手希望这些蚯蚓排成若干个队伍，初始时，每只蚯蚓各自排成一个仅有一只蚯蚓的队伍，该蚯蚓既在队首，也在队尾。 神刀手将会依次进...
【NOI】2017 蚯蚓排队（BZOJ 4943，LOJ 2303） 模拟+hash
weixin_34278190的博客
05-25 170
【题目】#2303. 「NOI2017」蚯蚓排队 【题意】给定n条长度不超过6的蚯蚓，初始各自在一个队伍。m次操作：1.将i号蚯蚓和j号蚯蚓的队伍合并（保证i为队尾，j为队首）。2.将i号蚯蚓和它后面的蚯蚓分离成两个队。3.询问：给定字符串S和正整数k，求f(每个长度为k的子串)的乘积。其中f(S)定义为蚯蚓在其队伍向后延伸k位组成的字符串等于S的蚯蚓个数。\(n \leq 2*10^5,m \l...
Luogu3823 NOI2017 蚯蚓排队
wennitao的博客
07-07 273
题面：luogu3823 题意：nnn只蚯蚓，每只蚯蚓有一个≤6\le 6≤6的长度，初始时每只蚯蚓一支队伍。给出mmm个操作。 三种操作： 将以iii结尾的队伍和以jjj开头的队伍合并，且iii的队伍在前。 将iii和iii后面一只蚯蚓处断开，分为两支队伍。 定义以第xxx只蚯蚓开始的长度为kkk的数字串为：从xxx开始的连续kkk只蚯蚓，它们的长度数字连成的字符串。定义给定一个字符串sss和...
LOJ#2303. 「NOI2017」蚯蚓排队
a2850363633的博客
06-20 186
$n \leq 200000$的$1 \leq a_i \leq 6$的蚯蚓，有三种操作：让一只队头蚯蚓接在一只队尾蚯蚓后面；让一队蚯蚓从某个蚯蚓后面断成两队；问：给个字符串，问他的。。算了你们直接看题吧 这什么沙雕题QAQ 所有询问的串只有$nk$种，把他们全丢进hash里面就好了。。注意双hash，一个用来当链表一个用来在链表里判重。 复杂度的话，只考虑合并是$nk$的...
NOI2017蚯蚓排队
dianlvgeng7763的博客
10-10 149
原题链接 发现 k<=50 ，在插入和删除时最多会影响不超过 k2 个串，用链表实现插入和删除,然后只需用哈希表维护每个长度不超过k的串的出现次数，哈希的话可以先用比较大的范围的值处理冲突，再映射到1e8的桶里统计。 考虑复杂度。 首先在删除时由于保证了 c<=1000 所以这部分复杂度是O(ck2)的。 插入时，如果插入操作很慢只有可能是连接两个长度不小于k的串，而...
[NOI2017]蚯蚓排队 题解
Hacheylight
01-05 336
题目就是要求我们能够支持三个操作： 1.把两个字符串合并 2.把两个字符串分裂 3.询问有多少一些个串是这么多串中出现了多少次 我们发现数据范围中的K很小，于是我们就考虑k与什么有关，也就是 我们发现第3中操作可以通过哈希表解决，（显而易见的吧） 那怎么维护哈希表？ 我们思考，两个字符串a,ba,ba,b合并会增加哪些哈希值？ 其实就是那些即有aaa的部分又有bbb部分的子串。 这样的串 ...
「NOI2017」泳池
Cold_Chair的博客
06-16 484
传送门. Code: #include<bits/stdc++.h> #define fo(i, x, y) for(int i = x, B = y; i <= B; i ++) #define ff(i, x, y) for(int i = x, B = y; i < B; i ++) #define fd(i, x, y) for(int i = x, ...
BZOJ4944: [Noi2017]泳池（线性递推）
DZYO的博客
05-11 505
传送门 题解： 首先转换问题为最终面积小等于kkk的个数。 注意到最终图形的底边会被禁止的地方分为不同段，每一段的最大面积都不能超过kkk，那么我们记fifif_i为底边长度为iii且他的上方的图形最大面积不超过kkk的概率，然后就是线性递推了。 考虑如何处理fifif_i，这一段底边一定会向上延伸至某个地方jjj被卡住，然后这个jjj位置相当于把这个图形切为两半，满足两边的高度大于jjj...
NOI2017分身术代码c++
最新发布
10-08
NOI2017（全国青少年信息学奥林匹克竞赛）中的“分身术”题目通常涉及计算机科学中的算法和数据结构，特别是动态规划或者贪心策略。在C++中，这类问题可能需要设计一个状态转移方程来解决。分身术的含义可能是指复制...
评论 1
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值