数论:唯一分解定理丶调和级数

H - Pairs Forming LCM
今天又偷学了一手唯一分解定理的小知识:
n=p${}_1$^e1 *p${}_2$^e2 *p${}_3$^e3 *…*p${}_k$^k.
设a,b为n的两个因子，则:
a=p${}_1$^a1 *p${}_2$^a2 *p${}_3$^a3 *…*p${}_k$^k.
b=p${}_1$^b1 *p${}_2$^b2 *p${}_3$^b3 *…*p${}_k$^k.
推论:
gcd(a,b)=p${}_1$^min(a1,b1) *p${}_2$^min(a2,b2) *p${}_3$^min(a3,b3) *…*p${}_k$^k.
lcm(a,b)=p${}_1$^max(a1,b1) *p${}_2$^max(a2,b2) *p${}_3$^max(a3,b3) *…*p${}_k$^k.
题目:

题解:
由上面可以知道:当lcm(a,b)==n时,max(a1,b1)==e1,max(a2,b2)==e2…max(ak,bk)==ek.
那么当ai == ei时,bi可以是[0,ei]，有ei+1种可能,bi == ei时,ai可以是[0,ei],有ei+1种可能,加起来就是2ei+2种可能.但ai == bi == ei时就重复了.所以是2ei+1种可能.但题目是a<=b.所以就是(2ei+1)/2+1种可能.
AC代码:

#include<iostream>
#include<cmath>
using namespace std;
#define ll long long
const ll N = 1e7+10;
bool v[N];
ll cnt,p[N];

void primes(ll n){
	for(ll i=2;i<=n;i++){
		if(v[i])continue;
		p[++cnt]=i;
		for(ll j=i;j<=n/i;j++)v[i*j]=true;
	}
}

int main(){
	primes(N); 
	int T;
	scanf("%d",&T);
	for(int tt=1;tt<=T;tt++){
		ll n,sum=1;
		scanf("%lld",&n);
		for(int i=1;i<=cnt&&p[i]*p[i]<=n;i++){
			if(n%p[i]==0){
				ll ccnt=0;
				while(n%p[i]==0)n/=p[i],ccnt++;
				sum*=(2*ccnt+1);
			}
		}
		if(n>1)sum*=3;
		printf("Case %d: %lld\n",tt,(sum+1)/2);
	}
	return 0;
} 

调和级数:
当n足够大时,sum(n)=ln(n)+C+1/(2n).
C为欧拉常数为0.57721566490153286060651209

AC代码:

#include<iostream>
#include<cmath>
using namespace std;
#define ll long long
const ll N = 1e4+5;
double C = 0.57721566490153286060651209;
double p[N];
ll num;

void pre(){
	num=1;
	for(ll i=1;i<=N;i++,num++){
		p[num]=p[num-1]+1.0/(1.0*i);
	}
}

int main(){
	pre(); 
	int T;
	scanf("%d",&T);
	for(int tt=1;tt<=T;tt++){
		ll n,sum=1;
		scanf("%lld",&n);
		if(n<=10000)
		printf("Case %d: %.10lf\n",tt,p[n]);
		else{
			double ans=log(n)+C+1.0/(2*n);
			printf("Case %d: %.10lf\n",tt,ans);
		}
	}
	return 0;
}