P6076-[JSOI2015]染色问题【组合数学,容斥】

正题

题目链接:https://www.luogu.com.cn/problem/P6076


题目大意

给出n∗mn*mnm的网格,ccc种颜色涂色要求

  1. 每个格子可以染色也可以不染
  2. 每一行每一列至少有一个格子被染
  3. 每个颜色至少用一次

1≤n,m,c≤4001\leq n,m,c\leq 4001n,m,c400


解题思路

一个比较简单的方法就是容斥,枚举有多少染色的和不染色的行列,和枚举使用的颜色个数
∑i=0c∑j=0n∑k=0m(ci)(nj)(mk)(i+1)j+k(−1)c+n+m−i−j−k\sum_{i=0}^c\sum_{j=0}^n\sum_{k=0}^m\binom{c}{i}\binom nj\binom mk(i+1)^{j+k}(-1)^{c+n+m-i-j-k}i=0cj=0nk=0m(ic)(jn)(km)(i+1)j+k(1)c+n+mijk

这样预处理就是O(nmc)O(nmc)O(nmc)的,但是可以做到更快。

fif_ifi表示最多染了iii种颜色的方案,那么久只需要满足第二个条件了。第二个条件可以用一个容斥搞定,考虑枚举多少行没染
fk=∑i=1n(−1)n−i((k+1)m−1)if_k=\sum_{i=1}^n(-1)^{n-i}((k+1)^{m}-1)^ifk=i=1n(1)ni((k+1)m1)i
这样预处理就可以做到O(nc)O(nc)O(nc)

这里写的是第一种,因为比较懒


code

#include<cstdio>
#include<cstring>
#include<algorithm>
#define ll long long
using namespace std;
const ll N=410,P=1e9+7;
ll n,m,c,C[N][N],pw[N*N],ans;
signed main()
{scanf("%lld%lld%lld",&n,&m,&c);C[0][0]=1;for(ll i=1;i<N;i++)for(ll j=0;j<N;j++)C[i][j]=(C[i-1][j]+(j?C[i-1][j-1]:0))%P; for(ll i=0;i<=c;i++){pw[0]=1;for(ll j=1;j<=n*m;j++)pw[j]=pw[j-1]*(i+1)%P;for(ll j=0;j<=n;j++)for(ll k=0;k<=m;k++){ll f=(c-i)+(n-j)+(m-k);if(f&1)f=-1;else f=1;(ans+=f*C[n][j]*C[m][k]%P*C[c][i]%P*pw[j*k]%P)%=P;}}printf("%lld\n",(ans+P)%P);return 0;
}

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/319146.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

E - Escape from the Island(最短路+dp)

E - Escape from the Island 大佬题解&#xff0c;码风真的爱了 状态表式&#xff1a;f(u,j)f(u,j)f(u,j)当前在uuu点&#xff0c;已经划了jjj步时离终点的最短距离 状态转移&#xff1a; 主动划一步min&#xff0c;转移到下一个点vvv f(u,j)f(v,j1)1,(u,v)∈E,(u,v)∈Ef(u,j…

国内开源社区巨作AspectCore-Framework入门

前些天和张队(善友),lemon(浩洋),斌哥(项斌)等MVP大咖一块儿吃饭,大家聊到了lemon名下的AOP这个项目,我这小白听得一脸懵逼,后面回来做了一下功课,查了下资料,在lemon的Github上把这个项目学习了一下,收获颇丰,让我这个没有接触过AOP的Coder叹为观止,陷入了对lemon的深深崇拜,在…

【数论】Product(P5221)

正题 P5221 题目大意 求 ∏i1n∏j1nlcm(i,j)gcd(i,j)\prod_{i1}^n\prod_{j1}^n\frac{lcm(i,j)}{gcd(i,j)}i1∏n​j1∏n​gcd(i,j)lcm(i,j)​ 解题思路 ∏i1n∏j1nlcm(i,j)gcd(i,j)\prod_{i1}^n\prod_{j1}^n\frac{lcm(i,j)}{gcd(i,j)}i1∏n​j1∏n​gcd(i,j)lcm(i,j)​ ∏i1n∏…

STL的可持久化数组

rope是一种可持久化数组&#xff0c;可持久化平衡树&#xff0c;采用块状链表实现 #include<ext/rope> using namespace __gnu_cxx; rope<int>a;rope就是一个用可持久化平衡树实现的“重型”string 库中模板计算基本和string一样简单 各种操作的复杂度都是O(log…

HarmonyOs4.0基础(一)

目录 一、HarmonyOs系统定义 1.1系统的技术特性(三大特征) 1.1.1、硬件互助、资源共享 1.1.2、一次开发、多端部署(面向开发者) 1.1.3、统一OS&#xff0c;弹性部署(支持多种API&#xff1a;ArkTs、JS、C/C、Java) 1.2、系统的技术架构 二、Harmony OS项目搭建 2.1、(D…

P4643-[国家集训队]阿狸和桃子的游戏【结论】

正题 题目链接:https://www.luogu.com.cn/problem/P4643 题目大意 给出nnn个点mmm条边的无向图&#xff0c;两个人轮流选择一个未被选择的点加入点集。 然后每个人的权值为选出的点的导出子图点权加边权和。 两个人都希望自己的权值减去对方的权值最大 求先手的权值减去后…

C. Code a Trie(Trie+dfs+贪心)

