哈罗斯-法雷 数列

1.定义

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

如：         F3：0/1,1/3,1/2,2/3,1/1

          F4：0/1,1/4,1/3,1/2,2/3,3/4,1/1

F5：0/1,1/5, 1/4, 1/3, 2/5, 1/2, 3/5, 2/3, 3/4, 4/5,1/1

2.性质

1）Fn的项数Xn=1+[p=1~n]φ(p),[]:1~nφ(p)的和，即：F1有两项：0/1,1/1；F3=1+1+1+2

例如：poj 3090

for(i=1;i<N;i++)
		a[i]=i;
	for(i=2;i<N;i++)//N
	{
		if(a[i]==i)
		{
			for(j=i;j<N;j+=i)
				a[j]=(a[j]/i)*(i-1);
		}
	}
	b[1]=3;
	for(i=2;i<N;i++)
		b[i]=b[i-1]+2*a[i];

2）Fn中任意连续两项a/b<c/d，都有bc-ad=1

用这个性质可将实数化分数，详情请戳百度：http://baike.baidu.com/view/7084262.htm

例如：hdu 4180

求离所给分数最近的一个分数，用扩展欧几里得算法解方程：bc-ad=1

#include<iostream>
#include<cstdio>
using namespace std;
__int64 x,y,r;
void ex_Eulid(__int64 a,__int64 b)
{
    if(b==0)
    {x=1;y=0;r=a;}
    else
    {
        ex_Eulid(b,a%b);
        __int64 t=x;
        x=y;
        y=t-a/b*y;
    }
}
int main()
{
    __int64 T,a,b,n1,n2;
    scanf("%I64d",&T);
    while(T--)
    {
        scanf("%I64d/%I64d",&a,&b);
        ex_Eulid(a,b);
        if(r!=1) 
            printf("%I64d/%I64d\n",a/r,b/r);
        else if(x==1)
            printf("1/%I64d\n",b-1);
        else
        {
            n1=(x%b+b)%b;
            n2=(-x%b+b)%b;
            if(n1>n2)
                printf("%I64d/%I64d\n",(-y%a+a)%a,n1);
            else
                printf("%I64d/%I64d\n",(y%a+a)%a,n2);
        }
    }
    return 0;
}