hdu 3244 原题链接:点击打开链接
题意:将4*(1+2+..+n)作为直角三角形一条直角边,要求该三角形的最小斜边,及此时三角形的周长-3
本原勾股数组
a^2+b^2=c^2
a=s*t;
b=(s^2-t^2)/2;
c=(s^2+t^2)/2;
已知一条直角边,当另一条直角边最小时,斜边最小,
b=2*(1+n)*n=(s^2-t^2)/2;-->s^2=4*n*(1+n)+t^2
a^2=s^2*t^2=(4*n*(1+n)+t^2)*t^2;
当t=1时,a最小,
s=2*n+1;
c=2*n*n+2*n+1;
周长=4*n*n+6*n+2;
code:
#include<stdio.h>
int main()
{
__int64 n;
while(scanf("%I64d",&n)!=EOF)
{
printf("%I64d %I64d\n",4*n*n+6*n-1,2*n*n+2*n+1);
}
return 0;
}
poj 1305 原题链接:点击打开链接
题意:求在1~n中有多少个互质三元组(a,b,c),a^2+b^2=c^2, a,b,c两两互质,再求有多少数与这些互质三元组没有任何关系(不是他们任何数的倍数)
a=s*t,经证明知 s与t是奇数,且互质,暴力枚举s,t
注意数据类型
code:
#include<stdio.h>
#include<math.h>
#include<string.h>
int flag[1000005];
__int64 gcd(__int64 a,__int64 b)
{
return b==0?a:gcd(b,a%b);
}
int main()
{
__int64 n,ans1,ans2,a,b,c,s,t,i;
while(scanf("%I64d",&n)!=EOF)
{
ans1=0;ans2=0;
memset(flag,0,sizeof(flag));
for(t=1;t<=n;t+=2)
for(s=t+2;s<=n;s+=2)
{
if(gcd(s,t)!=1)
continue;
a=s*t;
b=(s*s-t*t)/2;
c=(s*s+t*t)/2;
if(c>n)
break;
ans1++;
for(i=1;i*c<=n;i++)
{flag[i*a]=1;flag[i*b]=1;flag[i*c]=1;}
}
for(i=1;i<=n;i++)
if(!flag[i])
ans2++;
printf("%I64d %I64d\n",ans1,ans2);
}
return 0;
}