【洛谷P4709】信息传递（置换）（组合数学）（多项式Exp）

最新推荐文章于 2023-10-24 22:30:48 发布
最新推荐文章于 2023-10-24 22:30:48 发布
阅读量333 收藏 1
点赞数
分类专栏： 生成函数 多项式 置换群
版权声明：本文为博主原创文章，遵循 CC 4.0 BY 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/zxyoi_dreamer/article/details/97683487
版权
多项式 同时被 3 个专栏收录
85 篇文章 2 订阅
订阅专栏
生成函数
37 篇文章 1 订阅
订阅专栏
置换群
12 篇文章 0 订阅
订阅专栏

传送门

题解：

很早之前看到过这道题，当时连置换是什么都不知道。。。

首先考虑置换 g g g中的某一个长为 S S S的循环环，学过群论的应该知道，在置换自乘 n n n次之后，这个循环会变为 g c d ( n , S ) gcd(n,S) gcd(n,S)个循环。

考虑将 t t t个长为 S S S的循环拼接成一个，使得若干次自乘之后这个循环断裂为 t t t个长为 S S S的循环。翻过任何一本群论教材就知道，最后在同一循环中的元素是那些只与初始下标 % S \%S %S的值有关，所以这一部分的方案数为 S t − 1 ( t − 1 ) ! S^{t-1}(t-1)! St−1(t−1)!，其实就是考虑一个圆排列。

所有长度为 S S S的循环显然可以放在一起处理。

考虑 e x p exp exp的组合意义，显然我们对于每个长度求出拿 t t t个该长度的串来拼接成一个的方案数的 E G F EGF EGF，设为 A ( x ) A(x) A(x)，则拿出 c n t [ S ] cnt[S] cnt[S]个循环，拼接成若干个的方案数的生成函数为 exp ⁡ ( A ( x ) ) \exp(A(x)) exp(A(x))，复杂度 O ( n log ⁡ n ) O(n\log n) O(nlogn)，而且不好卡满。

然而这个玩意可以直接 D P DP DP，复杂度 O ( n 4 3 ) O(n^{\frac{4}{3}}) O(n34​)，同样，几乎卡不满

代码（DP）：

#include<bits/stdc++.h>
#define ll long long
#define re register
#define gc get_char
#define cs const

namespace IO{
	inline char get_char(){
		static cs int Rlen=1<<22|1;
		static char buf[Rlen],*p1,*p2;
		return (p1==p2)&&(p2=(p1=buf)+fread(buf,1,Rlen,stdin),p1==p2)?EOF:*p1++;
	}
	
	template<typename T>
	inline T get(){
		char c;
		while(!isdigit(c=gc()));T num=c^48;
		while(isdigit(c=gc()))num=(num+(num<<2)<<1)+(c^48);
		return num;
	}
	inline int getint(){return get<int>();}
}
using namespace IO;

using std::cerr;
using std::cout;

cs int mod=998244353;
inline int add(int a,int b){return (a+=b)>=mod?a-mod:a;}
inline int dec(int a,int b){return (a-=b)<0?a+mod:a;}
inline int mul(int a,int b){ll r=(ll)a*b;return r>=mod?r%mod:r;}
inline int power(int a,int b,int res=1){
	for(;b;b>>=1,a=mul(a,a))(b&1)&&(res=mul(res,a));
	return res;
}
inline void Inc(int &a,int b){(a+=b)>=mod&&(a-=mod);}
inline void Dec(int &a,int b){(a-=b)<0&&(a+=mod);}

inline int gcd(int a,int b){
#define ctz __builtin_ctz
	int shift=ctz(a|b);
	for(b>>=ctz(b);a;a-=b)if((a>>=ctz(a))<b)std::swap(a,b);
	return b<<=shift;
}

cs int N=1e5+5;

int n;

int fac[N],ifac[N];
inline void init(){
	fac[0]=fac[1]=ifac[0]=1;
	for(int re i=2;i<=n;++i)fac[i]=mul(fac[i-1],i);
	ifac[n]=power(fac[n],mod-2);
	for(int re i=n-1;i;--i)ifac[i]=mul(ifac[i+1],i+1);
}

