Gym103960B Fun with Stones

已于 2024-09-13 19:35:57 修改
阅读量996 收藏 14
点赞数 5
分类专栏： dp CodeForces 文章标签： 算法
于 2024-09-08 21:58:13 首次发布
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/m0_53714683/article/details/142033636
版权

并不是多困难或者有趣的题，写sol仅仅是因为觉得好笑（）。

题目大意

三堆石子 Nim 游戏，第 i i i 堆石子数量在 [ l i , r i ] [l_i , r_i] [li​,ri​] 中随机，求先手必胜的概率，对 1 0 9 + 7 10^9+7 109+7 取模。

l i , r i ≤ 1 0 9 l_i , r_i≤10^9 li​,ri​≤109。

题解

说人话就是求三个数异或起来不为0的概率，很容易想到数位dp，由于只有三个数，硬求就好了。

那么为什么我上面写到 “觉得好笑” 呢。

因为听说有个很高妙的写法十几分钟就可以写完，可我是sb我不会，然后很正经地写了个27类巨大分讨6k数位dp，更好笑的是写完居然一遍过了（）。事后知道有高妙写法的我：(ｉДｉ) 。

还是大概说一下，先算出异或和为0的方案数最后减一下除以总数即可。区间比较难求那就dp求出范围分别小于 x , y , z x,y,z x,y,z 的然后容斥一下即可。dp具体实现就是 f i , 0 / 1 , 0 / 1 , 0 / 1 f_{i,0/1,0/1,0/1} fi,0/1,0/1,0/1​ ，代表现在考虑到第 i i i 位、第一个数是否顶着上界、第二个数是否顶着上界、第三个数是否顶着上界的方案数。转移时分别枚举三个数在这一位上是 0/1，三个数的当前位置分别和上界当前位置比较大小再自己算一下这一位可以从哪些状态转移过来即可。

复杂度也许是 l o g log log 带巨大常数，总之能过。

Code

写得很蠢，但是因为很好笑所以还是贴上来吧，就当看个乐子。

#include<bits/stdc++.h>
using namespace std;
 
const int N=35,mod=1e9+7;
 
int l1,r1,l2,r2,l3,r3,ans,f[N][2][2][2];
 
int solve(int x,int a,int a1,int b,int b1,int c,int c1){
	int res=0;
	for(int i=a;i<=a1;i++)
		for(int j=b;j<=b1;j++)
			for(int k=c;k<=c1;k++)
				res=(res+f[x][i][j][k])%mod;
	return res;
}
 
