Loading

牛客多校第四场 H.Convolution 数论 推式子

牛客多校第四场 H.Convolution 数论 推式子

题意

定义运算

\[a \bigotimes b = \sum_{i}p_i^{|e_{a_i} - e_{b_i}|} \]

\(e_{x_i}\)表示\(x\)的第\(i\)个质因子的指数

给定\(a_1...a_n\)

\[b_i = \sum_{i \leq k,j \leq n,j \bigotimes k =i} a_jk^c \]

分析

现进行一些转化

有常用的等式

\[|a-b|= max(a,b) - min(a,b) \]

那么

\[a\bigotimes b = \sum_{i}p_i^{max(e_{a_i},e_{b_{i}}) - min(e_{a_i},e_{b_i})} \]

显然这就是\(\frac{[a,b]}{(a,b)} = \frac{ab}{(a,b)^2}\)

那么

\[b_i = \sum_{j=1}^n\sum_{k=1}^n[\frac{jk}{(j,k)^2} = i]a_j k^c\\ \]

考虑对称性构造\(x' = \frac{j}{(j,k)},y' = \frac{k}{(j,k)}\)

那么常用套路,枚举\(gcd = d\),\(j = dx',k = dy'\)

\[b_i = \sum_{x' | i} \sum_{d=1}^{min(\lfloor\frac{i}{x'}\rfloor,\lfloor\frac{i}{y'}\rfloor)}a_{dx'}(dy')^c \]

发现每次计算\(b_i\)中后式会多次重复计算

考虑计算

\[f_{x',m} = \sum_{d=1}^m a_{dx'}d^c \]

显然这可以\(O(1)\)递推,这样式子就只依赖\(a\) 计算贡献时可以额外乘上\(y\)的贡献

\[b_i = \sum_{x|i} f_{x',u} y'^c \]

代码

const int MOD = 998244353;

inline int mul(int a,int b){
	int res = (ll)a * b % MOD;
	if(res < 0) res += MOD;
	return res;
}

inline void add(int &a,int b){
	a += b;
	if(a >= MOD) a -= MOD;
}

inline void sub(int &a,int b){
	a -= b;
	if(a < 0) a += MOD;
}

const int maxn = 1e6 + 5;

int ksm(int a,int b = MOD - 2,int m = MOD){
	int ans = 1;
	int base = a;
	while(b){
		if(b & 1) ans = (ll)ans * base % m;
		base = (ll)base * ans % m;
		b >>= 1;
	}
	return ans;
}


int p[maxn];

int main(){
	int n = rd();
	int c = rd();
	for(int i = 1;i <= n;i++)
		p[i] = ksm(i,c);
	vector<int> v(n + 1);
	vector<int> ans(n + 1);
	for(int i = 1;i <= n;i++){
		v[i] = rd();
	}
	int res = 0;
	for(int i = 1;i <= n;i++){
		vector<int> f(n / i + 1);
		for(int j = 1;j * i <= n;j++)
			f[j] = f[j - 1],add(f[j],mul(v[i * j],p[j]));
		for(int j =1;i * j <= n;j++)
			if(__gcd(i,j) == 1){
				int u = n / max(i,j);
				add(ans[i * j],mul(f[u],p[j]));		
			}
	}
	for(int i = 1;i <= n;i++)
		res ^= ans[i];
	printf("%d",res);
}
posted @ 2021-08-06 18:38  MQFLLY  阅读(77)  评论(0编辑  收藏  举报