[NOI2015]寿司晚宴 --- 状压DP
[NOI2015]寿司晚宴
题目描述
为了庆祝NOI的成功开幕,主办方为大家准备了一场寿司晚宴。
小G和小W作为参加NOI的选手,也被邀请参加了寿司晚宴。
在晚宴上,主办方为大家提供了n−1种不同的寿司,编号1,2,3,⋯,n-1,其中第种寿司的美味度为i+1(即寿司的美味度为从2到n)。
现在小G和小W希望每人选一些寿司种类来品尝,他们规定一种品尝方案为不和谐的当且仅当:
小G品尝的寿司种类中存在一种美味度为x的寿司,小W品尝的寿司中存在一种美味度为y的寿司,而x与y不互质。
现在小G和小W希望统计一共有多少种和谐的品尝寿司的方案(对给定的正整数p取模)。
注意一个人可以不吃任何寿司
输入格式:
从文件dinner.in中读入数据。
输入文件的第1行包含2个正整数n,p中间用单个空格隔开,表示共有n种寿司,最终和谐的方案数要对p取模。
输出格式:
输出到文件dinner.out中。
输出一行包含1个整数,表示所求的方案模p的结果。
数据范围:
\(30pt : 2 <= n <= 30\)
\(50pt : 2 <= n <= 100\)
\(70pt : 2 <= n <= 200\)
\(100pt : 2 <= n <= 500\)
\(p <= 1000000000\)
嗯。。。。。。。。。怎么讲呢?虽然自己很快就想出来了,但是感觉自己还是很弱。
周六下午,一看这数据范围,\(O(n^{3})\)的\(DP\)
\(dp(i,j,k)\)表示枚举到第\(i\)个数字,第\(j\)个质数,第\(k\)个。。。。。哪来的三个维???
突然想到大佬提起过状压。
还能压什么?压因子啊。
看了下500以内的素数个数:94个。。。。。。不现实啊
突然想到\(\sqrt{500}=22\)内只有\(2,3,5,7,11,13,17,19\)这\(8\)个质数
\(dp(i,S)\)表示枚举到数字\(i\),因子状态为\(S\)
复杂度\(O(2^8*n)\)。。。。太小了吧。。。。。
不是两个集合吗?加一维:\(dp(i,S1,S2)\)表示枚举到数字\(i\),第一个集合的因子状态\(S1\),第二个\(S2\)
嗯。。。。。。大于\(\sqrt{n}\)的怎么办???
放置了一会去想其他题目,突然懂了大于\(\sqrt{n}\)的因子不会在一个一个数中出现两次,因此可以把每个大因子集合单独拿出来\(DP\)
???做完了???其实思路是差不多完了。
进入正题:
对于第一维根本没有必要保持,可以去掉。
设\(sta[i]\)表示\(i\)的因子状压出来的数字。
小因子转移方程:
\(dp(S1|sta[i],S2)=dp(S1|sta[i],S2)+dp(S1,S2)(sta[i]\;and\;S2==0)\)
\(dp(S1,S2|sta[i])=dp(S1,S2|sta[i])+dp(S1,S2)(sta[i]\;and\;S1==0)\)
而对于大因子,我们需要多开一维来记录大因子允许存在于第几个集合中
即\(f(0/1,S1,S2)\)
转移方程:(以下转移条件同上,略去)
\(f(0,S1|sta[i],S2)=f(0,S1|sta[i],S2)+f(0,S1,S2)\)
\(f(1,S1|sta[i],S2)=f(1,S1|sta[i],S2)+f(1,S1,S2)\)
而最后合并的时候,因为两个\(f\)数组中都考虑了大因子不存在的情况,因此要减去一个,减什么?
还有什么没有考虑大因子呢?
自然是\(dp(S1,S2)\)
大因子DP合并到小因子DP转移方程:
\(dp(S1,S2)=f(0,S1,S2)+f(1,S1,S2)-dp(S1,S2)\)
最后答案即为:
\(\sum_{i=0}^{2^{8}-1} \sum_{j=0}^{2^{8}-1} dp(i,j)\)
时间复杂度:\(O(n*2^{16})\)
补充:
模的是\(10^{10}\)是\(long\;\;long\)范围的,不能开\(int\)
卡常代码,不建议仿生,可以选择参考其他人的
[NOI2015]寿司晚宴 --- 状压DP的更多相关文章
- 【BZOJ4197】[Noi2015]寿司晚宴 状压DP+分解质因数
[BZOJ4197][Noi2015]寿司晚宴 Description 为了庆祝 NOI 的成功开幕,主办方为大家准备了一场寿司晚宴.小 G 和小 W 作为参加 NOI 的选手,也被邀请参加了寿司晚宴 ...
- B4197 [Noi2015]寿司晚宴 状压dp
这个题一开始想到了唯一分解定理,然后状压.但是显然数组开不下,后来想到每个数(n<500)大于19的素因子只可能有一个,所以直接单独存就行了. 然后正常状压dp就很好搞了. 题干: Descri ...
- bzoj4197 [Noi2015]寿司晚宴——状压DP
题目:https://www.lydsy.com/JudgeOnline/problem.php?id=4197 首先,两个人选的数都互质可以看作是一个人选了一个数,就相当于选了一个质因数集合,另一个 ...
- BZOJ 4197: [Noi2015]寿司晚宴 状压dp+质因数分解
挺神的一道题 ~ 由于两个人选的数字不能有互质的情况,所以说对于一个质因子来说,如果 1 选了,则 2 不能选任何整除该质因子的数. 然后,我们发现对于 1 ~ 500 的数字来说,只可能有一个大于 ...
- [NOI2015]寿司晚宴——状压dp
题目转化:将2~n的数分成两组,可以不选,使得这两组没有公共的质因子.求方案数. 选择了一个数,相当于选择了它的所有质因子. 30分: 发现,n<=30的时候,涉及到的质因子也就10个.2,3, ...
- 【BZOJ-4197】寿司晚宴 状压DP
4197: [Noi2015]寿司晚宴 Time Limit: 10 Sec Memory Limit: 512 MBSubmit: 694 Solved: 440[Submit][Status] ...
- BZOJ 4197 NOI 2015 寿司晚宴 状压DP
4197: [Noi2015]寿司晚宴 Time Limit: 10 Sec Memory Limit: 512 MBSubmit: 694 Solved: 440[Submit][Status] ...
- NOI 2015 寿司晚宴 (状压DP+分组背包)
题目大意:两个人从2~n中随意取几个数(不取也算作一种方案),被一个人取过的数不能被另一个人再取.两个人合法的取法是,其中一个人取的任何数必须与另一个人取的每一个数都互质,求所有合法的方案数 (数据范 ...
- [NOI2015][bzoj4197] 寿司晚宴 [状压dp+质因数]
题面 传送门 思路 首先,要让两个人选的数字全部互质,那么有一个显然的充要条件:甲选的数字的质因数集合和乙选的数字的质因数集合没有交集 30pt 这种情况下n<=30,也就是说可用的质数只有10 ...
随机推荐
- 【CodeForces】713 C. Sonya and Problem Wihtout a Legend
[题目]C. Sonya and Problem Wihtout a Legend [题意]给定n个数字,每次操作可以对一个数字±1,求最少操作次数使数列递增.n<=10^5. [算法]动态规划 ...
- Use JPath but not recursively loop a JObject to modify the values.
I am dealing with a Json file, I parsed it into jObject, I have another list which flattened the pro ...
- jQuery 页面加载初始化
jQuery 页面加载初始化的方法有3种 ,页面在加载的时候都会执行脚本,应该没什么区别,主要看习惯吧,本人觉得第二种方法最好,比较简洁. 第一种: $(document).ready(functio ...
- three.js_ "Failed to execute 'texImage2D' on 'WebGLRenderingContext': tainted canvases may not be loded."
这个报错是请求图片跨域了. 1.当我们使用thee.js的时候肯定会碰到利用各种请求去向服务器请求贴图. 2.假设我们获取的是图片在服务器上的路径然后我们用 加载贴图到这里都是没有问题当我们在贴图加载 ...
- WHY学习python?
1.python更容易上手 2.功能库很多,不用重复造轮子 3.能干的事情很多(网站开发,爬虫,自动化运维,数据分析,游戏开发,人工智能) 网站开发:豆瓣,知乎 网站框架:django (姜狗) py ...
- struts集合类型封装
1.list类型封装
- perl6正则 4: before / after 代码断言: <?{}> / <!{}>
<?before> <? befor XXX> 某字符在 xxx 之前 <?after > <?after XXX> 某字符之后有XXX 对应的取反分别 ...
- 14 - 函数参数检测-inspect模块
目录 1 python类型注解 2 函数定义的弊端 3 函数文档 4 函数注解 4.1 annotation属性 5 inspect模块 5.1 常用方法 5.2 signature类 5.3 par ...
- 016 sleep,wait,yield,join区别
1.线程通常有五种状态,创建,就绪,运行.阻塞和死亡状态.2.阻塞的情况又分为三种:(1).等待阻塞:运行的线程执行wait()方法,该线程会释放占用的所有资源,JVM会把该线程放入“等待池”中.进入 ...
- centos7安装完成后的一些配置
1.打开终端 输入 sudo yum -y update 先更新软件包 2.这是输入语言 应用程序->系统工具->设置->区域和语言->+ ->汉语(中国)-> ...