【HAOI2018】染色【反向二项式反演】【NTT卷积】

传送门

题意:NNN个位置染MMM种颜色,恰好出现SSS次的颜色数量恰好为kkk时的愉悦度为wkw_kwk,求所有方案的愉悦度之和。对100453580910045358091004535809取模。

N≤1e7N \leq 1e7N1e7,M≤1e5M \leq 1e5M1e5,S≤150S \leq 150S150

本题的恶心之处在于满足颜色数量要求恰好为kkk

试着求求至少为kkk的方案数。

一个显然的想法:钦点kkk种颜色刚好填SSS次,剩下位置随便填其他颜色。

即:

f(k)=(Mk)(NS)(N−SS)......(N−(k−1)SS)(M−k)N−kSf(k)=\binom{M}{k}\binom{N}{S}\binom{N-S}{S}......\binom{N-(k-1)S}{S}(M-k)^{N-kS}f(k)=(kM)(SN)(SNS)......(SN(k1)S)(Mk)NkS

=(Mk)N!S!(N−S)!(N−S)!S!(N−2S)!......[N−(k−1)S]!S!(N−kS)!(M−k)N−kS=\binom{M}{k}\frac{N!}{S!(N-S)!}\frac{(N-S)!}{S!(N-2S)!}......\frac{[N-(k-1)S]!}{S!(N-kS)!}(M-k)^{N-kS}=(kM)S!(NS)!N!S!(N2S)!(NS)!......S!(NkS)![N(k1)S]!(Mk)NkS

=(Mk)N!(M−k)N−kS(S!)k(N−kS)!=\binom{M}{k}\frac{N!(M-k)^{N-kS}}{(S!)^k(N-kS)!}=(kM)(S!)k(NkS)!N!(Mk)NkS

虽然很鬼畜,但是可以快速求出来。

但这个求出来是假的,因为剩下位置填其他颜色的时候一不小心就刚好填了SSS次,但这个在之后会重复计算。也就是说,一个方案会计算(实际刚好为S的颜色数k)\binom{实际刚好为S的颜色数}{k}(kS)次。所以似乎没有组合意义。

不过还是能用的。

设答案即恰好kkk种的方案数g(k)g(k)g(k)

为了描述方便,设n=min(n,N/S)n=min(n,N/S)n=min(n,N/S),即填SSS个的颜色最多多少个。

得到关系式

f(k)=∑i=kn(ik)g(i)f(k)=\sum_{i=k}^n\binom{i}{k}g(i)f(k)=i=kn(ki)g(i)

另一个方向的二项式反演?

g(k)=∑i=kn[i−k=0](ik)g(i)g(k)=\sum_{i=k}^n[i-k=0]\binom{i}{k}g(i)g(k)=i=kn[ik=0](ki)g(i)

g(k)=∑i=kn∑j=0i−k(−1)j(i−kj)(ik)g(i)g(k)=\sum_{i=k}^n\sum_{j=0}^{i-k}(-1)^j\binom{i-k}{j}\binom{i}{k}g(i)g(k)=i=knj=0ik(1)j(jik)(ki)g(i)

你会发现用之前的路子行不通了

为了找到这一步怎么推,你可以把结论代回去,然后你得到了这个东西:

(i−kj)(ik)=(ij+k)(j+kk)\binom{i-k}{j}\binom{i}{k}=\binom{i}{j+k}\binom{j+k}{k}(jik)(ki)=(j+ki)(kj+k)

理性证明:

(i−k)!j!(i−j−k)!i!k!(i−k)!=i!(j+k)!(i−j−k)!(j+k)!j!k!\frac{(i-k)!}{j!(i-j-k)!}\frac{i!}{k!(i-k)!}=\frac{i!}{(j+k)!(i-j-k)!}\frac{(j+k)!}{j!k!}j!(ijk)!(ik)!k!(ik)!i!=(j+k)!(ijk)!i!j!k!(j+k)!

感性证明:

iii个数选j+kj+kj+k个数再选kkk个数,等价于直接选kkk个再在剩下的选出jjj个作为中间商赚差价

g(k)=∑i=kn∑j=0i−k(−1)j(ij+k)(j+kk)g(i)g(k)=\sum_{i=k}^n\sum_{j=0}^{i-k}(-1)^j\binom{i}{j+k}\binom{j+k}{k}g(i)g(k)=i=knj=0ik(1)j(j+ki)(kj+k)g(i)

jjj加上kkk

g(k)=∑i=kn∑j=ki(−1)j−k(ij)(jk)g(i)g(k)=\sum_{i=k}^n\sum_{j=k}^{i}(-1)^{j-k}\binom{i}{j}\binom{j}{k}g(i)g(k)=i=knj=ki(1)jk(ji)(kj)g(i)

交换求和顺序

