Send a Table(歐拉函數埃式篩)

UVA10375

在這裏插入圖片描述
顯然,題目要求有多少二元組(x,y)(x,y)滿足1x,ynx,y1\leq x,y\leq n且x,y互素
根據歐拉函數ϕ\phi的定義,ϕ(n)\phi (n)爲小於n且與n互素的整數個數。
定義f(n)=i=2nϕ(i)f(n)=\sum\limits_{i=2}^{n}\phi(i)
那麼對於x<yx>yx< y或x>y,都有f(n)f(n)個整數對滿足條件,再加上(1,1)(1,1),答案爲2f(n)+12*f(n)+1

#include <iostream>
#include <algorithm>
#include <cmath>
#include <vector>
using namespace std;

const int MAXN = 50001;

int phi[MAXN + 1];
int sum[MAXN + 1];

//歐拉函數和埃式篩融合
void Initphi() {
	phi[1] = 1;
	for (int i = 2; i <= MAXN; ++i) {
		if (!phi[i]) {
			for (int j = i; j <= MAXN; j += i) {
				if (!phi[j]) {
					phi[j] = j;
				}
				//能進來到這的i都是素數
				phi[j] = phi[j] / i * (i - 1);
			}
		}
	}
	return;
}

void InitSum() {
	sum[2] = phi[2];
	for (int i = 3; i <= MAXN; ++i) {
		sum[i] += sum[i - 1] + phi[i];
	}
}

int main(){
	Initphi();
	InitSum();
	int n;
	while (cin >> n && n) {
		printf("%d\n", 2 * sum[n] + 1);
	}
	return 0;

}
發表評論
所有評論
還沒有人評論,想成為第一個評論的人麼? 請在上方評論欄輸入並且點擊發布.
相關文章