【题解】NOI2015寿司晚宴

  想好久啊+不敢写啊……但果然人还是应当勇敢自信,只有坚定地去尝试,才会知道最后的结果。1A真的太开心啦,不过好像我的做法还是比较复杂的样子……理解起来应该算是比较容易好懂的类型,大家可以参考一下思路~

  首先我们先考虑一下简单的30分算法:30以内的质数只有十个左右,可以利用状压表示出两个人所选择的集合,再通过寿司转移即可。之后的大数据呢?我们发现不能这样做是因为之后的质数越来越多,状压的空间就开不下了。

  这时要注意到一个性质:对于1~n内的每一个数而言,都可以分解成若干个<sqrt(n)的质数之积 || 在此基础之上再乘上一个 > sqrt(n)的质数。这说明什么?对于所有的>sqrt(n)的质数而言,我们选择一个寿司只可能选择其中的一个——换句话说,就是不同的大质数之间的决策是相互独立的。

  于是就有了如下算法:既然不同的大质数之间不会互相影响,我们就一个一个大质数来统计,之后再累加到一起即可。于是我们增加一维的状态,单独表示这一个大质数。0表示两个集合中均不含有这个大质数因子,1表示第一个人所选择的集合中含这个因子,2表示第二个人选择的集合中含有这个因子。不同的因子之间的转移将所有1&2的状态都加入0并清空1&2即可(对于新的质数来说,之前没有作出过相应的决策,所以是不含有该因子的)。

  网上代码很短,然而我莫名长……

#include <bits/stdc++.h>
using namespace std;
#define maxn 1000
#define maxt 260
#define int long long
int n, mod, S[maxn], CNST = (1 << 8) - 1;
int cnt, dp[3][maxt][maxt], num[maxn], mark[maxn];
int tot, P[maxn], cnp = 8, ans;
int pri[maxn] = {0, 2, 3, 5, 7, 11, 13, 17, 19};
map <int, int> Map;

int read()
{
    int x = 0;
    char c;
    c = getchar();
    while(c < ‘0‘ || c > ‘9‘) c = getchar();
    while(c >= ‘0‘ && c <= ‘9‘) x = x * 10 + c - ‘0‘, c = getchar();
    return x;
}

int up(int &x, int y)
{
    x += y;
    if(x >= mod) x -= mod;
}

void work(int i)
{
    for(int x = CNST; ~x; x --)
        for(int y = CNST; ~y; y --)
        {
            if(x & y) continue;
            int a = dp[1][x][y], b = dp[2][x][y];
            up(dp[1][x | num[i]][y], dp[0][x][y]);
            up(dp[2][x][y | num[i]], dp[0][x][y]);

            up(dp[1][x | num[i]][y], a);
            up(dp[2][x][y | num[i]], b);
        }
}

void DP(int k)
{
    for(int x = 0; x <= CNST; x ++)
        for(int y = 0; y <= CNST; y ++)
            if(x & y) continue;
            else
            {
                up(dp[0][x][y], dp[1][x][y]);
                up(dp[0][x][y], dp[2][x][y]);
                dp[1][x][y] = dp[2][x][y] = 0;
            }
    if(k) for(int i = k; i <= n + 1; i += k) work(i);
    else
    {
        for(int i = 1; i <= cnt; i ++)
            for(int x = CNST; ~x; x --)
                for(int y = CNST; ~y; y --)
                {
                    if(x & y) continue;
                    int k = S[i];
                    int a = dp[0][x][y];
                    up(dp[0][x | num[k]][y], a);
                    up(dp[0][x][y | num[k]], a);
                }
    }
}

signed main()
{
    n = read() - 1, mod = read();
    dp[0][0][0] = 1;
    for(int i = 2; i <= n + 1; i ++)
    {
        int k = i;
        for(int j = 1; j <= cnp; j ++)
        {
            if(!(k % pri[j])) num[i] |= (1 << (j - 1));
            while(!(k % pri[j])) k /= pri[j];
        }
        if(k != 1 && k != 0)
        {
            mark[i - 1] = k;
            if(!Map[k]) Map[k] = 1, P[++ tot] = k;
        }
        else S[++ cnt] = i;
    }
    for(int i = 0; i <= tot; i ++)
        DP(P[i]);
    for(int i = 0; i <= CNST; i ++)
        for(int j = 0; j <= CNST; j ++)
            if(i & j) continue;
            else
            {
                up(ans, dp[0][i][j]);
                up(ans, dp[1][i][j]);
                up(ans, dp[2][i][j]);
            }
    printf("%lld\n", ans);
    return 0;
}

原文地址:https://www.cnblogs.com/twilight-sx/p/8969685.html

时间: 2024-10-01 18:59:34

【题解】NOI2015寿司晚宴的相关文章

【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希

[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 的寿司,而

[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希望统计一共有多少种和谐的品尝寿司

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

[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-

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上时应