法雷数列-欧拉函数

一、法雷数列

法雷数列：对任意给定的一个自然数n，将分母小于等于n的不可约的真分数按升序排列，并且在第一个分数之前加上0/1，在最后一个分数之后加上1/1，这个序列称为n级法雷数列，以Fn表示。如F5为：1/5, 1/4, 1/3, 2/5, 1/2, 3/5, 2/3, 3/4, 4/5。

那么我们怎么求得法雷数列呢？我们可以总结出法雷数列的性质如下：

给你一个数n（比如5），则法雷数列的首尾分别为0/1，1/1，则法雷数列中间的那个元素必定为(0+1)/(1+1)=1/2，介于0/1和1/2之间的法雷数列的中间一个为(0+1)/(1+2)=1/3，介于1/2和1/1中间的一个味(1+1)/(1+2)=2/3，由此我们发现规律，假设起始上下界为a/b，c/d，则中间的必为(a+c)/(c+d)，可以一直分下去，直到c+d>n退出，这就如同我们用二分法一样，因此求法雷数列的一种方法如下：

#include <stdio.h>
int n;
void pre(int a,int b,int c,int d){//二分法
	if(b+d>n) return ;
	pre(a,b,a+c,b+d);
	printf("%d/%d,",a+c,b+d);
	pre(a+c,b+d,c,d);
}
int main(){
	while(~scanf("%d",&n)){
		printf("0/1,");
		pre(0,1,1,1);
		printf("1/1\n");
	}
	return 0;
}

练习题：http://acm.nyist.net/JudgeOnline/problem.php?pid=334，程序如上。之前没有看过这种方法，写了一个很挫的方法，几年一下：

#include<string.h>
#include<stdlib.h>
struct node{
    int x,y;
}set[100100];
int check(int x, int y){
    int i,j,mid=x/2;
    if(y%x==0) return 0;
    
    for(i=2;i<=mid;i++){
        if(y%i==0&&x%i==0) return 0;
    }
    return 1;
}
int cmp(const void *a, const void *b){
    double x=(*(struct node*)a).x*1.0/(*(struct node*)a).y;
    double y=(*(struct node*)b).x*1.0/(*(struct node*)b).y; 
    return x>y?1:-1;
}
int main(){
    int i,j,n,k;
    while(scanf("%d",&n)!=EOF){
        for(i=1,k=1;i<=n;i++){
            for(j=1;j<i;j++){
                if(j==1||check(j,i)){
                    set[k].x=j;
                    set[k++].y=i;
                }
            }
        }
        set[0].x=0,set[0].y=1;
        set[k].x=1,set[k++].y=1;
        qsort(set,k,sizeof(set[0]),cmp);
        for(i=0;i<k-1;i++){
            printf("%d/%d,",set[i].x,set[i].y);
        }
        printf("%d/%d\n",set[i].x,set[i].y);
    }
    return 0;
}

 

2、欧拉函数。

欧拉函数：在数论，对正整数n，欧拉函数是少于或等于n的数中与n互质的数的数目。

计算方法：φ(x)=x(1-1/p1)(1-1/p2)(1-1/p3)(1-1/p4)…..(1-1/pn),其中p1, p2……pn为x的所有质因数，x是不为0的整数。φ(1)=1（唯一和1互质的数就是1本身）。 (注意：每种质因数只一个。比如12=2*2*3 ，那么φ（12）=12*（1-1/2）*(1-1/3)=4。

你可能回问，这和法雷数列有啥关系？可以这样解释：给出一个数n，如果要求n级法雷数列的元素个数，则根据法雷数列元素的性质（不可约分数），n级法雷数列的个数=2到n之间的自然数欧拉值的和,。

看一道题：http://acm.nyist.net/JudgeOnline/problem.php?pid=335

那么我们想到用累加的方法计算n级法雷数列的个数，但是直接累加的话，很不幸，超时了。。。如下：

#include <stdio.h>
/*
int eular(int n){//第一种计算欧拉函数的方法，不太懂
	int ret=1,i;
	for(i=2;i*i<=n;i++)
	if(n%i==0){
		n/=i,ret*=i-1;
		while(n%i==0)
			n/=i,ret*=i;
	}
	if(n>1)
		ret*=n-1;
	return ret;
}
*/
int eular(int n){//这种事按照公式来的计算欧拉函数的方法
    double sum=n;
    int i;
    for(i=2;i<=n;i++){
        if(n%i==0){
            sum*=(i-1)*1.0/i;
            while(n%i==0){    
                n=n/i;
            }
        }
    }
    return (int)sum;
}
int a[1500001]={0,0,1};

int main(){
	int n,i;
	for(i=2;i<1500001;i++)
		a[i]=a[i-1]+eular(i);
	while(~scanf("%d",&n)){
		printf("%d\n",a[n]+1);
	}
	return 0;
}

 既然超时了，我们想能不能改进一下算法呢？根据上面我们求欧拉函数值得方法，如果我们能实现直到1到N之间哪些数是素数（筛法可以快速求得），则这个素数的任意倍数的欧拉函数公式里面都有这个质因子，所以我们可以这样做：

#include<stdio.h>
#define N 1500001
//__int64 prime[N]={0},eular[N],falei[N]={0};
long long prime[N]={0},eular[N],falei[N]={0};
int main(){
    //__int64 n,i,j;
	long long n,i,j;
    prime[0]=prime[1]=1;
    for(i=2;i<N;i++){//筛法求得素数表 
        if(prime[i]==0)
            for(j=i*2;j<N;j+=i){
                prime[j]=1;
            }
    }
    /*φ(x)=x(1-1/p1)(1-1/p2)(1-1/p3)(1-1/p4)..(1-1/pn),先准备好 x*/ 
    for(i=1;i<N;i++){
        eular[i]=i;
    }
    for(i=1;i<N;i++){
		if(prime[i]==0)
			for(j=i;j<N;j+=i)
    /*φ(x)=x(1-1/p1)(1-1/p2)(1-1/p3)(1-1/p4)..(1-1/pn),每找到一个j的质因子，则乘到 φ(x)上面*/
				eular[j]=eular[j]*(i-1)/i;
    }
    for(i=1,falei[0]=1;i<N;i++){//falei[n]=eular[1] + eular[1]+...eular[n] +1
        falei[i]=falei[i-1]+eular[i];
    }
    while(scanf("%lld",&n)!=EOF){ 
        printf("%lld\n",falei[n]);
    }
    system("pause");
    return 0;
}