P2396 yyy loves Maths VII

题目背景

yyy对某些数字有着情有独钟的喜爱,他叫他们为幸运数字;然而他作死太多,所以把自己讨厌的数字成为"厄运数字"

题目描述

一群同学在和yyy玩一个游戏

每次,他们会给yyy n张卡片,卡片上有数字,所有的数字都是"幸运数字",我们认为第i张卡片上数字是ai

每次yyy可以选择向前走ai步并且丢掉第i张卡片

当他手上没有卡片的时候他就赢了

但是呢,大家对"厄运数字"的位置布置下了陷阱,如果yyy停在这个格子上,那么他就输了

(注意:即使到了终点,但是这个位置是厄运数字,那么也输了)

现在,有些同学开始问:

yyy有多大的概率会赢呢?

大家觉得这是个好问题

有人立即让yyy写个程序

"电脑运行速度很快!24的阶乘也不过就620448401733239439360000,yyy你快写个程序来算一算"

yyy表示很无语,他表示他不想算概率,最多算算赢的方案数,而且是%1,000,000,007以后的值

大家都不会写程序,只好妥协

但是这时候yyy为难了,24!太大了,要跑好长时间.

他时间严重不够!需要你的帮助!

由于yyy人格分裂，某个数字可能既属于幸运数字又属于厄运数字。

输入输出格式

输入格式：

第一行n

下面一行n张卡片

第三行m 表示yyy的厄运数字个数(最多2个)

最后一行是m个厄运数字

输出格式：

方案数%1,000,000,007

输入输出样例

输入样例#1：

8

1 3 1 5 2 2 2 3

0

输出样例#1：

输入样例#2：

24

1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1

2

10 15

输出样例#2：

说明

数据范围:

10%的数据n<=10

50%的数据n<=23

100%的数据n<=24

sol：这样的数据范围，感觉不是爆搜，状压dp无疑

dp[z]表示状态为z时没有厄运数字的方案数

转移较易：枚举一个在集合z中的数字i，dp[z]+=dp[z^i]

注意判断Dis，即一个集合z的距离和Dis[z]为厄运数字，那么不能进行转移

要用lowbit进行帮助转移

lowbit(x)表示一个数在二进制意义下第一位非0的数位可以帮助枚举一个集合，比（1~n）要快

register帮助卡常

#include <bits/stdc++.h>

using namespace std;

typedef int ll;

inline ll read()

{

    ll s=;

    bool f=;

    char ch=' ';

    while(!isdigit(ch))

    {

        f|=(ch=='-'); ch=getchar();

    }

    while(isdigit(ch))

    {

        s=(s<<)+(s<<)+(ch^); ch=getchar();

    }

    return (f)?(-s):(s);

}

#define R(x) x=read()

inline void write(ll x)

{

    if(x<)

    {

        putchar('-'); x=-x;

    }

    if(x<)

    {

        putchar(x+''); return;

    }

    write(x/);

    putchar((x%)+'');

    return;

}

#define W(x) write(x),putchar(' ')

#define Wl(x) write(x),putchar('\n')

const int N=(<<)+,Mod=;

int n,m,B1,B2,Dis[N],dp[N];

#define lowbit(x) ((x)&(-x))

int main()

{

    register int i,j;

    R(n);

    for(i=;i<=n;i++) R(Dis[<<(i-)]);

    R(m);

    if(m>) R(B1); if(m>) R(B2);

    dp[]=;

    for(i=;i<(<<n);i++)

    {

        Dis[i]=Dis[i^(lowbit(i))]+Dis[lowbit(i)];

        if(Dis[i]==B1||Dis[i]==B2) continue;

        for(j=i;j;j^=lowbit(j))

        {

            dp[i]+=dp[i^lowbit(j)];

            dp[i]-=(dp[i]>=Mod)?Mod:;

        }

    }

    Wl(dp[(<<n)-]);

    return ;

}

/*

input

8

1 3 1 5 2 2 2 3

0

output

40320

input

24

1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1

2

10 15

output

0

*/

巴特西

洛谷P2396 yyy loves Maths VII

P2396 yyy loves Maths VII

题目背景

题目描述

输入输出格式

输入输出样例

说明

最新文章

热门文章