题目链接

题意,有m种颜色,给2n个位置染色,使左边n个和右边n个没有共同的颜色。

可以先递推求出恰用i种颜色染n个位置的方案数,然后枚举两边的染色数就可以了,代码很简单。

#include<bits/stdc++.h>
using namespace std;
typedef long long LL; //求C(a,b) ,其中,a,b大小大概不超过1e6级别,很多组数据
const LL mod=1e9+;
const LL M=2e3+;
LL fac[M+]; //阶乘
LL inv_of_fac[M+]; //阶乘的逆元
LL qpow(LL x,LL n) //求x^n%mod
{
LL ret=;
for(; n; n>>=)
{
if(n&) ret=ret*x%mod;
x=x*x%mod;
}
return ret;
}
void init()
{
fac[]=;
for(int i=; i<=M; i++)
fac[i]=fac[i-]*i%mod;
inv_of_fac[M]=qpow(fac[M],mod-);
for(int i=M-; i>=; i--)
inv_of_fac[i]=inv_of_fac[i+]*(i+)%mod;
}
LL C(LL a,LL b)
{
if(b<||a<b) return ;
return fac[a]*inv_of_fac[b]%mod*inv_of_fac[a-b]%mod;
}
//================================================================ int T,n,m;
LL cnt[]; //计数假如给定i种颜色,恰把i种颜色全部用上的方案数 void cal_cnt()
{
cnt[]=;
for(int i=;i<=m;i++)
{
cnt[i]=qpow(i,n);
for(int j=;j<i;j++)
cnt[i]=(cnt[i]-cnt[j]*C(i,j)%mod)%mod;
cnt[i]=(cnt[i]+mod)%mod;
}
} int main()
{
init();
scanf("%d",&T);
while(T--)
{
memset(cnt,,sizeof(cnt));
scanf("%d%d",&n,&m);
cal_cnt();
LL ans=;
for(int i=;i<=m;i++)
for(int j=;j<=m-i;j++)
ans=(ans+C(m,i)*cnt[i]%mod*C(m-i,j)%mod*cnt[j]%mod)%mod;
printf("%lld\n",(ans+mod)%mod);
}
}

最新文章

  1. 如何在Texstudio内加载语法检查词典?
  2. SpringMVC中servletFileUpload.parseRequest(request)解析为空获取不到数据问题
  3. 制作IOS企业版App网页扫描二维码下载安装
  4. CoreData总结
  5. gvim配置
  6. MVC控制器常用方法返回类型
  7. SQL SERVER 2008安装错误(is not a valid login or you do have permission)
  8. OpenGL: 渲染管线理论
  9. 下载APP 2个二维码合并到一个二维码
  10. jsp执行数据库查询并分页
  11. 10个SQL注入工具(转载)
  12. uva10791
  13. hdu 2149 Public Sale 简单博弈
  14. C++中运行外部程序
  15. OpenFlow协议1.0及1.3版本分析
  16. 设置Mac 终端走代理
  17. 记录下用axios遇到的问题
  18. AtCoder Regular Contest 082 (ARC082) E - ConvexScore 计算几何 计数
  19. CentOS6.5安装kafka-2.10-0.8.2(单机)
  20. 用Javascript,DHTML控制表格的某一列的显示与隐藏

热门文章

  1. Linux添加目录到环境变量以及添加Sublime Text到环境变量
  2. LeetCode 47——全排列 II
  3. 阶段1 语言基础+高级_1-3-Java语言高级_07-网络编程_第1节 网络通信概述_5_端口号
  4. iView 实战系列教程(21课时)_2.iView 实战教程之导航、路由、鉴权篇
  5. Delphi 跨单元进入(访问)类的私有成员,protected ,private部分
  6. Java ——修饰符 包 Bean
  7. Python笔记(二十二)_魔法方法_基本魔法方法
  8. mysql 主从 设置
  9. HDFS-HA高可用工作机制
  10. &lt;每日一题&gt; Day4:CodeForces-1042A.Benches(二分 + 排序)