[NOI2015]寿司晚宴——状压dp
题目转化:将2~n的数分成两组,可以不选,使得这两组没有公共的质因子。求方案数。
选择了一个数,相当于选择了它的所有质因子。
30分:
发现,n<=30的时候,涉及到的质因子也就10个。2,3,5,7,11,13,19,23,29
直接状压。f[i][A][B] 前i个数,第一个人的质因子选择状态A,第二人B,的方案数。(第一维可以滚动,当然,可以倒序循环直接省略)
每个数质因数分解,前八个质因子,压成二进制数,转移直接按位或。
100分:
质因子太多状压不了。
公理:一个数>=sqrt(n)的质因子最多只有一个。(显然啊,否则质因数分解就超过了n)
所以,对于最大质因子>=sqrt(n)的数,根据最大质因子分类,质因子相同的,要么都考虑进入A,要么都考虑进入B,要么都不选。
设g[0/1][A][B]表示这一阶段,数字都考虑进入A/B的方案数。初值和f一样。
最后,f[A][B]=g[0][A][B]+g[1][A][B]-f[A][B] 因为,所有数都不选的方案计算了两次,所以减一次。
对于最大质因子<sqrt(n)的数,直接根据上面的转移就好。
统计答案时,选择A&B=0 的方案数。
具体看代码:(luogu要开O2)
#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int N=+;
int n,mod;
struct node{
int has;
int z;
}num[N];
bool cmp(node a,node b){
return a.z<b.z;
}
ll ans;
int zh[N],cnt;
int ps[N],tot;
int vis[N];
bool in[N];
void sieve(){
for(int i=;i<=n;i++){
if(!vis[i]){
tot++;ps[tot]=i;
vis[i]=i;
}
for(int j=;j<=tot;j++){
if(i*ps[j]>n) break;
vis[i*ps[j]]=ps[j];
if(i%ps[j]==) break;
}
}
}
void sol(int x){
cnt=;
int kk=x;
for(int i=;i<=tot;i++){
if(x%ps[i]==){
zh[++cnt]=i;//prime's id
while(x%ps[i]==) x/=ps[i];
}
}
for(int i=;i<=cnt;i++){
if(zh[i]>) num[kk].z=zh[i];
else num[kk].has|=(<<zh[i]-);
}
if(!num[kk].z) num[kk].z=;
}
ll f[<<][<<];
ll g[][<<][<<];
int main()
{
scanf("%d%d",&n,&mod);
sieve();
for(int i=;i<=n;i++)sol(i);
sort(num+,num+n+,cmp);
f[][]=;
for(int i=;i<=n;i++){
if(num[i].z==){
for(int A=(<<)-;A>=;A--)
for(int B=(<<)-;B>=;B--){
f[A|num[i].has][B]=(f[A|num[i].has][B]+f[A][B])%mod;
f[A][B|num[i].has]=(f[A][B|num[i].has]+f[A][B])%mod;
}
}
else{
int zz=num[i].z;
memcpy(g[],f,sizeof f);
memcpy(g[],f,sizeof f);
while(num[i].z==zz&&i<=n){
for(int A=(<<)-;A>=;A--)
for(int B=(<<)-;B>=;B--){
g[][A|num[i].has][B]=(g[][A|num[i].has][B]+g[][A][B])%mod;
g[][A][B|num[i].has]=(g[][A][B|num[i].has]+g[][A][B])%mod;
}
i++;
}
i--;
for(int A=(<<)-;A>=;A--)
for(int B=(<<)-;B>=;B--){
f[A][B]=((g[][A][B]+g[][A][B])%mod-f[A][B])%mod;if(f[A][B]<) f[A][B]+=mod;
}
}
}
for(int A=;A<=(<<)-;A++)
for(int B=;B<=(<<)-;B++){
if((A&B)==) ans=(ans+f[A][B])%mod;
}
printf("%lld",ans);
return ;
}
总结:
这个题目其实很巧妙。
思路:
1.互质即没有公共质因子,所以选择一个数就是选择一些质因子。转化题意。
2.30分,发现质因子很小,可以直接状压。
3.100分,发现每个数>=sqrt的质因子最多只有一个。暴力分组考虑,每组内统计都放进A/B的方案。
关键的突破口还是题意的转化。以及第三点。
[NOI2015]寿司晚宴——状压dp的更多相关文章
- [NOI2015]寿司晚宴 --- 状压DP
[NOI2015]寿司晚宴 题目描述 为了庆祝NOI的成功开幕,主办方为大家准备了一场寿司晚宴. 小G和小W作为参加NOI的选手,也被邀请参加了寿司晚宴. 在晚宴上,主办方为大家提供了n−1种不同的寿 ...
- 【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 的数字来说,只可能有一个大于 ...
- 【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 ...
随机推荐
- linux下日志文件error监控报警脚本分享
即对日志文件中的error进行监控,当日志文件中出现error关键字时,即可报警!(grep -i error 不区分大小写进行搜索"error"关键字,但是会将包含error大小 ...
- Linux系统本地yum源环境配置记录
由于IDC的一些服务器没有外网,不能对外访问.所以打算部署一套内网的yum源环境,以供内网服务器使用.以下简单记录下操作过程: 1)下载centos6.9和centos7.3的镜像,并挂载 [root ...
- Python-元组-10
元祖 Why:对于容器型数据类型list,无论谁都可以对其增删改查,那么有一些重要的数据放在list中是不安全的,所以需要一种容器类的数据类型存放重要的数据,创建之初只能查看而不能增删改,这种数据类型 ...
- c++入门之关于cin,cout以及数据流的认识
- 个人作业Week2-代码复审(修改明确了要求)
代码复审 零,说在前面的话 大家完成了个人项目之后,都写了很多代码. 这些代码可能: 大括号换行/不换行 使用tab缩进/使用空格缩进 变量名函数名的定义很好/不好 每个函数都有详细的注释解释函数的功 ...
- SCRUM 12.18
明天就是编译课设的第二次中期考核了,大家都感到有一些压力. 所以我们决定今天减少一些工作量. 工作任务分配依旧如往常 成员 任务 彭林江 落实API 郝倩 研究遍历美团数据方法 牛强 落实意见反馈功能 ...
- beta阶段测试基本概况报告
文件地址 测试基本信息 Bitmap 测试 ...
- 《Linux内核设计与分析》第六周读书笔记——第三章
<Linux内核设计与实现>第六周读书笔记——第三章 20135301张忻估算学习时间:共2.5小时读书:2.0代码:0作业:0博客:0.5实际学习时间:共3.0小时读书:2.0代码:0作 ...
- 实训十二(stick的设定)
上篇我们介绍到人物主角的设定,其实人物是有工具使的,那就是——stick小棍. 信息的获取.起始位置.长度的加载.边界的判断.位置.长度重置是需要我们主要考虑的问题 信息获取上考虑的使什么时候加载st ...
- [福大软工] Z班 第11次成绩排行榜
注:本次成绩排行榜是针对团队Alpha阶段冲刺博客的得分统计 作业要求 http://www.cnblogs.com/easteast/p/7749181.html 评分细则 在团队Alpha阶段刚开 ...