[UOJ#129][BZOJ4197][Noi2015]寿司晚宴

试题描述

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

在晚宴上,主办方为大家提供了 n−1 种不同的寿司,编号 1,2,3,…,n−1,其中第 i 种寿司的美味度为 i+1 (即寿司的美味度为从 2 到 n)。
现在小 G 和小 W 希望每人选一些寿司种类来品尝,他们规定一种品尝方案为不和谐的当且仅当:小 G 品尝的寿司种类中存在一种美味度为 x 的寿司,小 W 品尝的寿司中存在一种美味度为 y 的寿司,而 x 与 y 不互质。
现在小 G 和小 W 希望统计一共有多少种和谐的品尝寿司的方案(对给定的正整数 p 取模)。注意一个人可以不吃任何寿司。

输入

输入文件的第 1 行包含 2 个正整数 n,p,中间用单个空格隔开,表示共有 n 种寿司,最终和谐的方案数要对 p 取模。

输出

输出一行包含 1 个整数,表示所求的方案模 p 的结果。

输入示例

 

输出示例


数据规模及约定

2≤n≤500

0<p≤1000000000

题解

果然还是第三题最厉害。

根据“互质”这个条件,我们发现每个数可以用它的质因数表示即可。然后发现一个性质:每个数最多有一个质因数大于 sqrt(500);而小于 sqrt(500) 的质数又只有 8 个。

根据以上性质,我们就可以进行状压 dp 了。

首先对 2 到 n 的每个数处理出两个值:big 最大质因数(若这个质因数小于 sqrt(500) 则 big = 1),s 这个数包含较小质因数的集合(易知 s 是一个 8 位的二进制);

现在我们再读一遍题,问题就转化成了给每个数进行三种决定(把它给小 G,给小 W,或都不给),满足小 G、小 W 拿到的数的质因子没有交集;

对于小于 sqrt(500) 部分的质因子显然可以状压解决,对于那些含有大于 sqrt(500) 的质因子的数我们把它们按 big 值分类,那么每一类中的数的主人必须是同一个,或者没有主人。

分析到这里,怎样 dp 应该很明了了吧。

#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <cctype>
#include <algorithm>
using namespace std; const int BufferSize = 1 << 16;
char buffer[BufferSize], *Head, *Tail;
inline char Getchar() {
if(Head == Tail) {
int l = fread(buffer, 1, BufferSize, stdin);
Tail = (Head = buffer) + l;
}
return *Head++;
}
int read() {
int x = 0, f = 1; char c = Getchar();
while(!isdigit(c)){ if(c == '-') f = -1; c = Getchar(); }
while(isdigit(c)){ x = x * 10 + c - '0'; c = Getchar(); }
return x * f;
} #define maxn 510
#define maxs 256 int n, MOD, f[maxn][3][maxs][maxs], prime[] = {2, 3, 5, 7, 11, 13, 17, 19}; struct Num {
int big, s; bool operator < (const Num& t) const { return big < t.big; } void getNum(int x) {
s = 0;
for(int i = 0; i < 8; i++) if(x % prime[i] == 0) {
s |= (1 << i);
while(x % prime[i] == 0) x /= prime[i];
}
big = x;
return ;
}
} ns[maxn]; int main() {
n = read(); MOD = read(); for(int i = 2; i <= n; i++) ns[i].getNum(i);
sort(ns + 2, ns + n + 1);
int all = (1 << 8) - 1;
f[1][0][0][0] = 1;
ns[1].big = 1;
for(int i = 1; i < n; i++)
for(int cho = 0; cho < 3; cho++)
for(int s1 = 0; s1 <= all; s1++)
for(int s2 = 0; s2 <= all; s2++) if(!(s1 & s2) && f[i][cho][s1][s2]) {
int tmp = f[i][cho][s1][s2];
if(ns[i+1].big == 1) {
(f[i+1][0][s1|ns[i+1].s][s2] += tmp) %= MOD;
(f[i+1][0][s1][s2|ns[i+1].s] += tmp) %= MOD;
(f[i+1][0][s1][s2] += tmp) %= MOD;
}
else if(ns[i+1].big == ns[i].big) {
if(cho == 0) {
(f[i+1][0][s1][s2] += tmp) %= MOD;
(f[i+1][1][s1|ns[i+1].s][s2] += tmp) %= MOD;
(f[i+1][2][s1][s2|ns[i+1].s] += tmp) %= MOD;
}
if(cho == 1) (f[i+1][1][s1|ns[i+1].s][s2] += tmp) %= MOD, (f[i+1][1][s1][s2] += tmp) %= MOD;
if(cho == 2) (f[i+1][2][s1][s2|ns[i+1].s] += tmp) %= MOD, (f[i+1][2][s1][s2] += tmp) %= MOD;
}
else {
(f[i+1][0][s1][s2] += tmp) %= MOD;
(f[i+1][1][s1|ns[i+1].s][s2] += tmp) %= MOD;
(f[i+1][2][s1][s2|ns[i+1].s] += tmp) %= MOD;
}
} int ans = 0;
for(int cho = 0; cho < 3; cho++)
for(int s1 = 0; s1 <= all; s1++)
for(int s2 = 0; s2 <= all; s2++)
if(!(s1 & s2)) (ans += f[n][cho][s1][s2]) %= MOD;
printf("%d\n", ans); return 0;
}