int work(int a,int b,int c){
	memset(f,0,sizeof(f));
	f[32][1][1][1]=1;
	for(int i=31,x,y,z;i>=0;i--){
		x=((a>>i)&1),y=((b>>i)&1),z=((c>>i)&1);
		for(int a1=0;a1<=1;a1++)
			for(int b1=0;b1<=1;b1++)
				for(int c1=0;c1<=1;c1++){
					if((a1^b1^c1)!=0) continue;
					if(a1<x){
						if(b1<y){
							if(c1<z) f[i][0][0][0]=(f[i][0][0][0]+solve(i+1,0,1,0,1,0,1))%mod;
							else if(c1==z){
								f[i][0][0][0]=(f[i][0][0][0]+solve(i+1,0,1,0,1,0,0))%mod;
								f[i][0][0][1]=(f[i][0][0][1]+solve(i+1,0,1,0,1,1,1))%mod;
							}
							else f[i][0][0][0]=(f[i][0][0][0]+solve(i+1,0,1,0,1,0,0))%mod;
						}
						else if(b1==y){
							if(c1<z){
								f[i][0][0][0]=(f[i][0][0][0]+solve(i+1,0,1,0,0,0,1))%mod;
								f[i][0][1][0]=(f[i][0][1][0]+solve(i+1,0,1,1,1,0,1))%mod;
							}
							else if(c1==z){
								f[i][0][0][0]=(f[i][0][0][0]+solve(i+1,0,1,0,0,0,0))%mod;
								f[i][0][1][0]=(f[i][0][1][0]+solve(i+1,0,1,1,1,0,0))%mod;
								f[i][0][0][1]=(f[i][0][0][1]+solve(i+1,0,1,0,0,1,1))%mod;
								f[i][0][1][1]=(f[i][0][1][1]+solve(i+1,0,1,1,1,1,1))%mod;
							}
							else{
								f[i][0][0][0]=(f[i][0][0][0]+solve(i+1,0,1,0,0,0,0))%mod;
								f[i][0][1][0]=(f[i][0][1][0]+solve(i+1,0,1,1,1,0,0))%mod;
							}
						}
						else{
							if(c1<z) f[i][0][0][0]=(f[i][0][0][0]+solve(i+1,0,1,0,0,0,1))%mod;
							else if(c1==z){
								f[i][0][0][0]=(f[i][0][0][0]+solve(i+1,0,1,0,0,0,0))%mod;
								f[i][0][0][1]=(f[i][0][0][1]+solve(i+1,0,1,0,0,1,1))%mod;
							}
							else f[i][0][0][0]=(f[i][0][0][0]+solve(i+1,0,1,0,0,0,0))%mod;
						}
					}
					else if(a1==x){
						if(b1<y){
							if(c1<z){
								f[i][0][0][0]=(f[i][0][0][0]+solve(i+1,0,0,0,1,0,1))%mod;
								f[i][1][0][0]=(f[i][1][0][0]+solve(i+1,1,1,0,1,0,1))%mod;
							}
							else if(c1==z){
								f[i][0][0][0]=(f[i][0][0][0]+solve(i+1,0,0,0,1,0,0))%mod;
								f[i][0][0][1]=(f[i][0][0][1]+solve(i+1,0,0,0,1,1,1))%mod;
								f[i][1][0][0]=(f[i][1][0][0]+solve(i+1,1,1,0,1,0,0))%mod;
								f[i][1][0][1]=(f[i][1][0][1]+solve(i+1,1,1,0,1,1,1))%mod;
							}
							else{
								f[i][0][0][0]=(f[i][0][0][0]+solve(i+1,0,0,0,1,0,0))%mod;
								f[i][1][0][0]=(f[i][1][0][0]+solve(i+1,1,1,0,1,0,0))%mod;
							}
						}
						else if(b1==y){
							if(c1<z){
								f[i][0][0][0]=(f[i][0][0][0]+solve(i+1,0,0,0,0,0,1))%mod;
								f[i][0][1][0]=(f[i][0][1][0]+solve(i+1,0,0,1,1,0,1))%mod;
								f[i][1][0][0]=(f[i][1][0][0]+solve(i+1,1,1,0,0,0,1))%mod;
								f[i][1][1][0]=(f[i][1][1][0]+solve(i+1,1,1,1,1,0,1))%mod;
							}
							else if(c1==z){
								f[i][0][0][0]=(f[i][0][0][0]+solve(i+1,0,0,0,0,0,0))%mod;
								f[i][0][1][0]=(f[i][0][1][0]+solve(i+1,0,0,1,1,0,0))%mod;
								f[i][0][0][1]=(f[i][0][0][1]+solve(i+1,0,0,0,0,1,1))%mod;
								f[i][0][1][1]=(f[i][0][1][1]+solve(i+1,0,0,1,1,1,1))%mod;
								f[i][1][0][0]=(f[i][1][0][0]+solve(i+1,1,1,0,0,0,0))%mod;
								f[i][1][1][0]=(f[i][1][1][0]+solve(i+1,1,1,1,1,0,0))%mod;
								f[i][1][0][1]=(f[i][1][0][1]+solve(i+1,1,1,0,0,1,1))%mod;
								f[i][1][1][1]=(f[i][1][1][1]+solve(i+1,1,1,1,1,1,1))%mod;
							}
							else{
								f[i][0][0][0]=(f[i][0][0][0]+solve(i+1,0,0,0,0,0,0))%mod;
								f[i][0][1][0]=(f[i][0][1][0]+solve(i+1,0,0,1,1,0,0))%mod;
								f[i][1][0][0]=(f[i][1][0][0]+solve(i+1,1,1,0,0,0,0))%mod;
								f[i][1][1][0]=(f[i][1][1][0]+solve(i+1,1,1,1,1,0,0))%mod;
							}
						}
						else{
							if(c1<z){
								f[i][0][0][0]=(f[i][0][0][0]+solve(i+1,0,0,0,0,0,1))%mod;
								f[i][1][0][0]=(f[i][1][0][0]+solve(i+1,1,1,0,0,0,1))%mod;
							}
							else if(c1==z){
								f[i][0][0][0]=(f[i][0][0][0]+solve(i+1,0,0,0,0,0,0))%mod;
								f[i][0][0][1]=(f[i][0][0][1]+solve(i+1,0,0,0,0,1,1))%mod;
								f[i][1][0][0]=(f[i][1][0][0]+solve(i+1,1,1,0,0,0,0))%mod;
								f[i][1][0][1]=(f[i][1][0][1]+solve(i+1,1,1,0,0,1,1))%mod;
							}
							else{
								f[i][0][0][0]=(f[i][0][0][0]+solve(i+1,0,0,0,0,0,0))%mod;
								f[i][1][0][0]=(f[i][1][0][0]+solve(i+1,1,1,0,0,0,0))%mod;
							}
						}
					}
					else{
						if(b1<y){
							if(c1<z) f[i][0][0][0]=(f[i][0][0][0]+solve(i+1,0,0,0,1,0,1))%mod;
							else if(c1==z){
								f[i][0][0][0]=(f[i][0][0][0]+solve(i+1,0,0,0,1,0,0))%mod;
								f[i][0][0][1]=(f[i][0][0][1]+solve(i+1,0,0,0,1,1,1))%mod;
							}
							else f[i][0][0][0]=(f[i][0][0][0]+solve(i+1,0,0,0,1,0,0))%mod;
						}
						else if(b1==y){
							if(c1<z){
								f[i][0][0][0]=(f[i][0][0][0]+solve(i+1,0,0,0,0,0,1))%mod;
								f[i][0][1][0]=(f[i][0][1][0]+solve(i+1,0,0,1,1,0,1))%mod;
							}
							else if(c1==z){
								f[i][0][0][0]=(f[i][0][0][0]+solve(i+1,0,0,0,0,0,0))%mod;
								f[i][0][1][0]=(f[i][0][1][0]+solve(i+1,0,0,1,1,0,0))%mod;
								f[i][0][0][1]=(f[i][0][0][1]+solve(i+1,0,0,0,0,1,1))%mod;
								f[i][0][1][1]=(f[i][0][1][1]+solve(i+1,0,0,1,1,1,1))%mod;
							}
							else{
								f[i][0][0][0]=(f[i][0][0][0]+solve(i+1,0,0,0,0,0,0))%mod;
								f[i][0][1][0]=(f[i][0][1][0]+solve(i+1,0,0,1,1,0,0))%mod;
							}
						}
						else{
							if(c1<z) f[i][0][0][0]=(f[i][0][0][0]+solve(i+1,0,0,0,0,0,1))%mod;
							else if(c1==z){
								f[i][0][0][0]=(f[i][0][0][0]+solve(i+1,0,0,0,0,0,0))%mod;
								f[i][0][0][1]=(f[i][0][0][1]+solve(i+1,0,0,0,0,1,1))%mod;
							}
							else f[i][0][0][0]=(f[i][0][0][0]+solve(i+1,0,0,0,0,0,0))%mod;
						}
					}
				}
	}
	return solve(0,0,1,0,1,0,1);
}
 
