【XSY1519】彩灯节 DP 数学 第二类斯特林数

最新推荐文章于 2021-09-28 10:30:51 发布
最新推荐文章于 2021-09-28 10:30:51 发布
阅读量492 收藏 2
点赞数
分类专栏： DP 数学&数论 数学&数论--第二类斯特林数
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/ez_yww/article/details/77368524
版权

题目大意

​　　有 n 盏灯，m个限制。每个限制 (x,y) 表示第 x 盏灯与第y盏灯之间必须且只能亮一盏。

​　　记一种情况 x 亮着的灯的数量为fx，求 ∑(fx)k

​　　 n≤200000,k≤100

题解

​　　我们先把整张图黑白染色。

​　　如果不是二分图就无解。

​　　我们发现两个不同的联通分量的灯的状态是没有关系的。

​　　我们可以考虑DP：

​　　 fi,j= 前 i 个联通分量中亮的彩灯个数的j次方和

​　　 ai,j= 第 i 个联通分量中亮的彩灯个数的j次方和

​　　根据二项式定理 (a+b)n=∑ni=0(ni)aibn−i 有：

fi,j=∑k=0j(jk)ai−1,k×fi−1,j−k

=∑k=0jj!k!(j−k)!ai−1,k×fi−1,j−k

=j!∑k=0jai−1,kk!×fi−1,j−k(j−k)!

​　　观察到 p=1004535809 是一个NTT模数（原根为 3 ），所以可以用NTT加速。

​　　没了

​　　时间复杂度：O(nklogk)

题解2

​　　FFT很明显会TLE

​　　有一个公式：

mn=∑k=0m(mk)S(n,k)k!

​　　后面两个东西很容易算，只用考虑第一个怎么求

​　　这不就是个组合数吗

​　　直接DP就行了。

​　　时间复杂度： O(nk)

代码

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cstdlib>
#include<ctime>
#include<utility>
using namespace std;
typedef long long ll;
typedef pair<int,int> pii;
ll p=1004535809;
ll fp(ll a,ll b)
{
    ll s=1;
    while(b)
    {
        if(b&1)
            s=s*a%p;
        a=a*a%p;
        b>>=1;
    }
    return s;
}
namespace ntt
{
    int N;
    ll w1[500010];
    ll w2[500010];
    int rev[500010];
    void get(int n)
    {
        N=1;
        while(N<n)
            N<<=1;
        int i;
        for(i=2;i<=N;i++)
        {
            w1[i]=fp(3,(p-1)/i);
            w2[i]=fp(w1[i],p-2);
        }
        rev[0]=0;
        for(i=1;i<N;i++)
            rev[i]=(rev[i>>1]>>1)|(i&1?N>>1:0);
    }
    void ntt(ll *a,int t)
    {
        int i,j,k;
        ll u,v,w,wn;
        for(i=0;i<N;i++)
            if(rev[i]<i)
                swap(a[i],a[rev[i]]);
        for(i=2;i<=N;i<<=1)
        {
            wn=t?w1[i]:w2[i];
            for(j=0;j<N;j+=i)
            {
                w=1;
                for(k=j;k<j+i/2;k++)
                {
                    u=a[k];
                    v=a[k+i/2]*w%p;
                    a[k]=(u+v)%p;
                    a[k+i/2]=(u-v+p)%p;
                    w=w*wn%p;
                }
            }
        }
        if(!t)
        {
            ll inv=fp(N,p-2);
            for(i=0;i<N;i++)
                a[i]=a[i]*inv%p;
        }
    }
}
struct list
{
    int v[500010];
    int t[500010];
    int h[200010];
    int n;
    list()
    {
        n=0;
        memset(h,0,sizeof h);
    }
    void add(int x,int y)
    {
        n++;
        v[n]=y;
        t[n]=h[x];
        h[x]=n;
    }
};
list l;
int vis[200010];
void failed()
{
    putchar('0');
    exit(0);
}
int s1,s2;
void dfs(int x,int fa,int v)
{
    vis[x]=v;
    if(v)
        s1++;
    else
        s2++;
    int i;
    for(i=l.h[x];i;i=l.t[i])
        if(l.v[i]!=fa)
        {
            if(vis[l.v[i]]==-1)
                dfs(l.v[i],x,v^1);
            else if(vis[l.v[i]]==vis[x])
                failed();
        }
}
ll f[500010];
ll a[500010];
int n,m,k;
ll fac[500010];
int main()
{
    freopen("bulb.in","r",stdin);
    freopen("bulb.out","w",stdout);
    scanf("%d%d%d",&n,&m,&k);
    ntt::get(2*(k+1));
    int i,j,x,y;
    fac[0]=1;
    for(i=1;i<=1000;i++)
        fac[i]=fac[i-1]*i%p;
    for(i=1;i<=m;i++)
    {
        scanf("%d%d",&x,&y);
        l.add(x,y);
        l.add(y,x);
    }
    memset(vis,-1,sizeof vis);
    f[0]=1;
    for(i=1;i<=n;i++)
        if(vis[i]==-1)
        {
            s1=s2=0;
            dfs(i,0,1);
            for(j=0;j<=k;j++)
                a[j]=(fp(s1,j)+fp(s2,j))%p;
            for(j=k+1;j<ntt::N;j++)
                a[j]=f[j]=0;
            for(j=0;j<ntt::N;j++)
            {
                ll inv=fp(fac[j],p-2);
                a[j]=a[j]*inv%p;
                f[j]=f[j]*inv%p;
            }
            ntt::ntt(a,1);
            ntt::ntt(f,1);
            for(j=0;j<ntt::N;j++)
                f[j]=f[j]*a[j]%p;
            ntt::ntt(f,0);
            for(j=0;j<=k;j++)
                f[j]=f[j]*fac[j]%p;
            // do something here ...
        }
    printf("%lld\n",f[k]);
    return 0;
}

