HDU 3939 毕达哥拉斯三元组

题意:求形成直角三角形切个边长小于L的个数,并且三边两两互素。

也就是求前L本原毕达哥拉斯三元组解的个数。令,其中m>n,。

由上式得z为直角三角形斜边且z<=L

#include <iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
using namespace std;
#define maxn 1000008
bool isprime[maxn];
int prime[maxn],nprime,phi[maxn];
void getprime()
{
    long long i,j;
    memset(isprime,1,sizeof(isprime));
    nprime=0;
    for(i=2; i<maxn; i++)
        if(isprime[i])
        {
            prime[nprime++]=i;
            for(j=i*i; j<maxn; j+=i)
                isprime[j]=0;
        }
}
void getphi()
{
    for(int i=1; i<maxn; i++) phi[i]=i;
    for(int i=2; i<maxn; i+=2) phi[i]>>=1;
    for(int i=3; i<maxn; i+=2)
        if(phi[i]==i)
            for(int j=i; j<maxn; j+=i)
                phi[j]=phi[j]/i*(i-1);
}
int check[35],num;
long long ans;
void getcheck(int x)
{
    num=0;
    int d=x;
    if(isprime[x])
    {
        check[num++]=x;
        return;
    }
    for(int i=0; prime[i]*prime[i]<=x; i++)
        if(d%prime[i]==0)
        {
            check[num++]=prime[i];
            while(d%prime[i]==0)d/=prime[i];
        }
    if(d>1&&isprime[d]) check[num++]=d;
}
void dfs(int k,int r,int s,int n)
{
    if(k==num)
    {
        if(r&1) ans-=n/s;
        else ans+=n/s;
        return;
    }
    dfs(k+1,r,s,n);
    dfs(k+1,r+1,s*check[k],n);
}
int main()
{
    int t;
    long long l;
    getprime();
    getphi();
    scanf("%d",&t);
    while(t--)
    {
        scanf("%I64d",&l);
        ans=0;
        int m=sqrt(double(l-1));
        for(int i=m; i>0; i--)
        {
            int p=sqrt(double(l-(long long)i*i));
            if(i&1)
            {
                getcheck(i);
                if(p<i) dfs(0,0,1,p>>1);
                else dfs(0,0,1,i>>1);
            }
            else
            {
                if(i<p) ans+=phi[i];
                else getcheck(i),dfs(0,0,1,p);
            }
        }
        printf("%I64d\n",ans);
    }
    return 0;
}
时间: 2024-11-05 14:21:41

HDU 3939 毕达哥拉斯三元组的相关文章

JScript 函数

js|jscript|函数 Microsoft Jscript 函数执行操作,也可以返回值.某些时候是计算或比较的结果.函数又被称为"全局方法". 一个函数中包含有几个操作.这样可使得代码更合理化.可以写一组语句并给其命名,然后通过调用它并传递其需要的信息来运行整组语句. 给函数传递信息可以把信息放在函数名称后面的圆括号中.传递给函数的信息称作参数.某些函数根本不带任何参数,而其他函数带一个或者多个参数.在某些函数中,参数的个数取决于如何使用该函数. Jscript 支持两种函数:一类

《数论概论(原书第4版)》一第2章 勾 股 数 组

第2章 勾 股 数 组 毕达哥拉斯定理(即勾股定理)是中学生"喜爱"的公式,它表明任一个直角三角形(如图21所示)的两条直角边长的平方和等于斜边长的平方.用公式表示就是 因为对数论(即自然数理论)感兴趣,所以我们会问是否存在毕达哥拉斯三角形,它的所有边长都是自然数.有许多这样的三角形,最著名的例子是边长为3,4,5的三角形(即中国古代数学家发现的"勾广三,股修四.径隅五").下面是前几个例子: 对这些毕达哥拉斯三元组(下称勾股数组)的研究在毕达哥拉斯时代以前很久就

hdu 1527

题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=1527 hint:威佐夫博弈 基本类似于模板 #include <iostream> #include <cmath> #include <cstdio> using namespace std; const double q = (1 + sqrt(5.0)) / 2.0; // 黄金分割数 int Wythoff(int a, int b) { if (a > b)

hdu 2551 竹青遍野

题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=2551 hint:就是读懂题就行了 #include <iostream> #include <cstdio> using namespace std; typedef long long LL; LL data[1005]; int main() { data[0]=0; for(int i=1; i<1005; i++) data[i]+=data[i-1]+i*i*i; LL

hdu 2054 A == B?

http://acm.hdu.edu.cn/showproblem.php?pid=2054 此题巨坑,刚开始我以为是简单的水题,就用strcmp过, but错了,后来经过我苦思冥想,结果还有几组数据 0.0 和 0,1.000和1.0 , 但是我不太确定前面的0是不是有作用我还是写了,但是有人过的时候,前面的0没考虑比如: 002和2可能是相等的,也可能是不想等的所以不用判断,只能说明hdu数据不是很强啊,嘿嘿 代码如下: #include <iostream> #include <c

hdu 4430 Yukari&#039;s Birthday

点击打开链接hdu 4430 思路:枚举r+二分k 分析: 1 题目要求的是找到一组最小的r*k,如果r*k相同那么就找r最小的. 2 很明显k>=2,根据n <= 10^12,那么可以知道r的最大值r<50,所以只要枚举枚举r的值,然后二分k的大小找到所有的解,存入一个结构体里面,然后在对结构体排序,那么这样就可以得到最后的ans 3 注意题目说了中心点最多一个蜡烛,所以写二分的时候应该注意判断的条件: 4 还有可能计算得到结果超了long long直接变成负数所以应该对或则个进行判断

hdu 1238 Substrings

点击打开链接hdu 1238 思路:kmp+暴力枚举子串 分析: 1 题目要求找到一个子串x,满足x或x的逆串是输入的n个字符串的子串,求最大的x,输出x的长度 2 题目的n最大100,每一个字符串的最大长度为100,那么暴力枚举子串就是o(n^2)才10000肯定是不会超时的,但是由于这里涉及到了逆串的问题,所以我们应该还要求出n个子串的逆串,然后在求最大的x. 代码: #include<iostream> #include<algorithm> #include<cstd

hdu 1857 Word Puzzle

点击打开链接hdu 1857 思路:字典树 分析: 1 题目要求的是给定的单词第一个字母在这个矩形里面的最小的坐标 2 矩形的最大500*500,单词的来源有三个方向,并且单词的起点和终点在矩形之内都是可能的.所以的如果利用枚举矩形之内的单词,那么肯定是超内存的 3 所以我们必须考虑另一种的方法就是对单词进行建字典树,那么我们只要去枚举单词的可能的起点,然后进行查找相应的单词是不是在树上,如果是的话就标记一下当前的坐标. 4 注意由于单词的来源有三个方向,但是因为要求的如果下相同的情况下要求坐标

hdu 1595 find the longest of the shortest

点击打开链接hdu 1595 思路:最短路+优先队列+Dijstra+枚举边 分析: 1 题目要求的是删掉一条边之和求出的最短路中的最大值. 2 很明显,肯定是要先求出原图的最短路并且记录父亲节点.现在我们可以想,如果要枚举所有的边,显然这个是不可能的实现的.所以我们仔细分析可以知道其实能够对最短路产生影响的就是原图最短路上的边,所以我们只需要去枚举删除最短路径上面边然后求最短路即可,最后得到ans 3 这一题的n <= 1000 , m<=n*(n-1)/2 , 刚开始我用的SPFA,然后就