int qpow(int x,int y){
	int res=1;
	for(;y;y>>=1,x=1ll*x*x%mod)
		if(y&1) res=1ll*res*x%mod;
	return res;
}
 
int main(){
	scanf("%d%d%d%d%d%d",&l1,&r1,&l2,&r2,&l3,&r3);
	ans=((work(r1,r2,r3)-work(l1-1,r2,r3)-work(r1,l2-1,r3)-work(r1,r2,l3-1)+work(l1-1,l2-1,r3)+work(l1-1,r2,l3-1)+work(r1,l2-1,l3-1)-work(l1-1,l2-1,l3-1))%mod+mod)%mod;
	int az=1ll*(r1-l1+1)*(r2-l2+1)%mod*(r3-l3+1)%mod;
	printf("%d",1ll*(az-ans+mod)%mod*qpow(az,mod-2)%mod);
	return 0;
}

看这幽默的代码不去做大模拟真的可惜了。

pengyixuan2007
关注
专栏目录
Hands-On Intelligent Agents with OpenAI Gym
11-03
Get started with OpenAI Gym and PyTorch for deep reinforcement learning Discover deep Q learning agents to solve discrete optimal control tasks Create custom learning environments for real-world ...
基于gym的q-learning强化学习实践
07-18
b. 在环境中执行这个动作，观察新状态和获得的奖励。 c. 更新Q表，使用Q-learning的更新公式：`Q(s, a) (s, a) + α * (r + γ * max(Q(s', a')) - Q(s, a))`，其中s是旧状态，a是执行的动作，s'是新状态，α是...
gym-0.13.1.tar.gz
07-29
《gym：开启强化学习之旅》 gym，全称为OpenAI Gym，是一个开源的Python库，用于开发和比较强化学习算法。这个“gym-0.13.1.tar.gz”文件是gym库的一个版本0.13.1的压缩包，包含了安装gym所需的所有组件，方便用户...
OpenAI Gym环境.md
02-10
OpenAI Gym是强化学习的一部分，强化学习的环境列出了可以在其上运行算法的Gym环境。通用包装的列出的环境。它们可以对代理与环境之间交换的数据执行预处理/后处理。包含与Gym环境兼容的代理列表。代理有助于针对...
gym-0.12.0所有环境
03-29
gym-0.12.0所有环境，800多个，比想象中的多。python-gym-0.12.0
＜使用生成式AI对四种冒泡排序实现形式分析解释的探讨整理＞
m0_57859086的博客
10-01 1130
/冒泡排序实现1——不常用（这个版本两层for循环分别是i++、j++）// 声明三个整型变量i、j和temp用于循环控制和交换元素for (i = 0;i < n - 1;i++) { // 外层循环，从数组的第一个元素开始遍历到倒数第二个元素j < n;j++) { // 内层循环，从当前外层循环位置的下一个元素开始遍历到最后一个元素if (A[i] > A[j]) { // 如果当前元素大于下一个元素// 将当前元素存储到临时变量temp中// 将下一个元素赋值给当前元素。
排序算法之归并排序
qq_64039411的博客
10-01 511
需要注意的是，在实际应用中，归并排序的空间复杂度较高，需要额外的空间来存储临时数组。如果对空间复杂度有要求，可以考虑使用其他排序算法。归并排序是一种分治算法，它将一个数组分成两个子数组，然后递归地对子数组进行排序，最后将两个有序的子数组合并成一个有序的数组。归并排序的时间复杂度是O(nlogn)，其中n是数组的长度。它是一种稳定的排序算法，适用于各种数据类型的排序。合并排序的关键在于合并操作。
DAY17||654.最大二叉树 |617.合并二叉树 |700.二叉搜索树中的搜索 |
2301_79865280的博客
09-29 1309
给定一个不含重复元素的整数数组。通过给定的数组构建最大二叉树，并且输出这个树的根节点。示例 ：构造二叉树一定要，因为要先找到中，才能分割左右子树。本题图解思路：本题做过上一题构建数组确实变得容易很多。代码优化方向有，不用每次找左右子树的时候构建新的数组的话效率就会高很多。（即直接用下标索引直接在原数组上操作）关于加不加if的问题，如果允许空节点进入递归，就不加。
算法：按既定顺序创建目标数组
AWSDN的博客
09-29 243
本文主要介绍按既定顺序创建目标数组问题
【韩顺平Java笔记】第6章：数组、排序和查找
随手写写
10-02 716
排序是将多个数据，依指定的顺序进行排列的过程。
代码随想录算法训练营day46
qq_56445436的博客
10-01 254
【代码】代码随想录算法训练营day46。
代码随想录打卡Day46
weixin_52151595的博客
09-28 469
今天是动态规划的最后一天，终于要结束折磨了！！！今天的两道题，第一道看视频做的，第二道自己AC的，开心！！！！
35.搜索插入位置
m0_59860403的博客
10-01 205
给定一个排序数组和一个目标值，在数组中找到目标值，并返回其索引。如果目标值不存在于数组中，返回它将会被按顺序插入的位置。请必须使用时间复杂度为O(log n)的算法。214。
TypeScript 算法手册【插入排序】
bobostudio1995的博客
10-01 1717
插入排序是一种简单直观的排序算法，类似于我们整理图书馆的书架，当你面对一排杂乱无章的书籍时，该如何整理它们?我们可能会这样做：从左到右一本一本来，拿起一本书，将它插入到已经按照特定顺序排列好的书籍中的正确位置，这就是插入排序的基本思想。i < len;i++) {j--;return arr;插入排序虽然在大规模数据排序中不如一些高级算法高效，它在某些特定场景下仍然有其独特的优势。当数据规模较小或者数据已经基本有序时，插入排序可能会比一些复杂的排序算法更快。
算法（食物链）
最新发布
yyssas的博客
10-04 446
240. 食物链动物王国中有三类动物 A,B,C𝐴,𝐵,𝐶，这三类动物的食物链构成了有趣的环形。A𝐴 吃 B𝐵，B𝐵 吃 C𝐶，C𝐶 吃 A𝐴。现有 N𝑁 个动物，以 1∼N1∼𝑁 编号。每个动物都是 A,B,C𝐴,𝐵,𝐶 中的一种，但是我们并不知道它到底是哪一种。有人用两种说法对这 N𝑁 个动物所构成的食物链关系进行描述：第一种说法是1 X Y，表示 X𝑋 和 Y𝑌 是同类。第二种说法是2 X Y，表示 X𝑋 吃 Y𝑌。
算法训练营打卡Day18
Zoe_igxgx的博客
09-29 924
示例 2: 输入: root = [3,5,1,6,2,0,8,null,null,7,4], p = 5, q = 4 输出: 5 解释: 节点 5 和节点 4 的最近公共祖先是节点 5。示例 1: 输入: root = [3,5,1,6,2,0,8,null,null,7,4], p = 5, q = 1 输出: 3 解释: 节点 5 和节点 1 的最近公共祖先是节点 3。例如，给定如下二叉树: root = [3,5,1,6,2,0,8,null,null,7,4]
python-鸡尾酒疗法/图像相似度/第n小的质数
2401_83954530的博客
10-02 901
该疗法的应用可以减少单一用药产生的抗药性，最大限度地抑制病毒的复制，使被破坏的机体免疫功能部分甚至全部恢复，从而延缓病程进展，延长患者生命，提高生活质量。假设鸡尾酒疗法的有效率为 x，新疗法的有效率为 y，如果 y−x 大于 5%，则效果更好，如果 x−y 大于 5%，则效果更差，否则称为效果差不多。其余 n 行每行两个整数，第一个整数是临床实验的总病例数（小于等于10000)，第二个疗效有效的病例数。第一行包含两个整数 m 和 n，表示图像的行数和列数，中间用单个空格隔开1≤m≤100,1≤n≤100。
算法笔记（六）——链表
m0_74317866的博客
10-03 387
new head。
【宽搜】2. leetcode 102 二叉树的层序遍历
gaogo
10-04 170
讲解《二叉树的层序遍历》
subprocess-exited-with-error gym
09-15
根据提供的引用内容，你遇到的问题是关于安装gym库时出现的subprocess-exited-with-error错误。根据引用和引用的建议，你可以尝试使用以下命令来安装gym库： pip install gym 如果出现错误提示，你可以尝试使用...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值