UOJ62 怎样跑的更快

关注

原创最新推荐文章于 2026-03-16 07:45:00 发布 573 阅读 AI 写同款· GEO 优化 ›

在Linux机器上部署UOJ评测系统前言部署要求开始部署 Ubuntu软件源更新 docker服务部署获取UOJ最新社区版安装UOJ 继续安装 UOJ相关操作部署GUI管理系统更多内容前言学校里想要部署一套OJ，然后就选了UOJ。但是由于服务器还要一段时间（永远）才能到，所以就先拿了一台电脑作为试验机进行内网部署测试。非常感谢以下人： vfleaking 无私开放了UOJ的源码阮行止维护社区版本UOJ源阅读详情

$bi=∑j=1n(i,j)c−d∗idjdzjb_i=\sum_{j=1}^n(i,j)^{c-d}*i^dj^dz_j$
$biid=∑j=1nzjjd(i,j)c−d\frac{b_i}{i^d}=\sum_{j=1}^nz_jj^d(i,j)^{c-d}$ 记 $(i,j)=x,Bi=biid,Zi=ziid,fi(d)=∑(i,j)=dZj,Fi(d)=∑d∣(i,j)Zj(i,j)=x,B_i=\frac{b_i}{i^d},Z_i=z_ii^d,f_i(d)=\sum_{(i,j)=d}Z_j,F_i(d)=\sum_{d|(i,j)}Z_j$
$biid=∑x∣inf(x)xc−d\frac{b_i}{i^d}=\sum_{x|i}^nf(x)x^{c-d}$ $biid=∑x∣inxc−d∑kF(kx)μ(k)\frac{b_i}{i^d}=\sum_{x|i}^nx^{c-d}\sum_{k}F(kx)\mu(k)$ $biid=∑kx∣inxc−dF(kx)μ(k)\frac{b_i}{i^d}=\sum_{kx|i}^nx^{c-d}F(kx)\mu(k)$ $B(i)=∑T∣inF(T)∑x∣Tμ(T/x)xc−dB(i)=\sum_{T|i}^nF(T)\sum_{x|T}\mu(T/x)x^{c-d}$ $令g(T)=μ∗idc−d(T)令g(T)=\mu*id^{c-d}(T)$ 得到 $B(i)=∑T∣iF(T)g(T)B(i)=\sum_{T|i}F(T)g(T)$ 那么 $F(T)g(T)=∑ij=Tμ(j)B(i)F(T)g(T)=\sum_{ij=T}\mu(j)B(i)$ 让后一路反推回去得到 $z$ 就行了

#include<stdio.h>
#include<string.h>
#include<algorithm>
#define N 100010
#define M 998244353
#define LL long long
using namespace std;
inline LL pow(LL x,LL k,int& s){
	while(k<0) k+=M-1;
	for(s=1;k;x=x*x%M,k>>=1) k&1?s=s*x%M:0;
}
int mu[N],c,d,n,q,inv[N];
int f[N],g[N],b[N],B[N],z[N],Z[N],F[N],xd[N],xcd[N],G[N];
int main(){
	scanf("%d%d%d%d",&n,&c,&d,&q);
	mu[1]=1;
	for(int i=1;i<=n;++i){
		pow(i,d,xd[i]);
		pow(i,c-d,xcd[i]);
		pow(xd[i],M-2,inv[i]);
		for(int j=i+i;j<=n;j+=i) mu[j]-=mu[i];
	}
	for(int i=1;i<=n;++i) if(mu[i])
		for(int j=1,k=i;k<=n;++j,k+=i)
			g[k]=(g[k]+mu[i]*xcd[j])%M;
	while(q--){
		for(int i=1;i<=n;++i){
			scanf("%d",b+i);
			B[i]=(LL)b[i]*inv[i]%M;
		}
		memset(G,0,n+1<<2);
		for(int i=1;i<=n;++i) if(mu[i])
			for(int j=1,k=i;k<=n;++j,k+=i)
				G[k]=(G[k]+mu[i]*B[j])%M;
		for(int i=1;i<=n;++i)
			if(!g[i]&&G[i]){
				puts("-1");
				goto end;
			} else{
				pow(g[i],M-2,F[i]); 
				F[i]=(LL)F[i]*G[i]%M;
			}
		memset(Z,0,n+1<<2);
		for(int i=1;i<=n;++i) if(mu[i])
			for(int j=1,k=i;k<=n;++j,k+=i)
				Z[j]=(Z[j]+mu[i]*F[k]%M)%M;
		for(int i=1;i<=n;++i){
			z[i]=(LL)Z[i]*inv[i]%M;
			printf("%d ",(z[i]+M)%M);
		}
		puts(""); end:;
	}
}