代码2

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cstdlib>
#include<ctime>
#include<utility>
using namespace std;
typedef long long ll;
typedef pair<int,int> pii;
ll p=1004535809;
ll fp(ll a,ll b)
{
    ll s=1;
    while(b)
    {
        if(b&1)
            s=s*a%p;
        a=a*a%p;
        b>>=1;
    }
    return s;
}
struct list
{
    int v[500010];
    int t[500010];
    int h[200010];
    int n;
    list()
    {
        n=0;
        memset(h,0,sizeof h);
    }
    void add(int x,int y)
    {
        n++;
        v[n]=y;
        t[n]=h[x];
        h[x]=n;
    }
};
list l;
int vis[200010];
void failed()
{
    putchar('0');
    exit(0);
}
int s1,s2;
void dfs(int x,int fa,int v)
{
    vis[x]=v;
    if(v)
        s1++;
    else
        s2++;
    int i;
    for(i=l.h[x];i;i=l.t[i])
        if(l.v[i]!=fa)
        {
            if(vis[l.v[i]]==-1)
                dfs(l.v[i],x,v^1);
            else if(vis[l.v[i]]==vis[x])
                failed();
        }
}
int n,m,k;
ll fac[500010];
ll f[510];
ll f2[510];
ll f1[510];
void dp1()
{
    int i;
    for(i=k;i>=1;i--)
        f1[i]=(f1[i]+f1[i-1])%p;
}
void dp2()
{
    int i;
    for(i=k;i>=1;i--)
        f2[i]=(f2[i]+f2[i-1])%p;
}
ll s[510][510];
int main()
{
//  freopen("bulb.in","r",stdin);
//  freopen("bulb4.out","w",stdout);
    scanf("%d%d%d",&n,&m,&k);
    int i,j,x,y;
    fac[0]=1;
    for(i=1;i<=1000;i++)
        fac[i]=fac[i-1]*i%p;
    s[0][0]=1;
    for(i=1;i<=100;i++)
    {
        s[i][0]=0;
        for(j=1;j<=100;j++)
            s[i][j]=(s[i-1][j-1]+s[i-1][j]*j%p)%p;
    }
    for(i=1;i<=m;i++)
    {
        scanf("%d%d",&x,&y);
        l.add(x,y);
        l.add(y,x);
    }
    memset(vis,-1,sizeof vis);
    f[0]=1;
    for(i=1;i<=n;i++)
        if(vis[i]==-1)
        {
            s1=s2=0;
            dfs(i,0,1);
            for(j=0;j<=k;j++)
                f1[j]=f2[j]=f[j];
            for(j=1;j<=s1;j++)
                dp1();
            for(j=1;j<=s2;j++)
                dp2();
            for(j=0;j<=k;j++)
                f[j]=(f1[j]+f2[j])%p;
        }
    ll ans=0;
    for(i=0;i<=k;i++)
        ans=(ans+f[i]*s[k][i]%p*fac[i]%p)%p;
    printf("%lld\n",ans);
    return 0;
}
ez_yww
关注
专栏目录
【XSY2551】球（DP）（组合数学）
ezoi_ly的博客
09-10 420
Description Snuke喜欢五颜六色的球，他一共有 n∗k 个。有 n 种颜色，每种颜色有 k 个球。现在他将所有球以任意顺序排成一排，然后将每种颜色的最左边的球染成一种新的颜色。求重新染色后的球的颜色序列有多少种，答案对1000000007取模。 Input 一行两个整数 n,k 。 Output 一行一个整数表示方案数。 Sample Input sample input 1: 2 ...
【XSY2536】【BZOJ2655】calc DP 数学 拉格朗日插值
ez_yww的博客
08-16 1431
题目大意​　　一个序列a1,…,ana_1,\ldots,a_n是合法的，当且仅当： ​　　长度为给定的nn。 ​　　a1,…,ana_1,\ldots,a_n都是[1,m][1,m]中的整数。 ​　　a1,…,ana_1,\ldots,a_n互不相等。 ​　　一个序列的值定义为它里面所有数的乘积，即a1×a2×⋯×ana_1\times a_2\times\cdots\times a_n。
【BZOJ】2159: Crash 的文明世界-第二类斯特林数
远行客
12-03 209
巧妙的转化
BZOJ 5093: 图的价值 第二类斯特林数 NTT
BlackJack
02-26 1812
5093: 图的价值Time Limit: 30 Sec  Memory Limit: 256 MBSubmit: 176  Solved: 89[Submit][Status][Discuss]Description“简单无向图”是指无重边、无自环的无向图（不一定连通）。一个带标号的图的价值定义为每个点度数的k次方的和。给定n和k，请计算所有n个点的带标号的简单无向图的价值之和。因为答案很大，请...
[BZOJ]4555 [TJOI2016&HEOI2016] 求和 第二类斯特林数 + NTT
Le Petit Prince
01-19 661
4555: [Tjoi2016&Heoi2016]求和 Time Limit: 40 Sec  Memory Limit: 128 MB Submit: 525  Solved: 418 [Submit][Status][Discuss] Description 在2016年，佳媛姐姐刚刚学习了第二类斯特林数，非常开心。 现在他想计算这样一个函数的值: S(i, j)表示第二
【BZOJ5093】图的价值（第二类斯特林数，NTT）
დ♂Hany01's Blog Space♂ღ
03-24 267
Description “简单无向图”是指无重边、无自环的无向图（不一定连通）。 一个带标号的图的价值定义为每个点度数的k次方的和。 给定n和k，请计算所有n个点的带标号的简单无向图的价值之和。 因为答案很大，请对998244353取模输出。 Solution 容易得到题目就是要求： n×2(n−1)(n−2)2∑i=0n−1(n−1i)ikn×2(n−1)(n−2)2∑i=0n−1...
第二类斯特林数合集
辣鸡葫芦娃
04-08 974
留坑。。HDU2643求一个n人参加的比赛，最终的排名情况有多少种，注意可以并列。思路1：直接DP。dp[i][j]表示i个人，总共有j种名次。dp[i][j] =j*(dp[i-1][j-1]+dp[i-1][j]) 如果第i个人单独占有了一个名次，则i-1个人占有j-1个名次，否则，i-1个人占有j个名次，然后枚举第i个人在哪一个名次。答案=sigma{dp[n][i]}(1&lt;=i&lt...
[XSY] 计数（DP，NTT，分治）
Edward The Bunny 的博客
03-23 298
计数 考虑转化题目，变为网格上有若干个点，要从(0,0)(0,0)(0,0)走到(n,an+1)(n,a_{n+1})(n,an+1​) ，每一步只能往右走一步或往上走一步，且若当前在(i,j)(i,j)(i,j) ，必须满足0≤j≤ai+10\leq j\leq a_{i+1}0≤j≤ai+1​，其中an+1=ana_{n+1}=a_nan+1​=an​ 。要对每个kkk求出有多少条恰走过kkk次形如(i,ai+1)→(i+1,ai+1)(i,a_{i+1})\to(i+1,a_{i+1})(i,ai+
[CSA35G][XSY3318]Counting Quests （DP）
Edward The Bunny 的博客
09-28 214
XSY3318 CSA35G 对于一个询问区间的集合SSS，求出每一个数被哪些区间覆盖了，记为SiS_iSi​。 要能保证猜出选中数，当且仅当每个数的SiS_iSi​互不相同。 考虑求出不满足要求的集合SSS的个数。 首先可以观察得到SiS_iSi​的一个性质：若a<b<c<da<b<c<da<b<c<d且Sa=Sc,Sb=SdS_a=S_c,S_b=S_dSa​=Sc​,Sb​=Sd​，则必有Sa=Sb=Sc=SdS_a=S_b=S_c=S_dSa​=
[XSY] PQ树（区间DP）
Edward The Bunny 的博客
03-09 451
PQ树 对于非叶子节点，它只能调换它的儿子的顺序，却不能调换整颗子树的顺序 合法的PQ树都能表示出1,2,3,…,n 这个排列 有以上两点可以得出：PQ树的每棵子树对应的一定都是一段连续的数字区间 考虑区间DP，设f(l,r)f(l,r)f(l,r)表示数字l~r构成一棵PQ树的方案数 数字l~r 能够构成一棵子树 当且仅当 数字l~r 在题目给出的所有排列中都在一段连续区间中，我们定义这样的 [l,r] 是“好的” 若[l,r]不是"好的"，f(l,r)=0f(l,r)=0f(l,r)=0 若[l,r]
HDOJ 题目2643 Rank(斯特林第二类数）
Jogging_Clown
11-04 778
Rank Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Others) Total Submission(s): 324    Accepted Submission(s): 175 Problem Description Recently in Teddy's hom
[第二类斯特林数 组合计数] 省选模拟赛 2 B. 两弹一星 missile
雯舞
03-18 1074
题目大意 一张无向图的权值定义为 xkx^k，其中 xx 是图中结构为树的连通块个数。给定 nn, kk, 求出所有 nn 个点带标号的的简单无向图的权值和，对 998244353998244353 取模。 令 xix_i 表示树联通块 i 是否存在 , 图的权值为 (∑xi)m(\sum x_i)^m。对于某 kk 个联通块，如果同时出现，那么贡献为 S(m,k)∗k!S(m, k) * k!
bzoj5093 图的价值【第二类斯特林数+NTT】
cdsszjj的博客
01-16 831
解题思路： 考虑直接枚举每个点连的边，图中其余点之间随便连，那么直接推式子： ans=n∗2C2n−1∑i=0n−1Cin−1ikans=n*2^{C_{n-1}^2}\sum\limits_{i=0}^{n-1}C_{n-1}^ii^k 不妨用nn表示n−1n-1，而iki^k可用第二类斯特林数表示，则有： ans=(n+1)∗2C2n∑i=0nCin∑j=0iSjkCjij!ans=(
[省选前题目整理][HDU 2512]一卡通大冒险(第二类斯特林数+贝尔数)
某沙茶的OI代码库
04-10 1144
题目链接http://acm.hdu.edu.cn/showproblem.php?pid=2512思路第一类斯特林数s2[i][j]=s_2[i][j]=将jj个互不相同的物品划分成jj个非空集合的方案数。 s2[i][j]=s2[i−1][j−1]+js2[i−1][j]s_2[i][j]=s_2[i-1][j-1]+js_2[i-1][j] 递推式的解释：对于第ii个物品有两种情况：
SDNU 1011（第二类stirling数）
Charles' blog
10-24 578
1011.盒子与球 Time Limit: 1000 MS    Memory Limit: 32768 KB Total Submission(s): 667    Accepted Submission(s): 165 Description 现有r个互不相同的盒子和n个互不相同的球，要将这n个球放入r个盒子中，且不允许有空盒子。则有多少种放
AtCoder Regular Contest 096 E - Everything on It 容斥+第二类斯特林数
wohenshuai
04-23 293
题意 给你一个n，一共有2^n种组合，你要选择一些组合，使得每个数都出现至少两次，答案模m，m是一个质数，n&lt;=5000 分析 这道题是计数题，很容易想到容斥 ∑i=0n(−1)i22n−i(ni)∑j=0i(ij)∑k=1jS(j,k)(2n−i)k∑i=0n(−1)i22n−i(ni)∑j=0i(ij)∑k=1jS(j,k)(2n−i)k\sum\limits_{i=0}^{n}...
HDU 6143 Killer Names （第二类斯特林数or容斥）
K键盘里的青春K
08-17 793
Killer Names Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 65536/65536 K (Java/Others) Total Submission(s): 235    Accepted Submission(s): 128 Problem Description > Galen Marek,
HDU 6143 Killer Names [第二类斯特林数]
ACTerminate的博客
08-17 516
HDU 6143 Killer Names
hdu 4746 Bell 中国剩余定理+矩阵乘法+第二类斯特林数 （2013网络赛）
Nicolas
10-02 1800
google   “Bell number”  能得出以下结论：(注意p是质数) 题目中给的mod不是质数怎么办？ 把它拆成多个质数(31, 37, 41, 43, 47)分别计算，然后用中国剩余定理合并一下。 n很大,我们先预处理出0--p的Bell数，然后矩阵乘法优化即可。 #include #include #include using namespace std
http://localhost/api/xsy-auth-security/loginUser/login 403 (Forbidden)
最新发布
07-28
- *2* [c/c++使用libcurl库做http客户端及封装(HTTP_GET和HTTP_POST)](https://blog.csdn.net/xsy29000/article/details/103181267)[target="_blank" data-report-click={"spm":"1018.2226.3001.9630","extra":{"utm...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值