C. Code a Trie 大佬题解&#xff0c;代码基本就是抄的 对于每一个值计算所有串的LCA&#xff0c;也就是最长公共前缀&#xff0c;将该节点&#xff08;Trie树的节点&#xff09;标记&#xff0c;对于这些字符串在LCA下面的点一定不存在&#xff08;如果存在他们不会返回相同…

swagger文档转换为WebApiClient声明式代码

1 swagger简介Swagger是一个规范且完整的框架&#xff0c;提供描述、生产、消费和可视化RESTful Web Service。其核心是使用json来规范描述RESTful接口&#xff0c;另外有提供UI来查看接口说明&#xff0c;并有一套生成不同语言的客户端调用代码生成器。1.1 对Api提供者自顶向下…

Musical Theme pku1743 (后缀数组)

Musical Theme(后缀数组) 题意&#xff1a; n个数&#xff0c;选取一段子序列&#xff0c;满足以下条件&#xff1a; 1.长度至少为5 2.在数列中其他位置出现过(允许转置) 3.与其他位置出现的不重叠 转置&#xff1a;将恒定的正或负值添加到子序列上 例如&#xff1a; n个数为…

各种有用的东西留言板

目录各种链接各种OJ各种有用网站模拟赛排名表格制作其他STLO2优化&#xff1a;快读fread黑科技各种链接 各种OJ ssl 洛谷 纪中 牛客 PTA Vijos JoyOI 51Nod 大视野 ssl比赛 NOI.AC AcWing USACO CodeVS 北京大学 浙江大学 Comet OJ 杭州电子科技大学 《算法竞赛进阶指南》题库…

P3180-[HAOI2016]地图【圆方树,莫队,分块】

正题 题目链接:https://www.luogu.com.cn/problem/P3180 题目大意 nnn个点mmm条边的一个仙人掌&#xff0c;有点权。 QQQ次询问给出op,x,yop,x,yop,x,y&#xff0c;封闭111到xxx号点的所有简单路径后xxx能到达的点的点权中&#xff0c;小于yyy且出现次数为奇数/偶数的权值数目…

M. Monster Hunter(树形dp)

M. Monster Hunter 才知道原来树形dp是三维的&#xff0c;一直没有学会过&#xff0c;感谢大佬的文章&#xff01;算法进阶—理解树形背包问题 状态表示&#xff1a;fi,k,j,{0/1}f_{i,k,j,\{0/1\}}fi,k,j,{0/1}​以iii为根的子树&#xff0c;考虑到第kkk个儿子时&#xff0c;…

Another Blog

说明&#xff1a; 此博客用于平时做题 但为了整理知识&#xff0c;用博客园来存知识 博客园网址&#xff1a; https://www.cnblogs.com/ssllyf

KubeCon+CloudNativeCon首秀中国!

2018年11月13-15日&#xff0c;全球顶级的Kubernetes官方技术论坛KubeConCloudNativeCon将首次登陆中国&#xff0c;此次活动由云原生计算基金会&#xff08;CNCF&#xff09;主办&#xff0c;在上海跨国采购会展中心隆重举行。KubeCon CloudNativeConKubeConCloudNativeCon 是…

可持久化(一)

参考博客 可持久化数据结构&#xff1a;可以保留每一个历史版本&#xff0c;若所有版本都既可以访问又可以修改&#xff0c;成为完全可持久化&#xff08;可以回滚到某个历史版本&#xff09; 时间线&#xff1a; 可持久化线段树 可持久化下标线段树 题目&#xff1a; 模板…

【模板】高斯消元

ACM模板 目录高斯消元解线性方程组异或方程组bitset优化异或方程组高斯消元解线性方程组 int a[N][N]输入矩阵&#xff0c;nnn行&#xff0c;n1n1n1列&#xff0c;下标从0开始 第n1n1n1列表示方程右边的值&#xff08;n行即n个方程&#xff0c;n列即n个未知数&#xff09; int…

各种模板(数据结构图论)

文章目录数据结构LCT线段树线段树分治树状数组图论Tarjan静态仙人掌最小生成树最短路-Floyd最短路-Dijkstra最短路-Bellman-Ford最短路-SPFA数据结构 LCT #include<cstdio> #include<cstring> #include<iostream> #include<algorithm> #define ll lo…

P6620-[省选联考2020A卷]组合数问题【组合数学,斯特林数】

正题 题目链接:https://www.luogu.com.cn/problem/P6620 题目大意 给出n,x,p,mn,x,p,mn,x,p,m和一个mmm次多项式fff求 ∑k0nf(k)xk(nk)\sum_{k0}^nf(k)\times x^k\times \binom{n}{k}k0∑n​f(k)xk(kn​) 答案对ppp取模。 1≤n≤109,1≤m≤10001\leq n\leq 10^9,1\leq m\leq …

ASP.NET Core中使用GraphQL - 第一章 Hello World

前言你是否已经厌倦了REST风格的API? 让我们来聊一下GraphQL。 GraphQL提供了一种声明式的方式从服务器拉取数据。你可以从GraphQL官网中了解到GraphQL的所有优点。在这一系列博客中&#xff0c;我将展示如何在ASP.NET Core中集成GraphQL, 并使用GraphQL作为你的API查询语言。…

L2-005 集合相似度 (25分)

L2-005 集合相似度 (25分) 给定两个整数集合&#xff0c;它们的相似度定义为&#xff1a;N​c/N​t​​ 100%。其中N​c 是两个集合都有的不相等整数的个数&#xff0c;N​t是两个集合一共有的不相等整数的个数。你的任务就是计算任意一对给定集合的相似度。 输入格式&#xf…