DAY 9

今日得分：5+35+0（忘删调试了。。。）（好在比前两天加起来的分高）

T1

题目大意：定义排列的乘法如下：对于h=f*g,h[i]=g[f[i]](1<=i<=n)。求长度为n，且满足P^(m+1)=P的排列个数，要求m为满足条件的最小正整数。n<=15000,m<=1e9。

题解：不难发现排列构成若干个置换环，满足要求的排列中，所有置换环长度的lcm等于m。考虑莫比乌斯反演（看到lcm就可以往这方面想了），问题变为求x为置换环长度的lcm的约数的排列个数。发现有贡献的x很少，处理出来之后暴力跑dp即可。注意极其卡常。

AC代码：

#include<iostream>
#include<cstdio>
#include<cstring>
#include<string>
#include<algorithm>
#include<queue>
#include<vector>
#include<cmath>
using namespace std;
inline int re_ad()
{
	int x=0,f=1;char ch=getchar();
	while(ch>'9'||ch<'0'){if(ch=='-')f=-1;ch=getchar();}
	while(ch>='0'&&ch<='9')x=x*10+(ch^48),ch=getchar();
	return x*f;
}
const int mo=998244353;
const long long mo2=1ll*mo*1000000000ll;
int n,m;
inline int ksm(int x,int y){int ret=1;while(y){if(y&1)ret=1ll*ret*x%mo;y>>=1;x=1ll*x*x%mo;}return ret;}
int pre[15010],ni[15010];
inline int C(int x,int y){if(x<y)return 0;return 1ll*pre[x]*ni[y]%mo*ni[x-y]%mo;}
inline int jia(int &x,int y){x+=y;if(x>=mo)x-=mo;}
vector<int> s,e[2010];
int d[50],a[50],cnt;
inline int F(int x,int y){return 1ll*pre[x]*ni[x-y]%mo;}
inline long long kc(long long x){return x>=mo2?x-mo2:x;}
long long ans,gx[15010][510],dp[15010];
void dfs(int x,int now)
{
	if(x>cnt){if(now>n)s.push_back(now);return;}
	dfs(x+1,now);
	for(int i=1;i<=a[x];++i)now*=d[x],dfs(x+1,now);
}
inline void work(register int x)
{
	register int i;
	for(i=2;i*i<=x;++i)
	{
	if(!(