先貼一個題解吧，最近懶得要死2333，可能是太弱的原因吧，總是扒題解，（甚至連題解都看不懂了），blog也沒更新，GG

http://blog.csdn.net/werkeytom_ftd/article/details/52527740

容斥原理真的很神奇233

1 #include <bits/stdc++.h> 2 #define LL long long 3 using namespace std; 4 5 const int maxn=500; 6 const int mod=1e9+7; 7 8 int fac[maxn],inv[maxn]; 9 void pre() 10 { 11 fac[0]=1; for (int i=1; i<=400; i++) fac[i]=(LL)fac[i-1]*i%mod; 12 inv[0]=inv[1]=1; 13 for (int i=2; i<=400; i++) inv[i]=(LL)(mod-mod/i)*inv[mod%i]%mod; 14 for (int i=2; i<=400; i++) inv[i]=(LL)inv[i]*inv[i-1]%mod; 15 } 16 int ksm(int x, int p) 17 { 18 int sum=1; 19 for (;p;p>>=1,x=(LL)x*x%mod) 20 if (p&1) sum=(LL)sum*x%mod; 21 return sum; 22 } 23 int C(int n, int m) 24 { 25 return (LL)fac[n]*inv[m]%mod*inv[n-m]%mod; 26 } 27 28 int n,m,p,ans; 29 int main() 30 { 31 cin>>n>>m>>p; pre(); 32 for (int i=0; i<=n; i++) 33 for (int k=0; k<=p; k++) 34 { 35 int qwq=(LL)C(n,i)*C(p,k)%mod; 36 int orz=ksm((1-ksm(k+1,i)+mod)%mod,m); 37 qwq=(LL)qwq*orz%mod; 38 if ((n+m+p-i-k)&1) qwq=-qwq; 39 ans=(ans+qwq)%mod; 40 } 41 printf("%d\n",(ans+mod)%mod); 42 return 0; 43 }

?

轉載于:https://www.cnblogs.com/ccd2333/p/6792629.html

總結

以上是生活随笔為你收集整理的bzoj 4487: [Jsoi2015]染色问题的全部內容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯，歡迎將生活随笔推薦給好友。