今年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. [用UpdateLayeredWindow实现任意异形窗口]

    前面提到,我们可以用SetWindowRgn或SetLayeredWindowAttributes实现不规则以及半透明的效果 对于SetWindowRgn,它通过一个Rgn来设置区域,这个Rgn一般可 ...

  2. JAVA FILE or I/O学习 - Desktop本地程序学习

    public class DesktopKnow { public void know() { try { Desktop.getDesktop().open(new File("C:\\P ...

  3. Dancing Stars on Me(判断正多边形)

    Dancing Stars on Me Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 262144/262144 K (Java/Ot ...

  4. 没有login页面

    "/"应用程序中的服务器错误. 无法找到资源. 说明:HTTP 404.您正在查找的资源(或者它的一个依赖项)可能已被移除,或其名称已更改,或暂时不可用.请检查以下 URL 并确保 ...

  5. C# Code Snip

    1.Tryf + TAB+TAB try { } finally { } 2.Prop+Tab+Tab public int MyProperty { get; set; } 3. #region + ...

  6. 简单字符串模式匹配算法的C++实现

    /* * simpleIndex.cpp * Author: Qiang Xiao * Time: 2015-07-13 */ #include<iostream> #include< ...

  7. UNIX/Linux进程间通信IPC---管道--全总结(实例入门)

    管道 一般,进程之间交换信息的方法只能是经由fork或exec传送打开文件,或者通过文件系统.而进程间相互通信还有其他技术——IPC(InterProcessCommunication) (因为不同的 ...

  8. Python之美[从菜鸟到高手]--urlparse源码分析

    urlparse是用来解析url格式的,url格式如下:protocol :// hostname[:port] / path / [;parameters][?query]#fragment,其中; ...

  9. Nginx 之六: Nginx服务器的反向代理功能

    一:Nginx作为正向代理服务器: 1.正向代理:代理(proxy)服务也可以称为是正向代理,指的是将服务器部署在公司的网关,代理公司内部员工上外网的请求,可以起到一定的安全作用和管理限制作用,正向代 ...

  10. python成长之路11

    一.线程: 创建线程有两种方式(本质是一样的,创建好线程之后,cpu调度创建好的线程时执行的其实是Thread的run()方法): import threading def f1(args):prin ...