codeforces 895C Square Subsets 状压dp

原创 2018年04月17日 00:40:12

题目大意:给你一个数组,问有多少个非空子集满足集合中所有数字的乘积是一个平方数。结果mod 1e9+7。

首先就可以想到把所有数出现的次数记录一下。

由于唯一分解定理可知若一个数是平方数,那么它唯一分解后质因子的指数一定是偶数。

而且数字都是小于70的,小于70的质数只有20个,于是就可以进行状压dp了。

第 i 位为1表示唯一分解后第 i 个质数的质数为奇数,否则为偶数。

接着考虑状态转移。

若数 k 没出现,即 cnt[ k ] =0,那么 dp[ k ][ j ]=dp[ k-1 ][ j ],

若 k 出现,且出现次数为 n ,从 n 个中选出现 m 次,显然若 m 为偶数不影响状态,m为奇数则影响。

因为 C (0,n) + C ( 1 ,n) +....+C (n,n)= 2^n, 所以 C (0,n) + C(2,n)+...=C ( 1,n ) +C(3 , n )+... = 2 ^ (n-1) 。

那么状态转移方程为

dp [i][j] = dp[i][j]+ dp[i-1][j] * 2 ^( n-1);

dp [i][ j ^ status[k] ] = dp[i][ j^status[k] + dp[i-1][j] * 2^(n-1);

由于内存问题所以要注意使用滚动数组。

最后要减去全部选的 1 不要忘记。

#include <bits/stdc++.h>

using namespace std;

typedef long long ll;

const int maxn=1e5+5;
const ll mod=1e9+7;
ll p[100005];
ll status[71];
ll dp[2][(1<<19)];
ll cnt[71];
ll prime[71];
bool vis[71];

void init()
{
    p[0]=1;
    memset(status,0,sizeof(status));
    memset(vis,false,sizeof(vis));
    for(int i=1;i<maxn;i++)
        p[i]=(p[i-1]*(ll)2)%mod;
    vis[0]=vis[1]=true;
    for(int i=2;i<71;i++)
    {
        if(!vis[i])
        {
            for(int j=2;i*j<71;j++)
                vis[i*j]=true;
        }
    }
    int cnt=0;
    for(int i=2;i<71;i++)
    {
        if(!vis[i])
            prime[cnt++]=i;
    }
    for(int i=1;i<71;i++)
    {
        int t=i;
        for(int j=0;j<cnt;j++)
        {
            int num=0;
            while(t%prime[j]==0)
            {
                num++;
                t/=prime[j];
            }
            if(num%2)
                status[i]+=(1<<j);
        }
    }
}

int main()
{
    int n,x;
    init();
    while(cin>>n)
    {
        memset(cnt,0,sizeof(cnt));
        for(int i=0;i<n;i++)
        {
            cin>>x;
            cnt[x]++;
        }
        memset(dp,0,sizeof(dp));
        dp[0][0]=1;
        for(int i=1;i<=70;i++)
        {
            int now=i&1,last=!now;
            for(int j=0;j<(1<<19);j++)
            {
                if(cnt[i])
                {
                    dp[now][j]=(dp[now][j]+dp[last][j]*p[cnt[i]-1]%mod)%mod;
                    dp[now][j^status[i]]=(dp[now][j^status[i]]+dp[last][j]*p[cnt[i]-1]%mod)%mod;
                }
                else dp[now][j]=dp[last][j];
            }
            for(int j=0;j<(1<<19);j++)
                dp[last][j]=0;
        }
        cout<<(dp[0][0]-1+mod)%mod<<endl;
    }
    return 0;
}



Codeforces895C. Square Subsets

C. Square Subsets time limit per test 4 seconds memory limit per test 256 megabytes ...
  • Jelly_acm
  • Jelly_acm
  • 2017-11-29 10:54:13
  • 495

895C - Square Subsets 状压DP

div2从没这么接近过rk100以内的水平。。 想到了奇偶性质模二的素因子表,想到了dp第一维逐渐添加数字,第二维枚举1 奈何第一维默认了是1-n,算下来复杂度1e10从而卡手做不下去 没完整做...
  • qq_34927456
  • qq_34927456
  • 2017-12-05 02:21:55
  • 52

Codeforces Round #448 (Div. 2) C. Square Subsets

C. Square Subsets time limit per test4 seconds memory limit per test256 megabytes inputstandard i...
  • zstu_zy
  • zstu_zy
  • 2017-11-27 12:14:11
  • 461

codeforces 895C Square Subsets

codeforces 895C Square Subsets
  • Rising_shit
  • Rising_shit
  • 2018-02-27 15:56:27
  • 40

codeforces 107C Arrangement (状压dp)

题意: 给出n个人和n个座位,给出m对限制,每对限制要求ai这个人的座位要排在bi这个人的前面。现在要求满足条件的第y-2001大的字典序的座位排列。 题解: 题目挺好的,处理方法很特别。对于这...
  • My_ACM_Dream
  • My_ACM_Dream
  • 2015-02-20 15:26:41
  • 618

895C - Square Subsets 状压DP + 离散化

C. Square Subsets time limit per test 4 seconds memory limit per test 256 megabytes ...
  • chudongfang2015
  • chudongfang2015
  • 2017-11-29 11:27:11
  • 124

Codeforces 16E Fish (状压dp+概率)

题意有n条鱼,他们相遇时会吃掉对方,给出他们相遇时双方获胜的概率,求这n条鱼最后剩下自己的概率。思路看范围就要考虑状压DP,dp[s]表示当前剩下的鱼的状态为s时的概率。 那么P(i吃掉j) = P...
  • Dreamon3
  • Dreamon3
  • 2016-05-08 19:00:05
  • 576

Codeforces 660E Different Subsets For All Tuples【组合数学】

看了官方题解+q神的讲解才懂。。。 智商问题。。 讲道理。。数学真的比脱单难啊。。。题目链接:http://codeforces.com/problemset/problem/660/E题意:给定...
  • Yukizzz
  • Yukizzz
  • 2016-04-12 12:33:26
  • 1502

CodeForces 895 C.Square Subsets(状压DP+线性基)

Description 给出nnn个数a1,...,ana1,...,ana_1,...,a_n,问从这nnn个数中选一非空子集使得乘积是完全平方数的方案数 Input 第一行一整数nnn表示数...
  • V5ZSQ
  • V5ZSQ
  • 2018-01-15 15:09:29
  • 114

CodeForces 449 D.Jzzhu and Numbers(状压DP+容斥原理)

Description 给出一个n个数a[i]的集合S,问S有多少个子集&起来是0 Input 第一行一整数n表示集合元素数量,之后n个整数a[i]表示该集合中的元素(1...
  • V5ZSQ
  • V5ZSQ
  • 2017-04-12 16:20:16
  • 312
收藏助手
不良信息举报
您举报文章:codeforces 895C Square Subsets 状压dp
举报原因:
原因补充:

(最多只允许输入30个字)