[UOJ#129][BZOJ4197][Noi2015]寿司晚宴的更多相关文章

  1. [BZOJ4197][Noi2015]寿司晚宴

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

  2. BZOJ4197 [Noi2015]寿司晚宴 【状压dp】

    题目链接 BZOJ4197 题解 两个人选的数都互质,意味着两个人选择了没有交集的质因子集合 容易想到将两个人所选的质因子集合作为状态\(dp\) \(n\)以内质数很多,但容易发现\(\sqrt{n ...

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

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

  4. UOJ 129/BZOJ 4197 寿司晚宴 状压DP

    //By SiriusRen #include <cstdio> #include <algorithm> using namespace std; ; struct Node ...

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

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

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

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

  7. BZOJ 4197: [Noi2015]寿司晚宴( dp )

    N^0.5以内的质数只有8个, dp(i, j, k)表示用了前i个大质数(>N^0.5), 2人选的质数(<=N^0.5)集合分别为j, k时的方案数. 转移时考虑当前的大质数p是给哪个 ...

  8. BZOJ_4197_[Noi2015]寿司晚宴_状态压缩动态规划

    BZOJ_4197_[Noi2015]寿司晚宴_状态压缩动态规划 Description 为了庆祝 NOI 的成功开幕,主办方为大家准备了一场寿司晚宴.小 G 和小 W 作为参加 NOI 的选手,也被 ...

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

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

随机推荐

  1. java代码(处理json串)

    package test; import com.alibaba.fastjson.JSON; import com.alibaba.fastjson.JSONObject; public class ...

  2. uiviewcontroller 键盘不遮挡信息

    //添加监听事件 [[NSNotificationCenter defaultCenter] addObserver:self selector:@selector(keyboardWillShow: ...

  3. 快学UiAutomator UiDevice API 详解

    一.按键使用 返回值 方法名 说明 boolean pressBack() 模拟短按返回back键 boolean pressDPadCenter() 模拟按轨迹球中点按键 boolean press ...

  4. python之list [ 列表 ]

    1. 列表是什么? list [ ] 逗号隔开 是一个容器 可以存放任意类型 列表 == 书包 书包里可以放水杯.衣服.袜子.钱包 钱包里可以放钱.身份证件,可以包套包 2. 列表能干什么? 存储大量 ...

  5. _IO_FILE

    hctf2017的babyprintf解法是house of orange,深入学习了一下,牵扯出许多知识,这里先进行第一步:_IO_FILE结构 0x00 _IO_FILE glibc-2.2.1\ ...

  6. 总结:JavaScript异步、事件循环与消息队列、微任务与宏任务

    本人正在努力学习前端,内容仅供参考.由于各种原因(不喜欢博客园的UI),大家可以移步我的github阅读体验更佳:传送门,喜欢就点个star咯,或者我的博客:https://blog.tangzhen ...

  7. (60)zabbix网络发现介绍Network Discovery

    网络发现简介 网络发现有什么用?网络发现怎么配置? 我们带着这两个问题开始我们的网络发现之旅. 比如小明有100台服务器,不想一台台主机去添加,能不能让zabbix自动添加主机呢,当然可以,网络发现便 ...

  8. Controller View 模式

    参考:https://blog.andrewray.me/the-reactjs-controller-view-pattern/ Flux参考:http://www.cnblogs.com/hell ...

  9. perl学习之:肯定匹配和否定匹配

    tr/ / / 替换操作符不支持正则表达式 也不具备双引号替换能力m/ /  s/ / / 都支持正则表达式,并且可以提供或限制双引号替换能力 $string = "25abc8" ...

  10. html块级元素和行级元素的区别和使用

    行内.块状元素区别: 1.行内元素与块级函数可以相互转换,通过修改display属性值来切换块级元素和行内元素,行内元素display:inline,块级元素display:block. 2.行内元素 ...