[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 的结果。

输入示例

100 100000000

输出示例

3107203

数据规模及约定

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;
}
时间: 2024-10-25 09:05:58

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

bzoj4197 [Noi2015]寿司晚宴

Description 为了庆祝 NOI 的成功开幕,主办方为大家准备了一场寿司晚宴.小 G 和小 W 作为参加 NOI 的选手,也被邀请参加了寿司晚宴. 在晚宴上,主办方为大家提供了 n?1 种不同的寿司,编号 1,2,3,-,n?1,其中第 i 种寿司的美味度为 i+1 (即寿司的美味度为从 2 到 n). 现在小 G 和小 W 希望每人选一些寿司种类来品尝,他们规定一种品尝方案为不和谐的当且仅当:小 G 品尝的寿司种类中存在一种美味度为 x 的寿司,小 W 品尝的寿司中存在一种美味度为 y

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

[BZOJ4197][Noi2015]寿司晚宴 Description 为了庆祝 NOI 的成功开幕,主办方为大家准备了一场寿司晚宴.小 G 和小 W 作为参加 NOI 的选手,也被邀请参加了寿司晚宴. 在晚宴上,主办方为大家提供了 n−1 种不同的寿司,编号 1,2,3,…,n−1,其中第 i 种寿司的美味度为 i+1 (即寿司的美味度为从 2 到 n). 现在小 G 和小 W 希望每人选一些寿司种类来品尝,他们规定一种品尝方案为不和谐的当且仅当:小 G 品尝的寿司种类中存在一种美味度为 x

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

N^0.5以内的质数只有8个, dp(i, j, k)表示用了前i个大质数(>N^0.5), 2人选的质数(<=N^0.5)集合分别为j, k时的方案数. 转移时考虑当前的大质数p是给哪个人即可. 时间复杂度O(N*2^16) ----------------------------------------------------------------- #include<cstdio> #include<cstring> #include<algorithm&

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

[NOI2015]寿司晚宴 题目描述 为了庆祝NOI的成功开幕,主办方为大家准备了一场寿司晚宴. 小G和小W作为参加NOI的选手,也被邀请参加了寿司晚宴. 在晚宴上,主办方为大家提供了n?1种不同的寿司,编号1,2,3,?,n-1,其中第种寿司的美味度为i+1(即寿司的美味度为从2到n). 现在小G和小W希望每人选一些寿司种类来品尝,他们规定一种品尝方案为不和谐的当且仅当: 小G品尝的寿司种类中存在一种美味度为x的寿司,小W品尝的寿司中存在一种美味度为y的寿司,而x与y不互质. 现在小G和小W希

[NOI2015]寿司晚宴

题目描述 为了庆祝NOI的成功开幕,主办方为大家准备了一场寿司晚宴.小G和小W作为参加NOI的选手,也被邀请参加了寿司晚宴. 在晚宴上,主办方为大家提供了n?1种不同的寿司,编号1,2,3,?,n-1,其中第种寿司的美味度为i+1(即寿司的美味度为从2到n). 现在小G和小W希望每人选一些寿司种类来品尝,他们规定一种品尝方案为不和谐的当且仅当:小G品尝的寿司种类中存在一种美味度为x的寿司,小W品尝的寿司中存在一种美味度为y的寿司,而x与y不互质. 现在小G和小W希望统计一共有多少种和谐的品尝寿司

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

Time Limit: 10 Sec  Memory Limit: 512 MB Description 为了庆祝 NOI 的成功开幕,主办方为大家准备了一场寿司晚宴.小 G 和小 W 作为参加 NOI 的选手,也被邀请参加了寿司晚宴. 在晚宴上,主办方为大家提供了 n?1 种不同的寿司,编号 1,2,3,-,n?1,其中第 i 种寿司的美味度为 i+1 (即寿司的美味度为从 2 到 n). 现在小 G 和小 W 希望每人选一些寿司种类来品尝,他们规定一种品尝方案为不和谐的当且仅当:小 G 品尝

【bzoj4197】[Noi2015]寿司晚宴 dp

因为每个数只有一个大于根号n的质因子,所以我们把每个数拆成一个大于根号n的质因子乘以一个数的形式,对于大于根号n的质因子相同的数,我们放到一起处理 dp[0/1][i][x][y]表示A/B选了当前的大质数,现在枚举到具有当前大质数的第i个数,之前A选中的集合为x,B选中的集合为y的方案数 dp[0/1][0][x][y]=f[i-1][x][y] dp[0][i][x][y]=dp[0][i-1][x][y]+dp[0][i-1][x-S][y] dp[1][i][x][y]=dp[1][i-

【题解】NOI2015寿司晚宴

想好久啊+不敢写啊--但果然人还是应当勇敢自信,只有坚定地去尝试,才会知道最后的结果.1A真的太开心啦,不过好像我的做法还是比较复杂的样子--理解起来应该算是比较容易好懂的类型,大家可以参考一下思路~ 首先我们先考虑一下简单的30分算法:30以内的质数只有十个左右,可以利用状压表示出两个人所选择的集合,再通过寿司转移即可.之后的大数据呢?我们发现不能这样做是因为之后的质数越来越多,状压的空间就开不下了. 这时要注意到一个性质:对于1~n内的每一个数而言,都可以分解成若干个<sqrt(n)的质数之

bzoj 4197: [Noi2015]寿司晚宴【状压dp】

一个数内可能多个的质因数只有小于根号n的,500内这样的数只有8个,所以考虑状压 把2~n的数处理出小于根号500的质因数集压成s,以及大质数p(没有就是1),然后按p排序 根据题目要求,拥有一个质因数的只能给一个人,所以排序后能给一个人的大质数就是一个区间 然后设f[s1][s2]为一人选s1,另一人选s2的方案数,注意这里的s只压了小于根号500的八个质数 设g[0/1][s1][s2]为一人选s1,另一人选s2的,当前枚举的大质数给小G/小W的方案数 正常转移即可 然后注意把g转到f上时应