题目转化:将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的更多相关文章

  1. [NOI2015]寿司晚宴 --- 状压DP

    [NOI2015]寿司晚宴 题目描述 为了庆祝NOI的成功开幕,主办方为大家准备了一场寿司晚宴. 小G和小W作为参加NOI的选手,也被邀请参加了寿司晚宴. 在晚宴上,主办方为大家提供了n−1种不同的寿 ...

  2. 【BZOJ4197】[Noi2015]寿司晚宴 状压DP+分解质因数

    [BZOJ4197][Noi2015]寿司晚宴 Description 为了庆祝 NOI 的成功开幕,主办方为大家准备了一场寿司晚宴.小 G 和小 W 作为参加 NOI 的选手,也被邀请参加了寿司晚宴 ...

  3. B4197 [Noi2015]寿司晚宴 状压dp

    这个题一开始想到了唯一分解定理,然后状压.但是显然数组开不下,后来想到每个数(n<500)大于19的素因子只可能有一个,所以直接单独存就行了. 然后正常状压dp就很好搞了. 题干: Descri ...

  4. bzoj4197 [Noi2015]寿司晚宴——状压DP

    题目:https://www.lydsy.com/JudgeOnline/problem.php?id=4197 首先,两个人选的数都互质可以看作是一个人选了一个数,就相当于选了一个质因数集合,另一个 ...

  5. BZOJ 4197: [Noi2015]寿司晚宴 状压dp+质因数分解

    挺神的一道题 ~ 由于两个人选的数字不能有互质的情况,所以说对于一个质因子来说,如果 1 选了,则 2 不能选任何整除该质因子的数. 然后,我们发现对于 1 ~ 500 的数字来说,只可能有一个大于 ...

  6. 【BZOJ-4197】寿司晚宴 状压DP

    4197: [Noi2015]寿司晚宴 Time Limit: 10 Sec  Memory Limit: 512 MBSubmit: 694  Solved: 440[Submit][Status] ...

  7. BZOJ 4197 NOI 2015 寿司晚宴 状压DP

    4197: [Noi2015]寿司晚宴 Time Limit: 10 Sec  Memory Limit: 512 MBSubmit: 694  Solved: 440[Submit][Status] ...

  8. NOI 2015 寿司晚宴 (状压DP+分组背包)

    题目大意:两个人从2~n中随意取几个数(不取也算作一种方案),被一个人取过的数不能被另一个人再取.两个人合法的取法是,其中一个人取的任何数必须与另一个人取的每一个数都互质,求所有合法的方案数 (数据范 ...

  9. [NOI2015][bzoj4197] 寿司晚宴 [状压dp+质因数]

    题面 传送门 思路 首先,要让两个人选的数字全部互质,那么有一个显然的充要条件:甲选的数字的质因数集合和乙选的数字的质因数集合没有交集 30pt 这种情况下n<=30,也就是说可用的质数只有10 ...

随机推荐

  1. Nginx反向代理的简单实现

    1)nginx的反向代理:proxy_pass2)nginx的负载均衡:upstream 下面是nginx的反向代理和负载均衡的实例: 负载机:A机器:103.110.186.8/192.168.1. ...

  2. cf946d 怎样逃最多的课dp

    来源:codeforces                                              D. Timetable Ivan is a student at Berland ...

  3. 【CV】ICCV2015_Describing Videos by Exploiting Temporal Structure

    Describing Videos by Exploiting Temporal Structure Note here: it's a learning note on the topic of v ...

  4. 《Linux内核设计与实现》课本第四章学习总结

    进程调度 4.1 多任务 多任务操作系统就是能同时并发的交互执行多个进程的操作系统. 多任务系统分为两种: 抢占式多任务:Linux提供了抢占式的多任务模式,由调度程序来决定什么时候停止一个进程的运行 ...

  5. Android 學習之旅!(1)

    就這樣就過去了一年加一個學期,現在是大二第二個學期而且是下半學期了,以前都是無所事事,沒事睡睡覺,打打遊戲就過去了,但是想到家境和以後的路,我還是決心自己找點東西學習下,以後出去還能有一技之長(雖然可 ...

  6. Java 一维数组作为参数和返回值

    一维数组作为参数: 传数组的引用 创建数组直接传,本质也是传数组的引用 传null public class Test { //数组作为参数时,可以传递3中形式 public void m1(int[ ...

  7. Alpha 答辩总结模板

    Alpha 答辩总结模板 每个小组提供一篇总结博客(组内共享,每个人都发布),包含: 本组(组名)所有成员(短学号,名,标注组长)(1分) 组内各成员贡献比例,如不提供,取平均分后组长得分减50% G ...

  8. vue容易混淆的点小记

    computed.methods及watch函数的差异 computed:基于依赖进行缓存,若依赖不变,则直接调用缓存(适用于性能开销比较大的时候) methods: 不管数据是否变更,都会进行计算( ...

  9. Eclipse 下载安装

    https://blog.csdn.net/qq_27200591/article/details/72823816(安装eclipse  copy) 第一步:下载eclipse,并安装. 下载链接: ...

  10. Node url模块

    const url = require("url");var u = "https://www.jd.com:443/ad/index?uname=qd";// ...