inline int C(int n,int m){return (n>=0&&m>=0&&n>=m)?mul(fac[n],mul(ifac[m],ifac[n-m])):0;}

std::vector<int> factor;

int to[N],cnt[N];
bool vis[N];

int dp[N],coef[N];

signed main(){
#ifdef zxyoi
	freopen("message.in","r",stdin);
#endif
	n=getint();init();
	for(int re i=1;i<=n;++i)to[i]=getint();
	for(int re i=1;i<=n;++i)if(!vis[i]){
		int len=0,x=i;
		while(!vis[x]){
			vis[x]=true;
			++len;x=to[x];
		}
		++cnt[len];
	}
	for(int re i=1;i*i<=n;++i)
	if(n%i==0){
		factor.push_back(i);
		if(i*i!=n)factor.push_back(n/i);
	}
	std::sort(factor.begin(),factor.end());
	int ans=1;
	for(int re i=1;i<=n&&ans;++i)if(cnt[i]){
		coef[1]=1;
		for(int re j=2;j<=cnt[i];++j)coef[j]=mul(coef[j-1],mul(i,j-1));
		memset(dp+1,0,sizeof(int)*cnt[i]);
		dp[0]=1;
		for(int re j=1;j<=cnt[i];++j){
			for(int re t:factor){
				if(t>j)break;
				if(gcd(n,i*t)!=t)continue;
				Inc(dp[j],mul(dp[j-t],mul(C(j-1,t-1),coef[t])));
			}
		}
		ans=mul(ans,dp[cnt[i]]);
	}
	cout<<ans<<"\n";
	return 0;
}

代码exp：

#include<bits/stdc++.h>
#define ll long long
#define re register
#define gc get_char
#define cs const

namespace IO{
	inline char get_char(){
		static cs int Rlen=1<<22|1;
		static char buf[Rlen],*p1,*p2;
		return (p1==p2)&&(p2=(p1=buf)+fread(buf,1,Rlen,stdin),p1==p2)?EOF:*p1++;
	}
	
	template<typename T>
	inline T get(){
		char c;
		while(!isdigit(c=gc()));T num=c^48;
		while(isdigit(c=gc()))num=(num+(num<<2)<<1)+(c^48);
		return num;
	}
	inline int getint(){return get<int>();}
}
using namespace IO;

using std::cerr;
using std::cout;

inline int gcd(int a,int b){
#define ctz __builtin_ctz
	int shift=ctz(a|b);
	for(b>>=ctz(b);a;a-=b)if((a>>=ctz(a))<b)std::swap(a,b);
	return b<<=shift;
}

cs int mod=998244353;
inline int add(int a,int b){return (a+=b)>=mod?a-mod:a;}
inline int dec(int a,int b){return (a-=b)<0?a+mod:a;}
inline int mul(int a,int b){ll r=(ll)a*b;return r>=mod?r%mod:r;}
inline int power(int a,int b,int res=1){
	for(;b;b>>=1,a=mul(a,a))(b&1)&&(res=mul(res,a));
	return res;
}
inline void Inc(int &a,int b){(a+=b)>=mod&&(a-=mod);}
inline void Dec(int &a,int b){(a-=b)<0&&(a+=mod);}

cs int SIZE=1<<19|1,bit=19;

typedef std::vector<int> Poly;

int r[SIZE],*w[bit+1];
inline void init_NTT(){
	for(int re i=1;i<=bit;++i)w[i]=new int[1<<i-1];
	int wn=power(3,(mod-1)>>bit);
	w[bit][0]=1;
	for(int re i=1;i<(1<<bit-1);++i)w[bit][i]=mul(w[bit][i-1],wn);
	for(int re i=bit-1;i;--i)
	for(int re j=0;j<(1<<i-1);++j)w[i][j]=w[i+1][j<<1];
}
inline void NTT(Poly &A,int len,int typ){
	for(int re i=0;i<len;++i)if(i<r[i])std::swap(A[i],A[r[i]]);
	for(int re i=1,t=1;i<len;i<<=1,++t)
	for(int re j=0;j<len;j+=i<<1)
	for(int re k=0;k<i;++k){
		int x=A[j+k],y=mul(A[j+k+i],w[t][k]);
		A[j+k]=add(x,y),A[j+k+i]=dec(x,y);
	}
	if(typ==-1){
		std::reverse(A.begin()+1,A.begin()+len);
		for(int re i=0,inv=power(len,mod-2);i<len;++i)A[i]=mul(A[i],inv);
	}
}
inline void init_rev(int len){
	for(int re i=0;i<len;++i)r[i]=r[i>>1]>>1|((i&1)?len>>1:0);
}

