[NOI Online #3 提高组]优秀子序列欧拉函数,SOS DP

This way

题意:

在这里插入图片描述

题解:

首先要了解的一个概念就是:枚举 [ 1 , 2 x ] [1,2^x] [1,2x]的所有状态的子状态的时间复杂度是 3 x 3^x 3x左右。那么我们这道题肯定不能for一遍n然后枚举所有当前点的补集。然后 发现a的数据范围是2e5,那么就可以先记录所有状态有多少种,然后再枚举所有状态的同时枚举补集进行转移。最后0要特别的运算一下因为它每一个取不取都不会对结果造成影响。

#include<bits/stdc++.h>
using namespace std;
const int maxn=(1<<18)+5;
bitset<maxn>notprime;
int phi[maxn],prime[maxn],cnt=0;
void pre(){
	phi[1]=1;
	for(int i=2;i<=maxn;i++){
        if(!notprime[i]){
            prime[++cnt]=i;
            phi[i]=i-1;
        }
        for(int j=1;j<=cnt&&prime[j]*i<=maxn;j++){
            notprime[i*prime[j]]=1;
            if(i%prime[j]==0){
                phi[i*prime[j]]=phi[i]*prime[j];
                break;
            }
            else phi[i*prime[j]]=phi[i]*(prime[j]-1);
        }
    }
}
const int N=(1<<18)+5;
int a;
#define ll long long
ll dp[N],num[N];
const ll mod=1e9+7;
int main(){
	pre();
	int n;
	scanf("%d",&n);
	for(int i=1;i<=n;i++)
        scanf("%d",&a),num[a]++;
    dp[0]=1;
    for(int i=1;i<(1<<18);i++){
        for(int j=i;;j=(j-1)&i){
            int res=i^j;
            if(res>j)break;
            dp[i]=(dp[i]+dp[res]*num[j])%mod;
        }
    }
    ll ans=1;
    for(int i=1;i<N-1;i++)
        ans=(ans+dp[i]*phi[i+1])%mod;
    for(int i=1;i<=num[0];i++)ans=ans*2%mod;
    printf("%lld\n",ans);
    return 0;
}
  • 0
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 1
    评论
评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值