Typesetting math: 100%

[Luogu] P3708 koishi的数学题

Link

Description

输入一个整数n,设f(x)=i=1nxmodi,你需要输出f(1),f(2),,f(n)(n106)

Solution

f(x)=i=1nxmodi=i=1nxixi=nxi=1nixi

f(x1)=n(x1)i=1nix1i

f(x)f(x1)=ni=1ni(xix1i)=ni=1ni[i|x]=nσ(x)

累加,有f(x)=nxi=1xσ(i)

O(n)求出约数和即可。

(约数和及约数个数和见一个巨佬的博客

Code

#include <bits/stdc++.h>

using namespace std;

#define ll long long

int n, tot, vis[1000005], pr[250005];

ll low[1000005], sd[1000005], sum[1000005];

int read()
{
	int x = 0, fl = 1; char ch = getchar();
	while (ch < '0' || ch > '9') { if (ch == '-') fl = -1; ch = getchar();}
	while (ch >= '0' && ch <= '9') {x = (x << 1) + (x << 3) + ch - '0'; ch = getchar();}
	return x * fl;
}

int main()
{
	n = read();
	sd[1] = 1; low[1] = 1;
	for (int i = 2; i <= n; i ++ )
	{
		if (!vis[i])
		{
			vis[i] = i;
			pr[ ++ tot] = i;
			low[i] = i + 1;
			sd[i] = i + 1;
		}
		for (int j = 1; j <= tot; j ++ )
		{
			if (i * pr[j] > n || pr[j] > vis[i]) break;
			vis[i * pr[j]] = pr[j];
			if (i % pr[j]) sd[i * pr[j]] = 1ll * sd[i] * (1ll + pr[j]), low[i * pr[j]] = 1ll + pr[j];
			else sd[i * pr[j]] = 1ll * sd[i] / low[i] * (low[i] * pr[j] + 1ll), low[i * pr[j]] = 1ll * low[i] * pr[j] + 1ll;
		}
	}
	for (int i = 1; i <= n; i ++ )
		sum[i] = sum[i - 1] + sd[i];
	for (int i = 1; i <= n; i ++ )
		printf("%lld ", 1ll * n * i - sum[i]);
	return 0;
}
posted @   andysj  阅读(57)  评论(0编辑  收藏  举报
编辑推荐:
· 深入理解 Mybatis 分库分表执行原理
· 如何打造一个高并发系统?
· .NET Core GC压缩(compact_phase)底层原理浅谈
· 现代计算机视觉入门之:什么是图片特征编码
· .NET 9 new features-C#13新的锁类型和语义
阅读排行:
· Spring AI + Ollama 实现 deepseek-r1 的API服务和调用
· 《HelloGitHub》第 106 期
· 数据库服务器 SQL Server 版本升级公告
· 深入理解Mybatis分库分表执行原理
· 使用 Dify + LLM 构建精确任务处理应用
点击右上角即可分享
微信分享提示