今年NOI确实是在下输了。最近想把当时不会做的题都写一下。

题意

从2到n(500)这些数字中,选若干分给A,若干分给B,满足不存在:A的某个数和B的某个数的GCD不等于1。

对于寿司晚宴这题,标准解答确实是有个神奇的DP。

算法

我们要关注的只是所有的质数。最简单的想法就是枚举A,B各获得哪些数。但是质数的数量实在比较多。然后有个技巧就是将小于\(\sqrt n\)的质数和大于\(\sqrt n\)的数分开处理。这样做的原因是一个数最多只能有一个大于\(\sqrt n\)的质因数。

这样的话,小于\(\sqrt n\)的质数就只有8个了。状态压缩DP或者容斥,就能在\(O(2^8 \times 2^8 n)\)内计算出,A获得的质数集合是x(状态压缩为一整数),B获得的质数集合是y时,有多少种方案,记为f(x, y)。注意,这里的方案并没有计算含有大于\(\sqrt n\)的质因数的数字

接下来,我们要把这些数字也算入答案。奇妙的DP就在这里体现了。

我们只要枚举大于\(\sqrt n\)的那些质数p,一个一个累加到f里,就可以得到最终的答案了。

设g(i, x, y)表示将\(ip\)这个数字分出去后(或者不分给任何人),A获得的质数集合是x,B获得的质数集合是y的方案数。那么,考虑下一个数\((i+1)p\),分给A(分给B同理):

\[g(i + 1, add(x, (i+1)p), y) = g(i, x, y) + f(x, y)
\]

add(x, num)表示将num这个数加进去后,x的变化。


我可能说得不是很清楚,但这毕竟我是打算自己看的。

代码

#include <bits/stdc++.h>
using namespace std; int n, MOD; const int prime[] = {2, 3, 5, 7, 11, 13, 17, 19}; void clear(int (*array)[256]) {
memset(array, 0, sizeof(int) * 256 * 256);
} void copy(int (*src)[256], int (*dest)[256]) {
memcpy(dest, src, sizeof(int) * 256 * 256);
} void add(int &x, int delta) {
if (delta >= MOD) delta -= MOD;
x += delta;
if (x >= MOD) x -= MOD;
} int main() {
scanf("%d%d", &n, &MOD); static int devide[503];
for (int j = 1; j <= n; j++) {
int ret = 0;
int num = j;
for (int i = 0; i < 8; i++) {
int x = prime[i];
while (num % x == 0) {
num /= x;
ret |= 1 << i;
}
}
if (num == 1) devide[j] = ret;
else devide[j] = ret ? -2 : -1;
} static int dp[2][256][256];
int (*cur)[256] = dp[0];
int (*next)[256] = dp[1]; cur[0][0] = 1;
for (int i = 2; i <= n; i++) {
int s = devide[i];
if (s < 0) continue;
copy(cur, next);
for (int a = 0; a < 256; a++)
for (int b = 0; b < 256; b++) {
int &x = cur[a][b];
if (x) {
if (! (s & b)) add(next[a | s][b], x);
if (! (s & a)) add(next[a][b | s], x);
}
}
swap(next, cur);
} static int f[256][256];
for (int i = 23; i <= n; i++) {
if (devide[i] != -1) continue;
clear(f);
for (int j = 1; j * i <= n; j++) {
int s = devide[j];
for (int a = 255; a >= 0; a--)
for (int b = 255; b >= 0; b--) {
if (! (s & b)) add(f[a | s][b], f[a][b] + cur[a][b]);
}
}
for (int a = 0; a < 256; a++)
for (int b = 0; b < 256; b++)
add(cur[a][b], f[a][b] + f[b][a]);
} int ans = 0;
for (int a = 0; a < 256; a++)
for (int b = 0; b < 256; b++)
add(ans, cur[a][b]);
printf("%d\n", ans); return 0;
}

NOI2015 寿司晚宴的更多相关文章

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

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

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

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

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

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

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

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

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

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

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

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

  7. BZOJ 4197: [Noi2015]寿司晚宴 状态压缩 + 01背包

    4197: [Noi2015]寿司晚宴 Time Limit: 10 Sec  Memory Limit: 512 MB Description 为了庆祝 NOI 的成功开幕,主办方为大家准备了一场寿 ...

  8. bzoj 4199 [NOI2015]寿司晚宴

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

  9. [BZOJ]4197: [Noi2015]寿司晚宴

    Time Limit: 10 Sec  Memory Limit: 512 MB Description 为了庆祝 NOI 的成功开幕,主办方为大家准备了一场寿司晚宴.小 G 和小 W 作为参加 NO ...

随机推荐

  1. 在TreeWidget中增加右键菜单功能 以及TreeWidget的基本用法

    TreeWidget 与 TreeView 中实现右键菜单稍有不同, TreeView 中是靠信号与槽 connect(ui->treeView,SIGNAL(customContextMenu ...

  2. WINDOWS硬件通知应用程序的常方法

    摘要:在目前流行的Windows操作系统中,设备驱动程序是操纵硬件的最底层软件接口.为了共享在设备驱动程序设计过程中的经验,给出设备驱动程序通知应用程序的5种方法,详细说明每种方法的原理和实现过程,并 ...

  3. SQL中采用Newtonsoft.Json处理json字符串

    原文 SQL中采用Newtonsoft.Json处理json字符串 使用环境: SQL Server2005; VS2010; 关于SQL中部署CLR程序集的方法,网上一搜一大把,需要了解的自行查阅, ...

  4. Spring-data-redis: 分布式队列

    Redis中list数据结构,具有"双端队列"的特性,同时redis具有持久数据的能力,因此redis实现分布式队列是非常安全可靠的.它类似于JMS中的"Queue&qu ...

  5. TCP/IP笔记 三.运输层(3)——TCP超时重传算法

    TCP 每发送一个报文段,就对这个报文段设置一次计时器.只要计时器设置的重传时间到但还没有收到确认,就要重传这一报文段 1. 平均往返时延RTT 往返时延:一个报文段发出的时间,以及收到相应的确认报文 ...

  6. 关于android源码中的APP编译时引用隐藏的API出现的问题

    今天在编译android源码中的计算器APP时发现,竟然无法使用系统隐藏的API,比如android.os.ServiceManager中的API,引用这个类时提示错误,记忆中在android源码中的 ...

  7. Drying(贪心)

    Drying Time Limit: 2000MS   Memory Limit: 65536K Total Submissions: 11512   Accepted: 2977 Descripti ...

  8. leetCode解题报告5道题(六)

    题目一: Longest Substring Without Repeating Characters Given a string, find the length of the longest s ...

  9. adb shell top

    PID:进程在系统中的ID CPU% - 当前瞬时所以使用CPU占用率 #THR - 程序当前所用的线程数 UID - 运行当前进程的用户id Name - 程序名称android.process.m ...

  10. spring的常用配置

    bean.xml <?xml version="1.0" encoding="UTF-8"?> <beans xmlns="http ...