g(k)=∑j=kn(−1)j−k(jk)∑i=jn(ij)g(i)g(k)=\sum_{j=k}^n(-1)^{j-k}\binom{j}{k}\sum_{i=j}^{n}\binom{i}{j}g(i)g(k)=j=kn(1)jk(kj)i=jn(ji)g(i)

g(k)=∑j=kn(−1)j−k(jk)f(j)g(k)=\sum_{j=k}^n(-1)^{j-k}\binom{j}{k}f(j)g(k)=j=kn(1)jk(kj)f(j)

g(k)=∑i=kn(−1)i−k(ik)f(i)g(k)=\sum_{i=k}^n(-1)^{i-k}\binom{i}{k}f(i)g(k)=i=kn(1)ik(ki)f(i)

回到之前的问题,我们现在要求所有g(k)g(k)g(k)

套路性的拆组合数

g(k)=∑i=kn(−1)i−ki!k!(i−k)!f(i)g(k)=\sum_{i=k}^n(-1)^{i-k}\frac{i!}{k!(i-k)!}f(i)g(k)=i=kn(1)ikk!(ik)!i!f(i)

k!g(k)=∑i=kn(−1)i−k(i−k)!i!f(i)k!g(k)=\sum_{i=k}^n\frac{(-1)^{i-k}}{(i-k)!}i!f(i)k!g(k)=i=kn(ik)!(1)iki!f(i)

出现了!卷积君!

翻一下i!f(i)i!f(i)i!f(i)

k!g(k)=∑i=kn(−1)i−k(i−k)!(n−i)!f(n−i)k!g(k)=\sum_{i=k}^n\frac{(-1)^{i-k}}{(i-k)!}(n-i)!f(n-i)k!g(k)=i=kn(ik)!(1)ik(ni)!f(ni)

卷出来是n−kn-knk,再翻一下除以k!k!k!就是答案。

#include <iostream>
#include <cstdio>
#include <cstring>
#include <cctype>
#include <algorithm>
#define MAXN 262144
#define MAXM 10000005
using namespace std;
const int MOD=1004535809;//479*2^21+1
typedef long long ll;
int fac[MAXM],finv[MAXM];
inline int qpow(int a,int p)
{int ans=1;while (p){if (p&1) ans=(ll)ans*a%MOD;a=(ll)a*a%MOD;p>>=1;}return ans;
}
#define inv(x) qpow(x,MOD-2)
inline int add(const int& x,const int& y){return x+y>=MOD? x+y-MOD:x+y;}
inline int dec(const int& x,const int& y){return x<y? x-y+MOD:x-y;}
int r[MAXN],rt[2][22];
inline void init(const int& l){for (int i=0;i<(1<<l);i++) r[i]=(r[i>>1]>>1)|((i&1)<<(l-1));}
void NTT(int* a,int l,int type)
{int lim=1<<l;for (int i=0;i<lim;i++) if (i<r[i]) swap(a[i],a[r[i]]);for (int L=0;L<l;L++){int mid=1<<L,len=mid<<1;int Wn=rt[type][L+1];for (int s=0;s<lim;s+=len)for (int k=0,w=1;k<mid;k++,w=(ll)w*Wn%MOD){int x=a[s+k],y=(ll)w*a[s+mid+k]%MOD;a[s+k]=add(x,y);a[s+mid+k]=dec(x,y);}}if (type){int t=inv(lim);for (int i=0;i<lim;i++) a[i]=(ll)a[i]*t%MOD;}
}
int f[MAXN],g[MAXN];
int main()
{rt[0][21]=qpow(3,479);rt[1][21]=inv(rt[0][21]);for (int i=20;i>=0;i--) {rt[0][i]=(ll)rt[0][i+1]*rt[0][i+1]%MOD;rt[1][i]=(ll)rt[1][i+1]*rt[1][i+1]%MOD;}int n,m,s;scanf("%d%d%d",&n,&m,&s);int lim=min(m,n/s),N=max(n,m);fac[0]=1;for (int i=1;i<=N;i++) fac[i]=(ll)fac[i-1]*i%MOD;finv[N]=inv(fac[N]);for (int i=N-1;i>=0;i--) finv[i]=(ll)finv[i+1]*(i+1)%MOD;int l=0;while ((1<<l)<=(lim<<1)) ++l;init(l);for (int i=0;i<=lim;i++) f[i]=(ll)fac[m]*fac[n]%MOD*qpow(m-i,n-s*i)%MOD*finv[m-i]%MOD*finv[n-s*i]%MOD*qpow(finv[s],i)%MOD;reverse(f,f+lim+1);for (int i=0;i<=lim;i++) g[i]=((i&1)? MOD-finv[i]:finv[i]);NTT(f,l,0);NTT(g,l,0);for (int i=0;i<(1<<l);i++) f[i]=(ll)f[i]*g[i]%MOD;NTT(f,l,1);reverse(f,f+lim+1);for (int i=0;i<=lim;i++) f[i]=(ll)f[i]*finv[i]%MOD;int ans=0;for (int i=0;i<=lim;i++){int w;scanf("%d",&w);ans=add(ans,(ll)w*f[i]%MOD);}
//	for (int i=0;i<=lim;i++)
//	{
//		int sum=0;
//		for (int j=i;j<=lim;j++) sum=add(sum,(ll)(((j-i)&1)? MOD-fac[j]:fac[j])*finv[j-i]%MOD*f[j]%MOD);
//		int w;
//		scanf("%d",&w);
//		ans=add(ans,(ll)w*finv[i]%MOD*sum%MOD);
//	}printf("%d\n",ans);return 0;
} 

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

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