int fac[SIZE],ifac[SIZE],inv[SIZE];
inline void init_inv(){
	fac[0]=fac[1]=ifac[0]=ifac[1]=inv[0]=inv[1]=1;
	for(int re i=2;i<SIZE;++i){
		fac[i]=mul(fac[i-1],i);
		inv[i]=mul(inv[mod%i],mod-mod/i);
		ifac[i]=mul(ifac[i-1],inv[i]);
	}
}

inline int C(int n,int m){return (n>=0&&m>=0&&n>=m)?mul(fac[n],mul(ifac[m],ifac[n-m])):0;}

inline Poly operator*(Poly a,Poly b){
	int deg=a.size()+b.size()-1,l=1;
	while(l<deg)l<<=1;init_rev(l);
	a.resize(l);NTT(a,l,1);
	b.resize(l);NTT(b,l,1);
	for(int re i=0;i<l;++i)a[i]=mul(a[i],b[i]);
	NTT(a,l,-1);a.resize(deg);
	return a;
}

inline Poly Deriv(Poly a){
	for(int re i=0;i+1<a.size();++i)a[i]=mul(a[i+1],i+1);
	a.pop_back();return a;
}

inline Poly Integ(Poly a){
	a.push_back(0);
	for(int re i=a.size()-1;i;--i)a[i]=mul(a[i-1],inv[i]);
	a[0]=0;return a;
}

inline Poly Inv(cs Poly &a,int lim){
	Poly c,b(1,power(a[0],mod-2));
	for(int re l=4;(l>>2)<lim;l<<=1){
		init_rev(l);
		c=a,c.resize(l>>1);
		c.resize(l),NTT(c,l,1);
		b.resize(l),NTT(b,l,1);
		for(int re i=0;i<l;++i)b[i]=mul(b[i],dec(2,mul(b[i],c[i])));
		NTT(b,l,-1);b.resize(l>>1);
	}b.resize(lim);
	return b;
}

inline Poly Ln(Poly a,int lim){
	a=Integ(Deriv(a)*Inv(a,lim));
	a.resize(lim);
	return a;
}

inline Poly Exp(cs Poly &a){
	Poly c,b(1,1);int n=a.size();
	for(int re i=2;(i>>1)<n;i<<=1){
		c=Ln(b,i);
		for(int re j=0;j<i;++j)c[j]=dec(j<n?a[j]:0,c[j]);
		c[0]=add(c[0],1);
		b=b*c;b.resize(i);
	}b.resize(n);
	return b;
}

cs int N=1e5+5;

int n;

std::vector<int> factor;

int to[N],cnt[N];
bool vis[N];

Poly f,g;

