[LOJ#6160][美团 CodeM 初赛 Round A][容斥][递推]二分图染色

最新推荐文章于 2021-07-16 17:51:39 发布

原创最新推荐文章于 2021-07-16 17:51:39 发布 · 704 阅读

0 ·

CC 4.0 BY-SA版权

文章标签：

#容斥 #递推

容斥原理同时被 2 个专栏收录

22 篇文章

订阅专栏

二分图

6 篇文章

订阅专栏

本文提供了一道关于组合数学的算法题目解决方案，使用了动态规划、组合数计算等技巧，适用于算法竞赛及面试准备。

题解

#include <cstdio>
#include <iostream>
#include <algorithm>

using namespace std;

typedef long long ll;

const int N=1e7+10,P=1e9+7;

int n,f[N],fac[N],inv[N];

inline ll C(int x,int y){
  return 1LL*fac[x]*inv[y]%P*inv[x-y]%P;
}

inline ll A(int x,int y){
  return 1LL*fac[x]*inv[x-y]%P;
}

int main(){
  scanf("%d",&n);
  f[0]=1; f[1]=2;
  for(int i=2;i<=n;i++)
    f[i]=(2LL*i*f[i-1]%P+P-1LL*(i-1)*(i-1)%P*f[i-2]%P)%P;
  fac[0]=inv[0]=inv[1]=1;
  for(int i=1;i<=n;i++) fac[i]=1LL*fac[i-1]*i%P;
  for(int i=2;i<=n;i++) inv[i]=1LL*(P-P/i)*inv[P%i]%P;
  for(int i=1;i<=n;i++) inv[i]=1LL*inv[i-1]*inv[i]%P;
  int ans=1LL*f[n]*f[n]%P;
  for(int i=1,p=-1;i<=n;i++,p=-p)
    ans=(ans+(P+1LL*p*f[n-i]*f[n-i]%P*C(n,i)%P*A(n,i)%P)%P)%P;
  printf("%d\n",ans);
  return 0;
}