相关文章

从严治码-别人在项目中下毒,我该怎么治?

01 从软考说起从4月份开始&#xff0c;由于备考《系统集成项目管理工程师》的原因&#xff0c;博客没有持续更新&#xff0c;在上半年考试结束之后&#xff0c;又对项目进行了一些收尾的工作。下面就这段时间的学习作一个记录和总结吧。在学习的过程中&#xff0c;提炼了一些自…

CodeCraft-20 (Div. 2) C. Primitive Primes 思维 + 数论

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 给你两个长度分别为n,mn,mn,m的多项式&#xff0c;将他们乘起来&#xff0c;问系数modp0\bmod p 0modp0的项的指数是多少&#xff0c;两个多项式所有项的系数gcd1gcd1gcd1。 n,m<1e6n,m<1e6n,m<1e6…

【LOJ166】拉格朗日插值2【拉格朗日插值】【NTT卷积】

传送门 题意&#xff1a;给定n,m,f(0),f(1),......,f(n)n,m,f(0),f(1),... ...,f(n)n,m,f(0),f(1),......,f(n)&#xff0c;求f(m),f(m1),......f(mn)f(m),f(m1),... ...f(mn)f(m),f(m1),......f(mn) 模998244353998244353998244353 n≤100000,m≤1e8,n<mn \leq 100000,m \…

学习MVVM设计模式后第一次用于生产

WPF的MVVM设计模式从winform转变到WPF的过程&#xff0c;难点主要还是在MVVM的设计模式。当然&#xff0c;如果依然采用winform的涉及方式&#xff0c;在每个控件背后绑定事件的方式运用在wpf中&#xff0c;依然可行&#xff0c;但是假如GUI改版&#xff0c;其背后绑定的特别为…

剑指 Offer 14- II. 剪绳子 II

给你一根长度为 n 的绳子&#xff0c;请把绳子剪成整数长度的 m 段&#xff08;m、n都是整数&#xff0c;n>1并且m>1&#xff09;&#xff0c;每段绳子的长度记为 k[0],k[1]...k[m - 1] 。请问 k[0]*k[1]*...*k[m - 1] 可能的最大乘积是多少&#xff1f;例如&#xff0c;…

Educational Codeforces Round 108 (Rated for Div. 2) D. Maximum Sum of Products 思维 + dp

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 给你两个长度为nnn的数组a,ba,ba,b&#xff0c;你可以至多反转一段连续区间&#xff0c;求∑i1nai∗bi\sum _{i1}^n a_i*b_i∑i1n​ai​∗bi​最大是多少。 n<5e3n<5e3n<5e3 思路&#xff1a; 首…

【CF1215E】Marbles【状压DP】

传送门 题意&#xff1a;给一个长为NNN的序列aaa&#xff0c;每次操作交换两个相邻位置&#xff0c;求最少操作次数使得所有相同的值连成一片。 N≤400000N \leq 400000N≤400000,ai≤20a_i \leq20ai​≤20 我们发现aia_iai​很小&#xff0c;盲猜单独考虑 我们重新确认一个…

netcore mvc快速开发系统(菜单,角色,权限[精确到按钮])开源

基于netcore2.0 mvc 开发的 快速搭建具有如下特色的后台管理系统用户管理菜单管理角色管理权限管理[精确到按钮]&#xff09;代码生成器代码克隆到本地 用vs2017或以上版本 打开工程。项目结构如下&#xff1a;找到DbModel下面的初始化db脚本里面包含4张表的schema和初始化数据…

剑指 Offer 25. 合并两个排序的链表

输入两个递增排序的链表&#xff0c;合并这两个链表并使新链表中的节点仍然是递增排序的。 示例1&#xff1a; 输入&#xff1a;1->2->4, 1->3->4 输出&#xff1a;1->1->2->3->4->4 限制&#xff1a; 0 < 链表长度 < 1000 思路&#xff1a;和…

