bzoj2863:愤怒的元首

传送门

考虑到这样一个性质,一个入度为0的点连一条边到一个DAG中,依然是一个DAG

于是设\(f(i)\)\(i\)个点组成的DAG方案数,

那么\(n\)个节点的DAG中至少有\(i\)个节点入度为\(0\)方案数为\(f(n-i)\binom{n}{i}2^{i(n-i)}\)

但是入度为\(0\)的点数为\(0\)时,方案数为\(0\),也就是
\[ \sum_{i=0}^{n}(-1)^if(n-i)\binom{n}{i}2^{i(n-i)}=0\\ f(n)=\sum_{i=1}^{n}(-1)^{i-1}f(n-i)\binom{n}{i}2^{i(n-i)} \]
代码:

#include<cstdio>
#include<iostream>
#include<algorithm>
using namespace std;
void read(int &x) {
    char ch; bool ok;
    for(ok=0,ch=getchar(); !isdigit(ch); ch=getchar()) if(ch=='-') ok=1;
    for(x=0; isdigit(ch); x=x*10+ch-'0',ch=getchar()); if(ok) x=-x;
}
#define rg register
const int maxn=3e3+10,mod=1e9+7;
int n,f[maxn],fac[maxn],inv[maxn];
int mul(int x,int y){return 1ll*x*y-1ll*x*y/mod*mod;}
int add(int x,int y){return x+y>=mod?x+y-mod:(x+y<0?x+y+mod:x+y);}
int mi(int a,int b)
{
    int ans=1;
    while(b)
    {
        if(b&1)ans=mul(a,ans);
        b>>=1,a=mul(a,a);
    }
    return ans;
}
int c(int x,int y){return mul(fac[x],mul(inv[y],inv[x-y]));}
int main()
{
    read(n),f[0]=f[1]=fac[0]=inv[0]=1;
    for(rg int i=1;i<=n;i++)fac[i]=mul(i,fac[i-1]);
    inv[n]=mi(fac[n],mod-2);
    for(rg int i=n-1;i;i--)inv[i]=mul(inv[i+1],i+1);
    for(rg int i=2;i<=n;i++)
        for(rg int j=0;j<i;j++)
            f[i]=add(f[i],mul(mul(mul(((i-j-1)&1?-1:1),f[j]),c(i,i-j)),mi(2,(i-j)*j)));
    printf("%d\n",f[n]);
}

转载于:https://www.cnblogs.com/lcxer/p/10691784.html

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值