[NOI2015][bzoj4197] 寿司晚宴 [状压dp+质因数]
题面
思路
首先,要让两个人选的数字全部互质,那么有一个显然的充要条件:甲选的数字的质因数集合和乙选的数字的质因数集合没有交集
30pt
这种情况下n<=30,也就是说可用的质数只有10个,我们可以开个状压搞一搞
设$dp[S_1][S_2]$表示甲选择的质因数集合是$S_1$,乙是$S_2$的总情况数,
对于每个2-n分解质因数,把每个质因数是否出现状压起来存下来,dp的时候从前往后扫
那么可以刷表法做一波$dp[i][S_1|k][S_2]+=dp[i-1][S_1][S_2]$$(k$ $bitand$ $S_2=0)$,或者$dp[i][S_1][S_2|k]+=dp[i-1][S_1][S_2]$$(k$ $bitand$ $S_1=0)$,其中k是当前数的质因数集合,$bitand$是位与
滚动数组优化一下,把第一维去掉,总效率是$O(2^{20}n)$
100pt
这时有什么变化了呢?没错,n到了500以后可用的质因数变多了,我们无法把它们全部都压进一个数里面了
注意到,一个小于500的数,最多只可能有1个比22大的质因子
所以我们可以把这个质因子单独拿出来记录一下(没有就记为0)
然后,我们把2-n这些数按照大质因子大小排序,这样令大质因子相同的数排在一起(也就是不能甲乙同时选的)
我们记录三个相同数组:$dp[S_1][S_2],f1[][],f2[][]$,因为小质因数只有8个,所以$0\leq S_1,S_2\leq 255$
对于每一段大质因子相同的数,我们在这一段开始的时候把dp的值赋给f1和f2,然后在这一段内部用刷表法推f1和f2
具体来说这么做
$f1[i][S_1|k][S_2]+=f1[i-1][S_1][S_2]$$(k$ $bitand$ $S_2=0)$
或者$f2[i][S_1][S_2|k]+=f2[i-1][S_1][S_2]$$(k$ $bitand$ $S_1=0)$
其中k是当前数的小质因数集合(只有8个二进制位了)
大家应当看出来了,f1表示的就是这个大质因子让第一个人选,f2就是这个大质因子让第二个人选
这一段数推完以后,再把f1f2合并到dp里面,$dp[S_1][S_2]=f1[S_1][S_2]+f2[S_1][S_2]-dp[S_1][S_2]$
这里减掉一个dp是因为两种情况会重复统计两个人都不选的情况(也就是原来的dp[S_1][S_2]的值),减掉即可
最后答案就是$dp[0-255][0-255]$的和了
总时间复杂度为$O(n2^{16})$(居然比30pt做法还快【滑稽】)
Code:
详细的实现见代码
// luogu-judger-enable-o2
#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#define ll long long
using namespace std;
inline int read(){
int re=0,flag=1;char ch=getchar();
while(ch>'9'||ch<'0'){
if(ch=='-') flag=-1;
ch=getchar();
}
while(ch>='0'&&ch<='9') re=(re<<1)+(re<<3)+ch-'0',ch=getchar();
return re*flag;
}
int n,MOD;
int p[10]={0,2,3,5,7,11,13,17,19,0};
struct node{
int val,big,S;//big就是大质因数,S是小质因数集合
void init(){
int i,tmp=val;big=-1;
for(i=1;i<=8;i++){//分解质因数
if(tmp%p[i]) continue;
S|=(1<<i-1);
while(tmp%p[i]==0) tmp/=p[i];
}
if(tmp!=1) big=tmp;
}
}a[510];
inline bool cmp(node l,node r){
return l.big<r.big;
}
int pl(int l,int r){
l+=r;
return l>=MOD?l-MOD:l;
}
int dp[300][300],f1[300][300],f2[300][300];
int main(){
n=read();MOD=read();int i,j,k;
for(i=2;i<=n;i++) a[i-1].val=i,a[i-1].init();
sort(a+1,a+n,cmp);
dp[0][0]=1;
for(i=1;i<n;i++){
if(i==1||a[i].big!=a[i-1].big||a[i].big==-1){
memcpy(f1,dp,sizeof(f1));
memcpy(f2,dp,sizeof(f2));
}
for(j=255;j>=0;j--){
for(k=255;k>=0;k--){//因为是滚动数组,所以一定要倒着推
if(j&k) continue;
if((a[i].S&j)==0) f2[j][k|a[i].S]=pl(f2[j][k|a[i].S],f2[j][k]);
if((a[i].S&k)==0) f1[j|a[i].S][k]=pl(f1[j|a[i].S][k],f1[j][k]);
}
}
if(i==n-1||a[i].big!=a[i+1].big||a[i].big==-1){
for(j=0;j<=255;j++){
for(k=0;k<=255;k++){
if(j&k) continue;
dp[j][k]=pl(f1[j][k],pl(f2[j][k],MOD-dp[j][k]));
}
}
}
}
ll ans=0;
for(j=0;j<=255;j++){
for(k=0;k<=255;k++){
if((j&k)==0&&dp[j][k]) ans=pl(ans,dp[j][k]);
}
}
cout<<ans;
}
[NOI2015][bzoj4197] 寿司晚宴 [状压dp+质因数]的更多相关文章
- BZOJ 4197: [Noi2015]寿司晚宴 状压dp+质因数分解
挺神的一道题 ~ 由于两个人选的数字不能有互质的情况,所以说对于一个质因子来说,如果 1 选了,则 2 不能选任何整除该质因子的数. 然后,我们发现对于 1 ~ 500 的数字来说,只可能有一个大于 ...
- NOI 2015 寿司晚宴 (状压DP+分组背包)
题目大意:两个人从2~n中随意取几个数(不取也算作一种方案),被一个人取过的数不能被另一个人再取.两个人合法的取法是,其中一个人取的任何数必须与另一个人取的每一个数都互质,求所有合法的方案数 (数据范 ...
- 【BZOJ4197】[Noi2015]寿司晚宴 状压DP+分解质因数
[BZOJ4197][Noi2015]寿司晚宴 Description 为了庆祝 NOI 的成功开幕,主办方为大家准备了一场寿司晚宴.小 G 和小 W 作为参加 NOI 的选手,也被邀请参加了寿司晚宴 ...
- bzoj4197 [Noi2015]寿司晚宴——状压DP
题目:https://www.lydsy.com/JudgeOnline/problem.php?id=4197 首先,两个人选的数都互质可以看作是一个人选了一个数,就相当于选了一个质因数集合,另一个 ...
- 【BZOJ-4197】寿司晚宴 状压DP
4197: [Noi2015]寿司晚宴 Time Limit: 10 Sec Memory Limit: 512 MBSubmit: 694 Solved: 440[Submit][Status] ...
- [NOI2015]寿司晚宴 --- 状压DP
[NOI2015]寿司晚宴 题目描述 为了庆祝NOI的成功开幕,主办方为大家准备了一场寿司晚宴. 小G和小W作为参加NOI的选手,也被邀请参加了寿司晚宴. 在晚宴上,主办方为大家提供了n−1种不同的寿 ...
- B4197 [Noi2015]寿司晚宴 状压dp
这个题一开始想到了唯一分解定理,然后状压.但是显然数组开不下,后来想到每个数(n<500)大于19的素因子只可能有一个,所以直接单独存就行了. 然后正常状压dp就很好搞了. 题干: Descri ...
- [NOI2015]寿司晚宴——状压dp
题目转化:将2~n的数分成两组,可以不选,使得这两组没有公共的质因子.求方案数. 选择了一个数,相当于选择了它的所有质因子. 30分: 发现,n<=30的时候,涉及到的质因子也就10个.2,3, ...
- BZOJ 4197 NOI 2015 寿司晚宴 状压DP
4197: [Noi2015]寿司晚宴 Time Limit: 10 Sec Memory Limit: 512 MBSubmit: 694 Solved: 440[Submit][Status] ...
随机推荐
- ReplaceChar
好吧,给个char的,替换单个字符.这样会快一些吧,这个是置换,连长度都不用了 bool ReplaceChar(char *str,const char src, const char dst){ ...
- hdu 2828 Buy Tickets
Buy Tickets Time Limit : 8000/4000ms (Java/Other) Memory Limit : 131072/65536K (Java/Other) Total ...
- phonegap二维码扫描插件
原文出处:http://rensanning.iteye.com/blog/2034026 谈谈我使用这个的体会吧; git地址 https://github.com/wildabeast/Barco ...
- Linux下vim操作的一些使用技巧
以下均为个人在编程时对vim编辑器的一些心得,大神请指点,新手可以看过来 1.多文本编辑 vim -On/-on filename_1 … filename_n 如上所示,在要编辑的文件名前加上“-O ...
- linux 开机自启动 Tomcat
1.修改脚本文件rc.local:vim /etc/rc.d/rc.local 这个脚本是使用者自定的开机启动程序,可以在里面添加想在系统启动之后执行的脚本或者脚本执行命令 2.添加如下内容: exp ...
- CentOS7 配置环境
1.安装CentOS 配置环境 (1)虚拟机中安装CentOS,进入系统使用yum命令不止正常执行…… 原因: 需要设置网卡激活 解决方法: vi /etc/sysconfig/network-scr ...
- [WC2002][洛谷P1578]奶牛浴场
洛谷题解里那个人可真是话多呢. 题目描述 由于John建造了牛场围栏,激起了奶牛的愤怒,奶牛的产奶量急剧减少.为了讨好奶牛,John决定在牛场中建造一个大型浴场.但是John的奶牛有一个奇怪的习惯,每 ...
- hadoop集群之Datenode无法启动解决办法
hadoop集群之Datenode无法启动解决办法 我们在启动hadoop集群的时候,通过jps查看进程,发现namenode RM和Secondary NameNode都有,但datanode没有启 ...
- es6实现简单模板编译
现在有各种框架,其中一个主要模块就是关于template.最火的vue.react等框架,在这一块上也是是下足了功夫.我也想写一个自己的模板编译工具,所以就做了个简单的实现,主要是使用es6的反引号编 ...
- 详解zabbix2.2.2安装部署(Server端篇)
今天开始安装zabbix.zabbix需要LNMP或者LAMP环境.环境的搭建不在本章范围内. LNMP环境配置 Linux安装:http://www.osyunwei.com/archives/10 ...