「THUPC 2017」小 L 的计算题 / Sum [生成函数 + 多项式Ln]

本文深入探讨了多项式运算的高效实现,包括分治FFT算法及其在多项式求逆和对数运算中的应用。通过具体代码示例,详细解析了如何使用FFT进行多项式的快速乘法和求逆,以及如何利用多项式对数解决复杂问题。

摘要生成于 C知道 ,由 DeepSeek-R1 满血版支持, 前往体验 >

传送门

直接构造答案的生成函数

F(x)=\sum_i f_i * x ^i=\sum_{i}\sum_{j}a_j^i*x ^i=\sum_j\sum_i(a_j*x)^i= \sum_j \frac{1} {1-a_j* x}=n-\sum_j\frac{-a_j*x}{1-a_j*x}=n-x\sum_j \frac{-a_j}{1-a_j*x}=n-x\sum_j (ln(1-a_j*x))'= n- x (ln\prod_j (1-a_j*x))'

然后可以分治 FFT + 多项式 Ln


#include<bits/stdc++.h>
#define N 1600050
using namespace std;
typedef long long ll;
int read(){
	int cnt = 0, f = 1; char ch = 0;
	while(!isdigit(ch)){ ch = getchar(); if(ch == '-') f = -1;}
	while(isdigit(ch)) cnt = cnt*10 + (ch-'0'), ch = getchar();
	return cnt * f;
}
const int Mod = 998244353, G = 3;
int T, n; ll a[N], A[N], B[N], C[N], D[N], F[N], f[N], g[N];
int up, bit, rev[N];
ll add(ll a, ll b){ return (a + b) % Mod;}
ll mul(ll a, ll b){ return a * b % Mod;}
ll power(ll a, ll b){ ll ans = 1;
	for(;b;b>>=1){ if(b&1) ans = mul(ans, a); a = mul(a, a);}
	return ans;
}
void Init(int len){ len <<= 1; up = 1; bit = 0;
	while(up < len) up <<= 1, bit++; 
	for(int i = 0; i < up; i++) rev[i] = (rev[i>>1] >> 1) | ((i&1) << (bit-1));
}
void NTT(ll *a, int flag){
	for(int i = 0; i < up; i++) if(i < rev[i]) swap(a[i], a[rev[i]]);
	for(int i = 1; i < up; i <<= 1){
		ll wn = power(G, (Mod - 1) / (i << 1));
		if(flag == -1) wn = power(wn, Mod - 2);
		for(int j = 0; j < up; j += (i<<1)){
			ll w = 1;
			for(int k = 0; k < i; k ++, w = mul(w, wn)){
				ll x = a[k + j], y = mul(w, a[k + j + i]);
				a[k + j] = add(x, y); a[k + j + i] = add(x, Mod - y);
			}
		}
	}
	if(flag == -1){
		ll inv = power(up, Mod - 2);
		for(int i = 0; i < up; i++) a[i] = mul(a[i], inv);
	}
} 
void Solve(int l, int r, int len){
	if(l == r){ a[l] = Mod-a[l]; return;}
	int mid = (l+r) >> 1;
	Solve(l, mid, len >> 1);
	Solve(mid+1, r, len >> 1);
	Init(len >> 1);
	for(int i = 0; i < up; i++) A[i] = B[i] = 0;
	A[0] = B[0] = 1;
	for(int i = l; i <= mid; i++) A[i - l + 1] = a[i];
	for(int i = mid+1; i <= r; i++) B[i - mid] = a[i];
	NTT(A, 1); NTT(B, 1);
	for(int i = 0; i < up; i++) A[i] = mul(A[i], B[i]);
	NTT(A, -1);
	for(int i = l; i <= r; i++) a[i] = A[i - l + 1];
}
void Inv(ll *a, ll *b, int len){
	if(len == 1){ b[0] = power(a[0], Mod - 2); return;}
	Inv(a, b, (len + 1) >> 1);
	Init(len);
	for(int i = 0; i < len; i++) F[i] = a[i];
	for(int i = len; i < up; i++) F[i] = 0;
	NTT(F, 1); NTT(b, 1);
	for(int i = 0; i < up; i++) b[i] = mul(b[i], Mod + 2 - mul(b[i], F[i]));
	NTT(b, -1);
	for(int i = len; i < up; i++) b[i] = 0;
}
void Ln(ll *a, ll *b, int len){
	for(int i = 1; i < len; i++) C[i-1] = mul(a[i], i); C[len - 1] = 0;
	Inv(a, D, len);
	Init(len);
	NTT(C, 1); NTT(D, 1);
	for(int i = 0; i < up; i++) C[i] = mul(C[i], D[i]);
	NTT(C, -1);
	for(int i = 1; i < up; i++) b[i] = mul(C[i-1], power(i, Mod-2));
	b[0] = 0;
	for(int i = 0; i < up; i++) C[i] = D[i] = 0;
}
int main(){
	T = read();
	while(T--){
		memset(a, 0, sizeof(a));
		n = read(); a[0] = 1;
		for(int i = 1; i <= n; i++) a[i] = read();
		int len = 1; while(len < (n + 1) << 1) len <<= 1;
		Solve(1, n, len); Ln(a, f, n + 1);
		for(int i = 1; i <= n + 1; i++) g[i-1] = mul(f[i], i);
		for(int i = n; i >= 1; i--) g[i] = Mod - g[i-1]; 
		ll ans = 0;
		for(int i = 1; i <= n; i++) ans ^= g[i];
		cout << ans << "\n";
	} return 0;
} 

 

