lightOJ 1366 Pair of Touching Circles

枚举圆心的相对坐标，看他们间的距离是否为整数。若是整数则枚举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;
}