UVA 10214(p339)----Trees in a Wood

最新推荐文章于 2024-11-27 23:00:53 发布

wang2147483647

最新推荐文章于 2024-11-27 23:00:53 发布

阅读量384

点赞数

CC 4.0 BY-SA版权

本文链接：https://blog.youkuaiyun.com/wang2147483647/article/details/50749353

算法竞赛入门经典（第二版）----数学概念与方法专栏收录该内容

39 篇文章

订阅专栏

本文介绍了一种基于Euler函数的算法实现方法，通过预处理得到最小质因数和欧拉函数值，并利用这些值解决特定数学问题。程序使用C++编写，包括了质因数分解、求最大公约数等实用函数。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

#include<iostream>
#include<cstdio>
using namespace std;
typedef unsigned long long ULL;
const int maxn=2010;
int phi[maxn+50];
int mindiv[maxn+50];
void prepare()
{
    for(int i=1; i<maxn; i++)
        mindiv[i]=i;
    for(int i=2; i*i<maxn; i++)
        if(mindiv[i]==i)
        {
            for(int j=i*i; j<maxn; j+=i)
                mindiv[j]=i;
        }
    phi[1]=1;
    for(int i=2; i<maxn; i++)
    {
        phi[i]=phi[i/mindiv[i]];
        if((i/mindiv[i])%mindiv[i]==0)
            phi[i]*=mindiv[i];
        else
            phi[i]*=mindiv[i]-1;
    }
}
int gcd(int a,int b)
{
    return b==0?a:gcd(b,a%b);
}
int main()
{
    int a,b;
    prepare();
    while(scanf("%d%d",&a,&b)==2&&(a+b))
    {
        ULL all=(ULL)4*a*b+2*(a+b);
        ULL ans=0;
        for(int i=1;i<=a;i++)
        {
            ans+=(b/i)*phi[i];
            int tmp=b%i;
            for(int j=1;j<=tmp;j++)
                if(gcd(i,j)==1) ans++;
        }
        ans*=4;ans+=4;
        printf("%.7lf\n",(double)ans/(double)all);
    }
    return 0;
}