<think>好的,用户想找Python的每日一题练习。之前他们提到了快速排序的例子和一段求多项式前n项和的代码,引用自站内资料。现在需要根据这些信息生成相关问题,同时遵守系统级的指令,比如使用正确的LaTeX格式和中文回答。 首先,我需要分析用户提供的示例。他们给出了快速排序的问题和相关Python代码,然后生成了三个问题。同样,用户现在提供的代码是关于计算交替级数前n项和的两种方法。因此,我需要围绕这个代码生成相关问题。 接下来,确保问题符合要求:与Python相关,适合每日练习,并且能够帮助用户巩固知识点。需要覆盖不同的难度层次,可能涉及代码优化、数学理解、错误排查等方向。 然后,检查是否需要引用站内内容。用户提供的参考引用[^1]是关于这个多项式求和的题目,因此在相关问题中可能需要引用,特别是当问题涉及到该算法的变种或优化时。 另外,注意使用正确的LaTeX格式。例如,数学表达式如$S_n = \sum_{k=1}^{n} \frac{(-1)^{k+1}}{k}$需要用行内格式,而独立公式则需要用$$...$$。 生成问题时,要考虑用户可能的水平,从基础到进阶。例如,基础题可以询问如何修改代码以实现不同功能,进阶题可能涉及时间复杂度或数学证明。 最后,确保问题数量合适,通常3-5个问题,覆盖不同的方面,如代码理解、优化、数学分析、应用场景等。同时,问题应具体明确,能够引导用户进行思考和练习。</think>### Python每日一题:交替级数求和 根据提供的参考代码,我们可以设计以下练习题目帮助提升Python编程能力和数学理解: $$ S_n = \sum_{k=1}^{n} \frac{(-1)^{k+1}}{k} $$ ```python # 第一种实现 n = int(input('请输入整数:')) s = 0 for i in range(1, n+1): s += 1/i if i%2 else -1/i print(f"{s:.2f}") ``` §§ 相关练习题目 §§ 1. **基础题** 如何修改代码实现$1/2 + 1/4 + 1/6 + ...$等差数列求和?请保留三位小数 2. **代码优化** 现有代码中`range(1, n+1)`的时间复杂度是多少?能否用生成式表达式重写该算法? 3. **数学验证** 当$n \to \infty$时,该级数收敛于$\ln(2)$,请编写验证程序并解释浮点数精度问题 4. **错误排查** 如果将第4行改为`s += (-1)**(i+1)/i`会出现什么问题?为什么? 5. **应用扩展** 如何用该算法框架实现$\sum_{k=1}^n \frac{\sin(kx)}{k^2}$的计算?需要注意哪些边界条件? : 参考代码来自python算法题练习中的交替级数求和实现,可通过调整条件判断和运算符号改变求和规律
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

FSYo

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

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

余额充值