打卡信奥刷题（3515）用C++实现信奥题 P10886 【MX-S3-T2】「FeOI Round 1」Journey

原创于 2026-08-18 08:51:59 发布 · 366 阅读

3 ·

本内容遵循CC 4.0 BY-SA版权协议

GEO检测

收录于

C++

P10886 【MX-S3-T2】「FeOI Round 1」Journey

题目背景

原题链接：https://oier.team/problems/S3B。

外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传

题目描述

小 W 最近正在学习某编程语言。

这个编程语言有个语句如下：

range(a,b,c)

这个语句表示一个序列，这个序列是这样的：

$[a,a+c,a+2c,\cdots,a+kc]$

其中， $k$ 是最大的满足 $a + k c < b$ 的非负整数 $k$ 。例如 range(1,7,2) 表示的序列为 $[1, 3, 5]$ 。

小 W 想问你一个问题：给你一个长为 $n$ 的序列 $g$ （下标从 $1$ 开始），求出下列式子的值，答案模 $10^9 +7$ 。

$\sum_{a=1}^{n}\sum_{b=a+1}^{n+1}\sum_{c=1}^{n}\sum\limits_{i\in \mathrm{range}(a,b,c)} g_i$

输入格式

为了加快读入速度，我们采用以下读入方法。

读入共一行五个非负整数 $n,A,B,C,g_n$ 。

其中 $n$ 为序列 $g$ 的长度， $A, B, C$ 为生成数据的参数， $g_n$ 为序列 $g$ 的第 $n$ 项。

对于序列 $g$ 的第 $i$ （ $1\leq i < n$ ）项，满足 $g_i \equiv Ag_{i+1}^2+Bg_{i+1}+C \pmod {10^9 + 7}$ 。

输出格式

共一行一个非负整数，为题目中所求的答案。

输入输出样例 #1

输入 #1

2 0 1 1 1

输出 #1

输入输出样例 #2

输入 #2

9 0 1 0 663

输出 #2

输入输出样例 #3

输入 #3

20 1 0 0 998244353

输出 #3

560706529

输入输出样例 #4

输入 #4

114514 17723 134 1045 233337

输出 #4

442762986

说明/提示

【样例解释 #1】

$g = [2, 1]$ （下标从 $1$ 开始）

令 $ans=\sum\limits_{i\in \mathrm{range}(a,b,c)} g_i$ 。

当 $a = 1, b = 2, c = 1$ 时， $an s = 2$ 。

当 $a = 1, b = 2, c = 2$ 时， $an s = 2$ 。

当 $a = 1, b = 3, c = 1$ 时， $an s = 2 + 1 = 3$ 。

当 $a = 1, b = 3, c = 2$ 时， $an s = 2$ 。

当 $a = 2, b = 3, c = 1$ 时， $an s = 1$ 。

当 $a = 2, b = 3, c = 2$ 时， $an s = 1$ 。

答案为 $\sum ans=2+2+3+2+1+1=11$ 。

【样例解释 #2】

该样例满足特殊性质 A。

【数据范围】

对于 $100\%$ 的数据： $1\leq n\leq 2\times 10^7$ ， $0\leq A,B,C,g_n<10^9 +7$ 。

测试点编号	$n =$	特殊性质
$1$	$1$	AB
$2$	$1000$	无
$3$	$1000$	无
$4$	$5000$	无
$5$	$5000$	无
$6$	$10^4$	A
$7$	$10^5$	A
$8$	$10^5$	B
$9$	$10^5$	无
$10$	$2\times 10^5$	A
$11$	$2\times 10^5$	B
$12$	$5\times 10^5$	无
$13$	$10^6$	A
$14$	$10^6$	B
$15$	$2\times 10^6$	无
$16$	$3\times 10^6$	无
$17$	$5\times 10^6$	A
$18$	$10^7$	A
$19$	$10^7$	B
$20$	$1.3\times 10^7$	无
$21$	$1.6\times 10^7$	无
$22$	$2\times 10^7$	A
$23$	$2\times 10^7$	B
$24$	$2\times 10^7$	无
$25$	$2\times 10^7$	无

特殊性质 A：保证所有 $g_i$ 相等。

特殊性质 B：保证只有 $g_n\neq 0$ 。

C++实现

#include<bits/stdc++.h>
using namespace std;
const int N=2e7+10,mod=1e9+7;
int n,vis[N],pr[N],m,t,a[N],f[N],b[N];
inline void get(int n){
	f[1]=1;
	for (int i=2;i<=n;i++){
		if (!vis[i]) pr[++m]=i,a[i]=1,f[i]=2;
		for (int j=1;i*pr[j]<=n;j++){
			int x=i*pr[j];
			vis[x]=1;
			if (i%pr[j]==0){
				a[x]=a[i]+1,f[x]=f[i]/a[x]*(a[x]+1);
				break;
			} 
			else a[x]=1,f[x]=f[i]*2;
		}
	}
	for (int i=1;i<=n;i++) f[i]+=f[i-1];
}
inline int read(){
	char c=getchar();int f=1;int ans=0;
	while(c<48||c>57) (c==45?f=-1:1),c=getchar();
	while(c>=48&&c<=57) ans=(ans<<1)+(ans<<3)+(c^48),c=getchar();
	return ans*f;
}
main(){
	int n=read(),A=read(),B=read(),C=read();b[n]=read();
	get(n);
	for (int i=n-1;i>0;i--)
		b[i]=(1ll*A*b[i+1]%mod*b[i+1]%mod+1ll*B*b[i+1]%mod+C)%mod;
	int ans=0;
	for (int i=1;i<=n;i++){
		int sum=0;
		sum=f[i-1];
		ans=(1ll*b[i]*(n+1-i)%mod*sum%mod+ans+1ll*b[i]*n%mod*(n+1-i)%mod)%mod;
	}
	cout <<ans;
	return 0;
}