POJ 1305&&FZU 1669

两道求毕达哥拉斯三元组的个数的题目，POJ那道是求N的范围内，满足毕达哥拉斯方程的数的组数和不满足的组数，FZU那道是求x+y+z<=L的情况下本原毕达哥斯拉三元组的组数。思路大同小异，首先假设x,y,z是一组三元组，把它们写成形如x=m*m-n*n y=2*m*n   z=m*m+n*n的形式 然后枚举m n暴力求解 三元组，POJ那道不是求本原三元组所以开个VIS标记本原三元组的K倍也是符合条件的三元组，然后未标记的就是范围内不满足条件的三元组的组数，而FZU那道直接枚举就好了，注意判断是否超过界限了就好，至于什么是毕达哥拉斯方程组，什么是本原三元组，可以看数论导引这本书里面的讲解，虽然各种不明觉厉，但是参照了一个学长的读书笔记还是看懂了一点点。这两道题都很水，没办法蒟蒻只会切水题，不像别人啥都会，好忧桑啊，赶脚马上就要退役了的即视感TVT，在有限的时间内再学点知识好了，提高智商，发散思维也是阔以的嘛TVT

POJ 1305

#include<cstdio>
#include<cstring>
#include<cmath>
#include<iostream>
#include<algorithm>
using namespace std;
const int maxn=1000000;
bool vis[maxn];
int n;
int gcd(int a,int b)
{
    if(b==0)
        return a;
    else return gcd(b,a%b);
}
void slove(int t)
{
    memset(vis,0,sizeof(vis));
    int ans1=0,ans2=0;
    int x,y,z;
    for(int n=1;n*n<=t;n++)
    {
        for(int m=n+1;m*m<=t;m++)
        {
            if(n*n+m*m>t)
            {
                break;
            }
            if(n%2!=m%2)
            {
                if(gcd(n,m)==1)
                {
                    x=m*m-n*n;
                    y=2*m*n;
                    z=m*m+n*n;
                    ans1++;
                    for(int i=1;;i++)
                    {
                        if(i*z>t)
                        {
                            break;
                        }
                        vis[i*x]=vis[i*y]=vis[i*z]=1;
                    }
                }
            }
        }
    }
    for(int i=1;i<=t;i++)
    {
        if(vis[i]==0)
        {
            ans2++;
        }
    }
    cout<<ans1<<" "<<ans2<<endl;
}

int main()
{
    while(cin>>n)
    {
        slove(n);
    }
    return 0;
}

FZU 1669

#include<cstdio>
#include<cstring>
#include<cmath>
#include<iostream>
#include<algorithm>
using namespace std;

const int maxn=2000000;
bool vis[maxn];
int n;
int gcd(int a,int b)
{
    if(b==0)
        return a;
    else
        return gcd(b,a%b);
}
void slove(int t)
{
    int ans=0;
    memset(vis,0,sizeof(vis));
    int x,y,z;

    for(int n=1;n*n<=t;n++)
    {
        for(int m=n+1;m*m<=t;m++)
        {
            if(m*m+n*n>t)
                break;
            if(n%2!=m%2)
            {
                if(gcd(n,m)==1)
                {
                    x=m*m-n*n;
                    y=2*m*n;
                    z=m*m+n*n;
                    for(int i=1; ;i++)
                    {
                        if(i*(x+y+z)>t)
                        {
                            break;
                        }
                        ans++;

                    }
                }
            }
        }
    }
    cout<<ans<<endl;
}
int main()
{
    while(cin>>n)
    {
        slove(n);
    }
    return 0;
}