枚举圆心的相对坐标,看他们间的距离是否为整数。若是整数则枚举1~(length-1) 的半径。看要多大的矩形能容纳这两个圆。然后在大矩形中找有几个这样的矩形。
貌似两个圆形不在水平线或者竖直线上的情况很少,列出来存在数组中,减少两次循环花的时间会更少。
#include <iostream>
#include <cstdio>
#include <cstring>
#include <queue>
#include <stack>
#include <algorithm>
#include <cmath>
#include <map>
#define LL long long
#define DB double
using namespace std;
const int N = 1009;
int dis(int x,int y)
{
int tmp = (int)sqrt(x*x+y*y+0.1);
if(tmp*tmp==x*x+y*y) return tmp;
return 0;
}
int v[N][N];
void init()
{
for(int i=0;i<=500;i++)
for(int j=0;j<=500;j++)
v[i][j] = dis(i,j);
}
int h,w;
int in(int x,int y)
{
return (h-x+1)*(w-y+1);
}
void solve(int T)
{
LL ans = 0;
for(int i=0;i<=h/2;i++)
for(int j=0;j<=w/2;j++)
if(v[i][j])
{
for(int k=1;k<v[i][j];k++)
{
int y1 = min(-k,j-(v[i][j]-k)),y2 = max(k,j+(v[i][j]-k));
int x1 = min(-k,i-(v[i][j]-k)),x2 = max(k,i+(v[i][j]-k));
int x = x2-x1,y=y2-y1;
if(x>h||y>w) continue;
int tmp = in(x,y);
if(i*j)
tmp<<=1;
ans += tmp;
}
}
printf("Case %d: %lld\n",T,ans);
}
int main()
{
#ifndef ONLINE_JUDGE
freopen("in.txt","r",stdin);
#endif
int cas,T=1;scanf("%d",&cas);
init();
while(cas--)
{
scanf("%d%d",&h,&w);
solve(T);T++;
}
return 0;
}