http://acm.hdu.edu.cn/showproblem.php?pid=1496
分析:a*x1^2+b*x2^2+1000000=-(c*x3^2+d*x4^2)+1000000(2*58*100*100)
#include<iostream>
#include<cstring>
#include<cstdio>
using namespace std;
const int NM=100;
int hash[2000005];
int main()
{
int i,j,ans,A,B,C,D;
int cc[NM+5];
for(i=1;i<=NM;i++)
cc[i]=i*i;
while(scanf("%d%d%d%d",&A,&B,&C,&D)!=EOF)
{
if((A>0&&B>0&&C>0&&D>0)||(A<0&&B<0&&C<0&&D<0))
{
printf("0\n");
continue;
}
memset(hash,0,sizeof(hash));
ans=0;
for(i=1;i<=NM;i++)
for(j=1;j<=NM;j++)
hash[A*cc[i]+B*cc[j]+1000000]++;
for(i=1;i<=NM;i++)
for(j=1;j<=NM;j++)
ans+=hash[-(C*cc[i]+D*cc[j])+1000000];
printf("%d\n",ans*16);
}
return 0;
}
外:冲突法减小内存消耗
#include<stdio.h>
#include<memory.h>
#define MAX 50021
int f[MAX],g[MAX];
int hash(int k)
{
int t=k%MAX;
if(t<0)
t+=MAX;
while(f[t]!=0&&g[t]!=k)
t=(t+1)%MAX;
return t;
}
int main()
{
int a,b,c,d,p,i,j,s,n,t[101];
for(i=1;i<=100;i++)
t[i]=i*i;
while(scanf("%d%d%d%d",&a,&b,&c,&d)>0)
{
if(a>0&&b>0&&c>0&&d>0||a<0&&b<0&&c<0&&d<0)
{
printf("0\n");
continue;
}
memset(f,0,sizeof(f));
n=0;
for(i=1;i<=100;i++)
for(j=1;j<=100;j++)
{
s=a*t[i]+b*t[j];
if(s==9)
int t1=0;
p=hash(s);
g[p]=s;
f[p]++;
}
for(i=1;i<=100;i++)
for(j=1;j<=100;j++)
{
s=-(c*t[i]+d*t[j]);
p=hash(s);
n+=f[p];
}
printf("%d\n",n*16);
}
}