Codeforces Round #720 (Div. 2) C. Nastia and a Hidden Permutation 交互

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 给你一个序列ppp长度nnn&#xff0c;每次可以执行两个种询问&#xff1a; t1max(min(x,pi),min(x1,pj))t1\ \ max(min(x,p_i),min(x1,p_j))t1 max(min(x,pi​),min(x1,pj​)) t2min(max(x,pi),max(x1,pj))t…

【NOIP2018】赛道修建【二分】【树形dp】【multiset】【贪心】

传送门 题意&#xff1a;给一棵带边权的树&#xff0c;求MMM条没有公共边的路径的最小长度的最大值。 N≤50000N \leq50000N≤50000 抛开NOIP不谈&#xff0c;其实这题本身出的很好 显然二分 问题转化成了“最多可以找多少条长度不小于kkk的路径” 递归处理 对于每个结点…

asp.net core 系列之Startup

这篇文章简单记录 ASP.NET Core中 &#xff0c;startup类的一些使用。一.前言在 Startup类中&#xff0c;一般有两个方法&#xff1a;ConfigureServices 方法: 用来配置应用的 service 。 Configure 方法&#xff1a;创建应用的请求处理管道它们都在应用启动时&#xff0c;被AS…

可持久化Splay 学习笔记

可持久化Splay是怎么回事呢&#xff1f;Splay相信大家都很熟悉&#xff0c;但是可持久化Splay是怎么回事呢&#xff0c;下面就让小编带大家一起了解吧。   可持久化Splay&#xff0c;其实就是将Splay持久化一下&#xff0c;大家可能会很惊讶Splay怎么可以持久化呢&#xff1f…

译 | .NET Core 基础架构进化之路(一)

原文&#xff1a;Matt Mitchell翻译&#xff1a;Edi Wang随着 .NET Core 3.0 Preview 6 的推出&#xff0c;我们认为简要了解一下我们基础设施系统的历史以及过去一年左右所做的重大改进会很有用。如果您对构建基础结构感兴趣&#xff0c;或者想要了解我们如何构建与 .NET Core…

【洛谷P4705】玩游戏【二项式定理】【NTT卷积】【生成函数】【分治NTT】【函数求导】【多项式对数】

传送门 题意&#xff1a;给定长度为N,MN,MN,M的序列a,ba,ba,b和ttt&#xff0c;随机选取x∈[1,N],y∈[1,M]x \in[1,N],y\in[1,M]x∈[1,N],y∈[1,M],对于i1,2,...,t,i 1,2,...,t,i1,2,...,t,求(axby)i(a_xb_y)^i(ax​by​)i的期望 N,M,t≤100000N,M,t \leq100000N,M,t≤100000 …

剑指 Offer 27. 二叉树的镜像

思路&#xff1a;递归 /*** Definition for a binary tree node.* struct TreeNode {* int val;* TreeNode *left;* TreeNode *right;* TreeNode(int x) : val(x), left(NULL), right(NULL) {}* };*/ class Solution { public:TreeNode* mirrorTree(TreeNode*…

Codeforces Round #626 (Div. 2) D. Present 按位贡献 + 快排新姿势

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 给你一个长度为nnn的序列aaa&#xff0c;让你计算 n≤4e5,a≤1e7n\le 4e5,a\le 1e7n≤4e5,a≤1e7 思路&#xff1a; 首先这个式子是n2n^2n2的&#xff0c;显然不能直接算&#xff0c;并且异或没有分配律&…

.NET开发框架(一)-框架介绍与视频演示

本文主要介绍一套基于.NET CORE的SPA高并发、高可用的开发框架.我们暂且称它为&#xff1a;&#xff08;让你懂.NET&#xff09;开发框架。以此为主线&#xff0c;陆续编写教程&#xff0c;讲述如何构建高并发、高可用的框架。&#xff08;欢迎转载与分享&#xff09;它标准化了…

【CF700E】Cool Slogans【后缀自动机】【可持久化线段树合并】【树上倍增】

传送门 题意&#xff1a;给定字符串SSS&#xff0c;求一堆字符串s1,s2,s3,...,sks_1,s_2,s_3,...,s_ks1​,s2​,s3​,...,sk​&#xff0c;满足s1s_1s1​是SSS的子串&#xff0c;且sis_isi​在si−1s_{i-1}si−1​中至少出现两次&#xff0c;最大化kkk ∣S∣≤200000|S| \leq …

P1377 [TJOI2011]树的序 笛卡尔树优化建树

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 给你一棵二叉树的生成序列&#xff0c;让你输出一个字典序最小的序列&#xff0c;使其生成的二叉树与原来的二叉树相同。 思路&#xff1a; 首先想到暴力建树&#xff0c;让后输出先序遍历即可&#xff0c…