UVa Problem 10233 Dermuba Triangle （Dermuba 三角）

// Dermuba Triangle （Dermuba 三角）
// PC/UVa IDs: 111207/10233, Popularity: C, Success rate: high Level: 2
// Verdict: Accepted
// Submission Date: 2011-11-01
// UVa Run Time: 0.016s
//
// 版权所有（C）2011，邱秋。metaphysis # yeah dot net
//
// [解题方法]
// 以编号为 0 的房子中心为原点计算其他房子中心的相对坐标，使用平面坐标系的距离公式计算其距离即可。

#include <iostream>
#include <iomanip>
#include <cmath>

using namespace std;

struct point
{
double x;
double y;
};

double calDistance(unsigned n, unsigned m)
{
point nn, mm;

// 计算编号为 n 的房子的横坐标和纵坐标，start 为 n 所在行的起始房子编号。
long long start = (long long)sqrt(n);
nn.x = (n - start * start - start) / 2.0;
nn.y = 0.0;
if (start)
{
nn.y = (start - 1) * (sqrt(3) / 2.0);
nn.y += ((start + n) & 1 ? (sqrt(3) / 3.0) : (sqrt(3) / 2.0));
}

// 计算编号为 m 的房子的横坐标和纵坐标。
start = (long long)sqrt(m);
mm.x = (m - start * start - start) / 2.0;
mm.y = 0.0;
if (start)
{
mm.y = (start - 1) * (sqrt(3) / 2.0);
mm.y += ((start + m) & 1 ? (sqrt(3) / 3.0) : (sqrt(3) / 2.0));
}

// 计算距离。
return sqrt(pow(nn.x - mm.x, 2) + pow(nn.y - mm.y, 2));
}

int main(int ac, char *av[])
{
unsigned n, m;

cout.precision(3);
cout.setf(ios::fixed | ios::showpoint);

while (cin >> n >> m)
cout << calDistance(n, m) << endl;

return 0;
}