BZOJ 4197: [Noi2015]寿司晚宴 状压dp+质因数分解
挺神的一道题 ~
由于两个人选的数字不能有互质的情况,所以说对于一个质因子来说,如果 1 选了,则 2 不能选任何整除该质因子的数.
然后,我们发现对于 1 ~ 500 的数字来说,只可能有一个大于 $\sqrt 500$ 的质因子(两个的话乘积就超过 500 了)
而不大于 $\sqrt 500$ 的质因子总数只有 8 个,所以可以对这 8 个质因子状压.
我们先假设所有数字都 $\eqslant 30$,即所有质因子都 $leqslant \sqrt 500$.
定义状态 dp[i][j] 表示第一个人选的质因子集合为 $i$,第二个人选的质因子集合为 $j$.
那么直接更新 $dp[i|sta][j]+=dp[i][j]$ (sta 与 k 无交集)
对于第二个维度同理.
由于我们用滚动数组,所以还要定义 f1[i][j], f2[i][j] 表示第一个人/第二个人 和 第二个人/第一个人 状态下的方案数.
这个转移弄完之后,我们再考虑加入大于 $\sqrt 500$ 的质因子的情况:
令这个质因子为 $x$,那么会有一段数字都包含这个 $x$,显然这些数中只能让 1 个人去选择这个质因子.
我们沿用上面那个 f1[i][j], f2[i][j] 来更新,然后对于这一段的开始位置直接令 f1=f2=dp,其余的正常 dp。
处理完这段的时候还要减去两个数都不选的情况.
最后 dp'[i][j]=f1[i][j]+f2[i][j]-dp[i][j].
code:
#include <bits/stdc++.h>
#define N 703
#define LL long long
#define setIO(s) freopen(s".in","r",stdin)
using namespace std;
int n;
int poww[N];
LL ans,mod;
LL dp[N][N],f1[N][N],f2[N][N];
int p[13]={0,2,3,5,7,11,13,17,19,23};
struct data
{
int num,bi,sta;
}a[N];
bool cmp(data a,data b)
{
return a.bi<b.bi;
}
void init(int num)
{
int i,j,v=num,bi=0;
for(i=1;i<=8;++i)
{
if(v%p[i]==0)
{
a[num].sta|=poww[i];
while(v%p[i]==0) v/=p[i];
}
}
if(v>1) bi=v;
a[num].bi=bi;
a[num].num=num;
}
int main()
{
// setIO("input");
int i,j;
scanf("%d%lld",&n,&mod);
for(i=1;i<=12;++i) poww[i]=(1<<(i-1));
for(i=2;i<=n;++i) init(i);
sort(a+2,a+1+n,cmp);
dp[0][0]=1ll;
for(i=2;i<=n;++i)
{
if(i==2||!a[i].bi||a[i].bi!=a[i-1].bi)
{
memcpy(f1,dp,sizeof(dp));
memcpy(f2,dp,sizeof(dp));
}
for(j=255;j>=0;--j)
{
for(int k=255;k>=0;--k)
{
if((j&a[i].sta)==0) f1[j][k|a[i].sta]=(f1[j][k|a[i].sta]+f1[j][k])%mod;
if((k&a[i].sta)==0) f2[j|a[i].sta][k]=(f2[j|a[i].sta][k]+f2[j][k])%mod;
}
}
if(i==n||!a[i].bi||a[i].bi!=a[i+1].bi)
{
for(j=0;j<=255;++j)
{
for(int k=0;k<=255;++k)
{
dp[j][k]=(f1[j][k]+f2[j][k]-dp[j][k]+mod)%mod;
}
}
}
}
LL ans=0ll;
for(i=0;i<=255;++i) for(j=0;j<=255;++j) if((i&j)==0) (ans+=dp[i][j])%=mod;
printf("%lld\n",ans);
return 0;
}
BZOJ 4197: [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 首先,两个人选的数都互质可以看作是一个人选了一个数,就相当于选了一个质因数集合,另一个 ...
- [NOI2015][bzoj4197] 寿司晚宴 [状压dp+质因数]
题面 传送门 思路 首先,要让两个人选的数字全部互质,那么有一个显然的充要条件:甲选的数字的质因数集合和乙选的数字的质因数集合没有交集 30pt 这种情况下n<=30,也就是说可用的质数只有10 ...
- [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] ...
- BZOJ 4197: [Noi2015]寿司晚宴( dp )
N^0.5以内的质数只有8个, dp(i, j, k)表示用了前i个大质数(>N^0.5), 2人选的质数(<=N^0.5)集合分别为j, k时的方案数. 转移时考虑当前的大质数p是给哪个 ...
- 【BZOJ-4197】寿司晚宴 状压DP
4197: [Noi2015]寿司晚宴 Time Limit: 10 Sec Memory Limit: 512 MBSubmit: 694 Solved: 440[Submit][Status] ...
随机推荐
- 网格搜索与K近邻中更多的超参数
目录 网格搜索与K近邻中更多的超参数 一.knn网格搜索超参寻优 二.更多距离的定义 1.向量空间余弦相似度 2.调整余弦相似度 3.皮尔森相关系数 4.杰卡德相似系数 网格搜索与K近邻中更多的超参数 ...
- 华为云实战开发】5.如何快速创建免费Git代码仓库【华为云技术分享】
1 文章目的 本文主要帮助已经掌握或者想要掌握Git的开发者,如何更好的应用Git,以及更好的将Git与DevCloud结合应用. 2 概述 2.1 版本控制系统介绍 从狭义上来说,版本控制系统是软件 ...
- IDEA远程调试Ambari Server
1.配置端口 Ambari Server默认配置了服务端的debug参数,端口为5005.如果要修改端口,可以在/usr/sbin/ambari_server_main.py文件中对应地方修改,直接改 ...
- ArcGIS Engine中C#开发不能引用ESRI.ArcGIS.AxControls问题
问题:ArcGIS Engine中C#开发不能引用ESRI.ArcGIS.AxControls问题 解决方法:将这里的特定版本改成“False”即可.
- kube-state-metrics 详解
原文:https://mp.weixin.qq.com/s/176eyFBknzdA5wpiJrxDSg 概述 已经有了 cadvisor.heapster.metric-server,几乎容器运行的 ...
- 2019 魔域java面试笔试题 (含面试题解析)
本人5年开发经验.18年年底开始跑路找工作,在互联网寒冬下成功拿到阿里巴巴.今日头条.魔域等公司offer,岗位是Java后端开发,因为发展原因最终选择去了魔域,入职一年时间了,也成为了面试官,之 ...
- RHEL6搭建网络yum源软件仓库
RHEL的更新包只对注册用户生效,所以需要自己手动改成Centos的更新包 一.查看rhel本身的yum安装包 rpm -qa | grep yum 二.卸载这些软件包 rpm -qa | grep ...
- 【转载】C#中double.TryParse方法和double.Parse方法的异同之处
在C#编程过程中,double.TryParse方法和double.Parse方法都可以将字符串string转换为double类型,但两者还是有区别,最重要的区别在于double.TryParse方法 ...
- fileinput 配置项大全,从源码中翻出了很多属性,没那么多时间一一验证,特发出来给大家参考参考
fileinput 配置项大全,从源码中翻出了很多属性,没那么多时间一一验证,特发出来给大家参考参考 fileinput 配置项大全 option 属性名 属性类型 描述说明 默认值 language ...
- I2C总线
PHILIPS公司开发的两线式串行总线 GPIO模拟i2c驱动中有自己的一套传输算法.GPIO模拟I2C是要占用CPU资源的,而用I2C芯片是不占CPU资源的 特点 接口线少,控制方式简单,器件封装形 ...