bzoj 4487: [Jsoi2015]染色问题

先贴一个题解吧，最近懒得要死2333，可能是太弱的原因吧，总是扒题解，（甚至连题解都看不懂了），blog也没更新，GG

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

容斥原理真的很神奇233

 #include <bits/stdc++.h>
 #define LL long long
 using namespace std;

 const int maxn=;
 const int mod=1e9+;

 int fac[maxn],inv[maxn];
 void pre()
 {
     fac[]=; for (int i=; i<=; i++) fac[i]=(LL)fac[i-]*i%mod;
     inv[]=inv[]=;
     for (int i=; i<=; i++) inv[i]=(LL)(mod-mod/i)*inv[mod%i]%mod;
     for (int i=; i<=; i++) inv[i]=(LL)inv[i]*inv[i-]%mod;
 }
 int ksm(int x, int p)
 {
     int sum=;
     for (;p;p>>=,x=(LL)x*x%mod)
         if (p&) sum=(LL)sum*x%mod;
     return sum;
 }
 int C(int n, int m)
 {
     return (LL)fac[n]*inv[m]%mod*inv[n-m]%mod;
 }

 int n,m,p,ans;
 int main()
 {
     cin>>n>>m>>p; pre();
     for (int i=; i<=n; i++)
         for (int k=; k<=p; k++)
         {
             int qwq=(LL)C(n,i)*C(p,k)%mod;
             int orz=ksm((-ksm(k+,i)+mod)%mod,m);
             qwq=(LL)qwq*orz%mod;
             if ((n+m+p-i-k)&) qwq=-qwq;
             ans=(ans+qwq)%mod;
         }
     printf("%d\n",(ans+mod)%mod);
     return ;
 }