signed main(){
	init_NTT();
	init_inv();
#ifdef zxyoi
	freopen("message.in","r",stdin);
#endif
	n=getint();
	for(int re i=1;i<=n;++i)to[i]=getint();
	for(int re i=1;i<=n;++i)if(!vis[i]){
		int len=0,x=i;
		while(!vis[x]){
			vis[x]=true;
			++len;x=to[x];
		}
		++cnt[len];
	}
	for(int re i=1;i*i<=n;++i)
	if(n%i==0){
		factor.push_back(i);
		if(i*i!=n)factor.push_back(n/i);
	}
	std::sort(factor.begin(),factor.end());
	int ans=1;
	for(int re i=1;i<=n&&ans;++i)if(cnt[i]){
		f=Poly(cnt[i]+1,0);
		for(int re j:factor){
			if(j>cnt[i])break;
			if(gcd(j*i,n)!=j)continue;
			f[j]=power(i,j-1,inv[j]);
		}
		g=Exp(f);
		ans=mul(ans,mul(g[cnt[i]],fac[cnt[i]]));
	}
	cout<<ans<<"\n";
	return 0;
}
zxyoi_dreamer
关注
  • 0
    点赞
  • 踩
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
多项式exp的正确姿势
zyy大蒟蒻的博客
04-16 1367
等到晚上应该可以更新当日题解。 D1T1 设DP状态fi,0/1,jf_{i,0/1,j}fi,0/1,j​表示前iii个数字，总共有jjj个选择了A,上一个选择了A/B是否合法，直接暴力复杂度O(n2)O(n^2)O(n2)。 观察性质，可以证明对于每个iii，其合法的肯定为一个连续段。证明略。 因此转移只需要记录左右端点。复杂度O(n)O(n)O(n) #include<bits/std...
luogu P4726 【模板】多项式指数函数（多项式 exp）
ACoder
04-14 253
形式幂级数求exp
多项式exp学习笔记
QQ82272760的博客
03-22 451
文章目录问题前置知识思路代码 问题 现在有一个多项式 A(x)A(x)A(x)，保证其常数项为 000，要求多项式 B(x)=eA(x)B(x)=e^{A(x)}B(x)=eA(x)。 前置知识 泰勒多项式 对于一个任何一个函数 f(x)f(x)f(x)，在点 x0x_0x0​ 附近都有f(x)=∑i=0n1i!f(i)(x0)(x−x0)i+Rn(x)f(x)=\sum_{i=0}^n\frac{1}{i!}f^{(i)}(x_0)(x-x_0)^i+R_n(x)f(x)=i=0∑n​i!1​f(i)(x
Exp的组合意义与Euler变换
00ffcc的博客
08-31 1187
ExpExpExp的组合意义与EulerEulerEuler变换 ExpExpExp的组合意义 根据泰勒展开,exp⁡(f(X))=∑i=0∞fi(X)i!\exp(f(X))=\sum_{i=0}^{\infty}\frac{f^i(X)}{i!}exp(f(X))=∑i=0∞​i!fi(X)​ 设F(X)=exp⁡(f(X))F(X)=\exp(f(X))F(X)=exp(f(X)) 注意f0=0f_0=0f0​=0 组合意义 fff表示集合内部的分配方案数关于集合大小的EGF,FFF表示有标号小球分
多项式全家桶——Part.4 多项式ln、exp、快速幂
RainbowCrown_CGH's blog
09-03 680
这拿头来学……
组合数学第五版课后习题答案，英文版，详解。
10-27
组合数学是数学的一个重要分支，主要研究有限集合中对象的排列、组合以及各种计数问题。这门学科在计算机科学中扮演着至关重要的角色，特别是在算法设计、概率论、图论、编码理论等领域都有广泛应用。这里提供的...
组合数学机械化通用程序库软件V10用户手册.docx
11-21
2. 对称函数实际：对称函数实际是代数组合学中的一个重要研讨范围，它主要研讨对称群和对称多项式的代数性质和组合性质，在数学的其他分支和数学物理中有宽广的运用，是一个遭到普遍关注的研讨方向。 3. 组合双射...
MANET 中基于二次置换多项式的滑动窗口网络编码 算法
02-26
针对移动自组网络（mobile Ad hoc network，MANET）的数据传输机制和数据吞吐量问题，提出了一种 MANET 中基于 二 次 置 换 多 项 式 的 滑 动 窗 口 网 络 编 码（quadratic permutation polynomials-based sliding...
组合数学 组合数 整数划分 递推 方程 多项式定理
08-03
讲解组合数学的知识，大连理工应用数学系讲义 对学习组数学的朋友有很大的帮助 第二章 多项式定理及其应用
信息安全数学基础-许春香-第五章 多项式环与有限域.ppt
08-08
信息安全数学基础-多项式环与有限域 本章节主要讨论了多项式环与有限域的概念和性质，在信息安全数学基础中扮演着重要的角色。 一、多项式环的定义 在域F上，一个一元多项式可以表示为f(x) = anxn + an-1xn-1 +…...
多项式指数函数（多项式 exp）
dygxczn的博客
05-06 564
已知gx，你要求出fx≡egxmodxn。ln。
多项式exp
gigo_64 sink in my world
08-08 392
关于多项式exp 求 则 则 设 使用牛顿迭代来求G的0点。 根据泰勒展开，，牛顿迭代就是只取i=0和1的两项，虽然不精确但可以逐渐缩小。 而多项式倍增的时候，从平方开始在%x^n的意义下，都是同余的，所以对于多项式，牛顿迭代得到的是准确值。 令b0为B的倍增下层的多项式。 如果在0点展开，则在递归倍增的时候， 所以 因为 所以 所以递归，然后每次多项式ln，再卷起来。 ...
【模板】多项式指数函数（多项式 exp）
m0_74821245的博客
12-11 431
int r[N];{ll b=1;{x=x*x%mod;t>>=1;}}{i
EGF中多项式exp的组合意义
最新发布
sophilex的博客
10-24 329
仔细理解一下，众所周知，一个大小为n的排列一定可以拆成若干个环，每一个环内部的排列数就是一个圆排列的方案数，所以大小为n的排列的方案数就是把。EGF一般用来处理多重集的排列问题，在其上可以定义多项式的exp运算，在处理一类问题的时候有独特的作用。分成若干个非空集合，每一个集合的圆排列方案数之积，这与我们上面讲到的exp的组合意义相符合。的个数，这是一个只跟集合大小有关的值），那么根据生成函数的定义，我们不难得到下式。是因为这k个集合是无序的，而原本的多个多项式卷积显然是有序的。考虑大小为n的排列的个数是。
一道多项式exp模板题
Mys_C_K的博客
07-28 465
多项式什么的真难写QwQ 题目大意：已知： f(x)=∏ni=1(1+aix)g(x)=∏mi=1(1+bix)h(x)=∏ni=1∏mj=1(1+aibjx)f(x)=∏i=1n(1+aix)g(x)=∏i=1m(1+bix)h(x)=∏i=1n∏j=1m(1+aibjx)f(x)=\prod_{i=1}^n(1+a_ix)\\g(x)=\prod_{i=1}^m(1+b_ix)\\h(x)...
【校内模拟】子树问题（组合数学DP）（多项式exp）
zxyoi_dreamer的博客（不定期诈尸）
11-02 296
简要题意： 请你对满足如下限制的树计数： 根节点深度为111，最大节点深度在L−RL-RL−R之间（分别回答） 节点的标号满足堆性质 给定正整数集合{a}\{a\}{a}，没有任何一个节点的子树大小在集合中。 n≤500n\leq 500n≤500 题解： 考虑子树拼接可以得到一个比较好懂的O(n3)O(n^3)O(n3)组合数学DP。 然而我对这个模型实在是太熟悉了，于是考场上想都没想直接去...
多项式全家桶（三）：多项式的ln，exp
周道-Althen的博客
04-03 1012
      &ThinSpace;\ \ \ \ \ \ \,      对数和指数是很重要的东西了，复杂的多项式都和他们有关系，所以说掌握很重要，这里不建议光背板子，因为这两个板子都有致命的限制，而在实际操作的时候，可以通过一些方法绕过这个限制直接求解，这个就很重要了。...
[带标号DAG计数 容斥原理 多项式求逆 多项式求ln] COGS 有标号的DAG计数系列
雯舞
02-24 1050
题目链接 传送门I 传送门II 传送门III 传送门IV 题解戳这里 大致解法和思想是同这里一样的DAG1找出一个入度为00的点 那么答案为C1n∗21∗(n−1)∗fn−1C_n^1*2^{1*(n-1)}*f_{n-1} 但是一张DAG里会有多个入度为00的点会重复计数 考虑容斥 得到 fn=∑k=1n−1k−1∗Ckn∗2k∗(n−k)∗fn−kf_n=\sum_{k=1}^n -1
Matlab实验一传递函数表示方法.pptx
11-12
在Matlab中，传递函数是控制系统理论中常用的一种数学模型，用于描述系统输出与输入之间的关系。本实验主要围绕以下几个知识点展开： 1. **传递函数的因子形式与多项式形式**： - 因子形式传递函数通常表示为直接...

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值