bzoj4197 [Noi2015]寿司晚宴（数学+状压dp）

最新推荐文章于 2020-03-15 13:19:10 发布
最新推荐文章于 2020-03-15 13:19:10 发布
阅读量189 收藏 1
点赞数
分类专栏： bzoj ---------数学--------- 素数 ---------dp----------- 状压dp
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/Icefox_zhx/article/details/80800654
版权

首先考虑n<=30怎么做，没几个质因子，状压掉
然后n<=100，状压着有点费劲，各种优化，比如*2>n的质因子肯定是贡献3倍答案，*3>n的质因子可以讨论一下也不压进状态里，然后就可以卡过去了
这启发了我们正解的思路！对于大质因子我们一起转移，然后就把这个大质因子这一位扔掉。
我们只需要状压< 500−−−√ 500 的那么8个质因子就好了。
然后复杂度 O(38n) O ( 3 8 n ) 

#include <bits/stdc++.h>
using namespace std;
#define ll long long
#define inf 0x3f3f3f3f
#define N 510
#define k1 1117
#define eps 1e-8
inline char gc(){
    static char buf[1<<16],*S,*T;
    if(S==T){T=(S=buf)+fread(buf,1,1<<16,stdin);if(T==S) return EOF;}
    return *S++;
}
inline ll read(){
    ll x=0,f=1;char ch=gc();
    while(ch<'0'||ch>'9'){if(ch=='-') f=-1;ch=gc();}
    while(ch>='0'&&ch<='9') x=x*10+ch-'0',ch=gc();
    return x*f;
}
int n,prime[N],tot=0;
ll mod;
bool notprime[N];
inline ll mul(ll x,ll y){ll res=x*y-(ll)((long double)x/mod*y+eps)*mod;res%=mod;return res<0?res+mod:res;}
inline ll ksm(ll x,int k){
    ll res=1;for(;k;k>>=1,x=mul(x,x)) if(k&1) res=mul(res,x);return res;
}
inline void getprime(){
    notprime[1]=1;
    for(int i=2;i<=n;++i){
        if(!notprime[i]) prime[++tot]=i;
        for(int j=1;prime[j]*i<=n;++j){
            notprime[prime[j]*i]=1;
            if(i%prime[j]==0) break;
        }
    }
}
namespace sol1{
    bool used[N];
    struct Hash_table{
        ll val;int key1,key2,next;
    };
    inline void inc(ll &x,ll y){x+=y;x%=mod;}
    struct Icefox{
        Hash_table data[20000];
        int h[k1][k1],num;
        inline void ins(int key1,int key2,ll val){
            int x=key1%k1,y=key2%k1;
            for(int i=h[x][y];i;i=data[i].next)
                if(data[i].key1==key1&&data[i].key2==key2){inc(data[i].val,val);return;}
            data[++num].key1=key1;data[num].key2=key2;data[num].val=val;data[num].next=h[x][y];h[x][y]=num;
        }inline void clear(){
            for(int i=1;i<=num;++i){
                int x=data[i].key1%k1,y=data[i].key2%k1;h[x][y]=0;
            }num=0;
        }
    }f[2];
    int a[N],bin[20];
    inline void gao(){
        bin[0]=1;int m=0,p=0;
        for(int i=1;i<=tot;++i) if(prime[i]*2<=n) ++m;
        for(int i=1;i<=m;++i) bin[i]=bin[i-1]<<1;
        for(int i=2;i<=n;++i)
            for(int j=1;j<=m;++j)
                if(i%prime[j]==0) a[i]|=bin[j-1];
        f[p].ins(0,0,1);
        for(int j=m;j>=1;--j){
            for(int x=prime[j];x<=n;x+=prime[j]){
                if(used[x]) continue;used[x]=1;f[p^1].clear();
                for(int i=1;i<=f[p].num;++i){
                    int s1=f[p].data[i].key1,s2=f[p].data[i].key2;ll val=f[p].data[i].val;
                    if((a[x]&s1)==0) f[p^1].ins(s1,s2|a[x],val);
                    if((a[x]&s2)==0) f[p^1].ins(s1|a[x],s2,val);
                    f[p^1].ins(s1,s2,val);
                }p^=1;
            }
        }ll ans=0;
        for(int i=1;i<=f[p].num;++i) inc(ans,f[p].data[i].val);
        printf("%lld\n",mul(ans,ksm(3,tot-m)));
    }
}
namespace sol2{
    bool used[N],mk[N];
    struct Hash_table{
        ll val;int key1,key2,next;
    };
    inline void inc(ll &x,ll y){x+=y;x%=mod;}
    struct Icefox{
        Hash_table data[180000];
        int h[k1][k1],num;
        inline void ins(int key1,int key2,ll val){
            int x=key1%k1,y=key2%k1;
            for(int i=h[x][y];i;i=data[i].next)
                if(data[i].key1==key1&&data[i].key2==key2){inc(data[i].val,val);return;}
            data[++num].key1=key1;data[num].key2=key2;data[num].val=val;data[num].next=h[x][y];h[x][y]=num;
        }inline void clear(){
            for(int i=1;i<=num;++i){
                int x=data[i].key1%k1,y=data[i].key2%k1;h[x][y]=0;
            }num=0;
        }
    }f[2];
    int a[N],bin[20];
    inline void gao(){
        bin[0]=1;int m=0,p=0,n1=0,n2=0;
        for(int i=1;i<=tot;++i){
            if(prime[i]*3<=n) ++m;
            else if(prime[i]*2<=n) mk[prime[i]*2]=1,++n1;
            else ++n2;
        }
        for(int i=1;i<=m;++i) bin[i]=bin[i-1]<<1;
        for(int i=2;i<=n;++i)
            for(int j=1;j<=m;++j)
                if(i%prime[j]==0) a[i]|=bin[j-1];
        f[p].ins(0,0,1);
        for(int j=m;j>=1;--j){
            for(int x=prime[j];x<=n;x+=prime[j]){
                if(used[x]) continue;used[x]=1;f[p^1].clear();
                for(int i=1;i<=f[p].num;++i){
                    int s1=f[p].data[i].key1,s2=f[p].data[i].key2;ll val=f[p].data[i].val;
                    if((a[x]&s1)==0) f[p^1].ins(s1,s2|a[x],val*(mk[x]?2:1));
                    if((a[x]&s2)==0) f[p^1].ins(s1|a[x],s2,val*(mk[x]?2:1));
                    f[p^1].ins(s1,s2,val*(mk[x]?3:1));
                }p^=1;
            }
        }ll ans=0;
        for(int i=1;i<=f[p].num;++i) inc(ans,f[p].data[i].val);
        printf("%lld\n",mul(ans,ksm(3,n2)));
    }
}
namespace sol3{
    int a[N],bin[10];ll f[2][520][520];
    bool used[N];
    inline void inc(ll &x,ll y){x+=y;x%=mod;}
    inline void gao(){
        bin[0]=1;int m=8,p=0;
        for(int i=1;i<=m+1;++i) bin[i]=bin[i-1]<<1;
        for(int i=2;i<=n;++i)
            for(int j=1;j<=m;++j)
                if(i%prime[j]==0) a[i]|=bin[j-1];
        f[p][0][0]=1;
        for(int j=tot;j>8;--j){
            for(int i=1;;i++){
                if(prime[j]*i>n) break;int x=prime[j]*i;used[x]=1;memset(f[p^1],0,sizeof(f[p^1]));a[x]|=bin[m];
                for(int s1=0;s1<bin[m+1];++s1){
                    int st=bin[m+1]-1-s1;
                    for(int s2=st;s2;s2=(s2-1)&st){
                        if((a[x]&s1)==0) inc(f[p^1][s1][s2|a[x]],f[p][s1][s2]);
                        if((a[x]&s2)==0) inc(f[p^1][s1|a[x]][s2],f[p][s1][s2]);
                        inc(f[p^1][s1][s2],f[p][s1][s2]);
                    }int s2=0;
                    if((a[x]&s1)==0) inc(f[p^1][s1][s2|a[x]],f[p][s1][s2]);
                    if((a[x]&s2)==0) inc(f[p^1][s1|a[x]][s2],f[p][s1][s2]);
                    inc(f[p^1][s1][s2],f[p][s1][s2]);
                }p^=1;
            }for(int s1=0;s1<bin[m];++s1){
                    int st=bin[m]-1-s1;
                    for(int s2=st;s2;s2=(s2-1)&st){
                        inc(f[p][s1][s2],f[p][s1|bin[m]][s2]),f[p][s1|bin[m]][s2]=0;
                        inc(f[p][s1][s2],f[p][s1][s2|bin[m]]),f[p][s1][s2|bin[m]]=0;
                    }int s2=0;
                    inc(f[p][s1][s2],f[p][s1|bin[m]][s2]),f[p][s1|bin[m]][s2]=0;
                    inc(f[p][s1][s2],f[p][s1][s2|bin[m]]),f[p][s1][s2|bin[m]]=0;
            }
        }
        for(int j=1;j<=8;++j){  
            for(int x=prime[j];x<=n;x+=prime[j]){
                if(used[x]) continue;used[x]=1;memset(f[p^1],0,sizeof(f[p^1]));
                for(int s1=0;s1<bin[m];++s1){
                    int st=bin[m]-1-s1;
                    for(int s2=st;s2;s2=(s2-1)&st){
                        if((a[x]&s1)==0) inc(f[p^1][s1][s2|a[x]],f[p][s1][s2]);
                        if((a[x]&s2)==0) inc(f[p^1][s1|a[x]][s2],f[p][s1][s2]);
                        inc(f[p^1][s1][s2],f[p][s1][s2]);
                    }int s2=0;
                    if((a[x]&s1)==0) inc(f[p^1][s1][s2|a[x]],f[p][s1][s2]);
                    if((a[x]&s2)==0) inc(f[p^1][s1|a[x]][s2],f[p][s1][s2]);
                    inc(f[p^1][s1][s2],f[p][s1][s2]);
                }p^=1;
            }
        }ll ans=0;
        for(int s1=0;s1<bin[m];++s1){
            int st=bin[m]-1-s1;
            for(int s2=st;s2;s2=(s2-1)&st) inc(ans,f[p][s1][s2]);inc(ans,f[p][s1][0]);
        }printf("%lld\n",ans);
    }
}
int main(){
//  freopen("a.in","r",stdin);
    n=read();mod=read();getprime();
    if(n<=50) sol1::gao();
    else if(n<=100) sol2::gao();
    else sol3::gao();
    return 0;
}
Icefox_zhx
关注
专栏目录
【BZOJ 3067】Hyperdrome 哈希+状压
pbihao的博客
02-05 476
首先假设内存足够大的话可以用f[i]表示i这个状态出现了多少次，状态是指异或前缀和，那么根据题意，当且只有在i==0或者二进制表示下只有一个1的情况下是合法状态，所以每加一个字符，就更新当前前缀和，又由于A^B==C->A^C==B 那么可以预处理出所有合法状态，然后分别和当前状态异或后得到需要的状态，加入答案。 但是内存不够，不能用数组，amp貌似是一个不错的选择，呵呵，又T又mle，那么h
bzoj3886[Usaco2015 Jan]Moovie Mooving 状压DP
老臣
07-20 408
其实是很简单的状压..但是我太菜了没有想到。 一开始的DP方程就列错了，直接列了个求最终答案的。。 设f[i]表示状态为i时能看到的最远时间，转移显然，二分一下起始点就好了= =#include<cstdio> #include<algorithm> #include<cstring> #define fo(i,a,b) for(int i=a;i<=b;i++) #define fd(i,a,
【dp】【NOI 2015】【bzoj 4197】寿司晚宴
morestep的专栏
08-13 2912
4197: [Noi2015]寿司晚宴Time Limit: 10 Sec Memory Limit: 512 MB Submit: 58 Solved: 46 Description为了庆祝 NOI 的成功开幕，主办方为大家准备了一场寿司晚宴。小 G 和小 W 作为参加 NOI 的选手，也被邀请参加了寿司晚宴。在晚宴上,主办方为大家提供了 n−1 种不同的寿司，编号 1,2,3,…,n−1，其
bzoj4197
weixin_34124577的博客
07-10 82
这题现场想的思路方向都是对的，但限于现场和实力因素没能A 很显然我们会想到质因数的选取 如果某个质数p被W选了，那G就不能选含有质因子p的数 因此我们不难想到状压质数的选取情况，令f[i][j]为w质数选取状态为i，g质数选取状态j的方案数 但是500以内质数太多了怎么办？我们考虑大质数能不能分开考虑 考虑到sqrt(500)以外的质数，他们在每个数中最多出现一次，因此我们只要按这些质数...
[bzoj4197][Noi2015]寿司晚宴【dp】
D_Vanisher的博客
01-23 258
【题目描述】 Description 为了庆祝 NOI 的成功开幕，主办方为大家准备了一场寿司晚宴。小 G 和小 W 作为参加 NOI 的选手，也被邀请参加了寿司晚宴。 在晚宴上,主办方为大家提供了 n−1 种不同的寿司，编号 1,2,3,…,n−1，其中第 i 种寿司的美味度为 i+1 （即寿司的美味度为从 2 到 n）。 现在小 G 和小 W 希望每人选一些寿司种类来品尝，他们
BZOJ 4197 寿司晚宴 （状压dp）
w4149
10-18 311
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 （
BZOJ 4197 Noi2015 寿司晚宴
Fuxey
05-29 835
分层图背包
BZOJ 4169 LＭC的游戏 （博弈+树dp）
Mogician_Evian的博客
03-08 553
4169 LMC的游戏 【题目描述】 RHL 有一天看到 lmc 在玩一个游戏。 “愚蠢的人类哟,what are you doing”,RHL 说。 “我在玩一个游戏。现在这里有一个有 n 个结点的有根树,其中有 m 个叶子结点。这 m个叶子从 1 到 m 分别被给予了一个号码,每个叶子的号码都是独一无二的。一开始根节点有一个棋子,两个玩家每次行动将棋子移动到当前节点的一个...
bzoj1231-混乱的奶牛(状压dp)
Chase的博客
05-28 277
题意：bzoj1231 分析： 给定一段序列，让你随便排序，问有多少种排序使得相邻元素之差大于mmm。 显然是个状压dpdpdp，状态转移方程： dp[i∣(1&lt;&lt;k)][k]+=dp[i][j]dp[i|(1&lt;&lt;k)][k]+=dp[i][j]dp[i∣(1<<k)][k]+=dp[i][j]， 其中k∉i,j∈i,abs(a[...
【状压DP】状态压缩动态规划入门超详解
繁凡さん的博客
03-15 9870
状压DP一、概述1.状态压缩2.使用条件3.状压DP二、位运算三、例题引入入门例题【例1】填满棋盘入门例题【例二】玉米地 感觉好多讲状压DP的博客都有点乱，我就结合各路大佬的博客，加上我自己的理解，总结出一篇博客来，供初学者参考 一、概述 1.状态压缩 状态压缩就是使用某种方法，简明扼要地以最小代价来表示某种状态，通常是用一串01数字（二进制数）来表示各个点的状态。这就要求使用状态压缩的对象的点的...
【BZOJ4197】[Noi2015]寿司晚宴【状压DP】【背包】
BraketBN
04-08 665
【题目链接】 一开始取了前10个质数状压，然后就TLE了...改成前8个就没事了。 /* Pigonometry */ #include #include #include using namespace std; typedef long long LL; const int maxn = 505, maxk = 505; int n, p, prime[maxn], cn
【BZOJ4197】寿司晚宴，状态压缩DP
xiaoyimi的博客
06-29 469
term final
BZOJ4197 [Noi2015]寿司晚宴
neither_nor
08-19 564
注意到对于任意一个数，他与别的数是否互质取决于他俩是否有相同的质因数，而对于任意的x，x的大于根号x的质因数最多有一个，而n 所以我们可以考虑对每个人取的小于根号500的质数的状态状压，然后枚举每个数，进行DP f[i][j]表示第一个人小于根号五百的质数的状态为i，第二个人为j的方案数 枚举当前数，设当前数为x，x去掉所有小于根号五百的质因子之后的数为y 对于所有y相等但不为1的数（每个
bzoj4197 [NOI2015] [状压dp] 寿司晚宴
某炒鸡无底蒟蒻的博客
10-30 325
Description 为了庆祝 NOI 的成功开幕，主办方为大家准备了一场寿司晚宴。小 G 和小 W 作为参加 NOI 的选手，也被邀请参加了寿司晚宴。 在晚宴上,主办方为大家提供了 n−1 种不同的寿司，编号 1,2,3,…,n−1，其中第 i 种寿司的美味度为 i+1 （即寿司的美味度为从 2 到 n）。 现在小 G 和小 W 希望每人选一些寿司种类来品尝，他们规定一种
【bzoj4197】[Noi2015]寿司晚宴 状压DP
DQS的树洞
03-13 1618
Description为了庆祝 NOI 的成功开幕，主办方为大家准备了一场寿司晚宴。小 G 和小 W 作为参加 NOI 的选手，也被邀请参加了寿司晚宴。在晚宴上,主办方为大家提供了 n−1 种不同的寿司，编号 1,2,3,…,n−1，其中第 i 种寿司的美味度为 i+1 （即寿司的美味度为从 2 到 n）。 现在小 G 和小 W 希望每人选一些寿司种类来品尝，他们规定一种品尝方案为不和谐的当且仅当
BZOJ4197: [Noi2015]寿司晚宴 状压DP
リーインカーネイション
04-11 481
题意：把2~n这些数分成两个可空集合，要求两个集合中任取一对数都必须互质，求方案数 2≤n≤500 比较显然的是可以把每个数分解质因数，相当于两个集合中不能存在同一质因数。 但是质因数个数可能很多，并不好设状态。 注意到每个数中最多只有一个大于等于根号n的质因子，而小于根号n的质数只有8个，那么将所有数按照大质因子分类，规定每一类最多只能被一个集合选，就可以把前八个质数状压来表示状态了。
【bzoj4197】[Noi2015]寿司晚宴
wang_sj的博客
05-24 473
Description为了庆祝 NOI 的成功开幕，主办方为大家准备了一场寿司晚宴。小 G 和小 W 作为参加 NOI 的选手，也被邀请参加了寿司晚宴。在晚宴上,主办方为大家提供了 n−1 种不同的寿司，编号 1,2,3,…,n−1，其中第 i 种寿司的美味度为 i+1 （即寿司的美味度为从 2 到 n）。 现在小 G 和小 W 希望每人选一些寿司种类来品尝，他们规定一种品尝方案为不和谐的当且仅当
【bzoj4197】[Noi2015]寿司晚宴 dp
Oxer的专栏
05-10 1566
因为每个数只有一个大于根号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][
NOIP模拟9.23（TYVJ NOIP2017模拟赛D1）
Icefox的博客
09-27 708
T1 回形遍历 T2 排列 T3 近似排列计数
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值