NC53079 Forsaken喜欢数论 线性筛 每个数字的最小质因子和

链接:https://ac.nowcoder.com/acm/problem/53079
来源:牛客网

题目描述
Forsaken有一个有趣的数论函数。对于任意一个数xxx,f(x)f(x)f(x)会返回xxx的最小质因子。如果这个数没有最小质因子,那么就返回0。
现在给定任意一个nnn,Forsaken想知道∑i=1nf(i)\sum_{i = 1}^{n}{f(i)}∑i=1n​f(i)的值。

输入描述:

一个整数nnn。

输出描述:

一个整数代表上面的求和式的值。

示例1
输入
复制

4

输出
复制

7

备注:

1≤n≤3e71 \leq n \leq 3e71≤n≤3e7


线性筛塞一下就好了

#define int long long
int n, m, Q, K;

int prime[MAXN], pn = 0, tmin[MAXN], N;
ll ans = 0;
bool vis[MAXN];

void get_prime() {
	vis[0] = vis[1] = true;
	for(int i=2; i<n; i++) {
		if(!vis[i]) { 
			prime[pn++] = i;
			tmin[i] = i;
		}
		for(int j=0; j<pn && prime[j]*i<n; j++) {
			int p = prime[j]*i;
			vis[p] = true;
			tmin[p] = prime[j];
			// if(!tmin[p] && !vis[i]) tmin[p] = i;
			if(i%prime[j] == 0) break;
		}
		ans += tmin[i];
		// show(tmin[i]);
	}
}

signed main() {
#ifdef debug
	freopen("test.txt", "r", stdin);
	clock_t stime = clock();
#endif
	read(n);
	n ++;
	get_prime();
	printf("%lld\n", ans);
	// forarr(prime, 0, 10);




#ifdef debug
	clock_t etime = clock();
	printf("rum time: %lf 秒\n",(double) (etime-stime)/CLOCKS_PER_SEC);
#endif 
	return 0;
}

根据提供的引用内容,可以看出Forsaken喜欢数论,并且在数论方面有一定的研究。以下是一些与数论相关的Forsaken的例子: 1. 线性筛法求质数 ```python def get_primes(n): is_prime = [True] * (n + 1) primes = [] for i in range(2, n + 1): if is_prime[i]: primes.append(i) for j in range(len(primes)): if i * primes[j] > n: break is_prime[i * primes[j]] = False if i % primes[j] == 0: break return primes ``` 2. 分组算法 ```python def group_students(students, k): groups = [[] for _ in range(k)] for i in range(len(students)): groups[i % k].append(students[i]) return groups ``` 3. 正方形计数 ```python def count_squares(n): count = 0 for i in range(1, int(n ** 0.5) + 1): count += (n // i - i + 1) return count ``` 4. 三维数点 ```python def count_3d_points(points): bit = [0] * (max(points, key=lambda x: max(x)) + 1) count = 0 for i in range(len(points)): x, y, z = points[i] count += i - sum(bit[j] for j in range(x)) for j in range(x, len(bit)): bit[j] += 1 return count ``` 5. 独一无二的最小生成树 ```python def unique_minimum_spanning_tree(n, edges): edges.sort(key=lambda x: x[2]) parent = list(range(n)) rank = [0] * n def find(x): if parent[x] != x: parent[x] = find(parent[x]) return parent[x] def union(x, y): px, py = find(x), find(y) if px == py: return False if rank[px] < rank[py]: parent[px] = py elif rank[px] > rank[py]: parent[py] = px else: parent[px] = py rank[py] += 1 return True mst = [] for u, v, w in edges: if union(u, v): mst.append((u, v, w)) return mst --相关问题--:
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

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

余额充值