

Posted by EHXM. Posted in " 경험/알고리즘 "2009/10/26 01:56
With a growing desire for modernization in our increasingly larger cities comes a need for new street designs. Chris is one of the unfortunate city planners responsible for these designs. Each year the demands keep increasing, and this year he has even been asked to design a completely new city.
More work is not something Chris needs right now, since like any good bureaucrat, he is extremely lazy. Given that this is a character trait he has in common with most computer scientists it should come as no surprise that one of his closest friends, Paul, is in fact a computer scientist. And it was Paul who suggested the brilliant idea that has made Chris a hero among his peers: Fractal Streets! By using a Hilbert curve, he can easily fill in rectangular plots of arbitrary size with very little work.
A Hilbert curve of order 1 consists of one “cup”. In a Hilbert curve of order 2 that cup is replaced by four smaller but identical cups and three connecting roads. In a Hilbert curve of order 3 those four cups are in turn replaced by four identical but still smaller cups and three connecting roads, etc. At each corner of a cup a driveway (with mailbox) is placed for a house, with a simple successive numbering. The house in the top left corner has number 1, and the distance between two adjacent houses is 10 m.
The situation is shown graphically in figure 2. As you can see the Fractal Streets concept successfully eliminates the need for boring street grids, while still requiring very little effort from our bureaucrats.
(a) Order 1 (b) Order 2 (c) Order 3
Figure 2: Hilbert curves of order 1, 2 and 3, with the location of the houses indicated.
As a token of their gratitude, several mayors have offered Chris a house in one of the many new neighborhoods built with his own new scheme. Chris would now like to know which of these offerings will get him a house closest to the local city planning office (of course each of these new neighborhoods has one). Luckily he will not have to actually drive along the street, because his new company “car” is one of those new flying cars. This high-tech vehicle allows him to travel in a straight line from his driveway to the driveway of his new office. Can you write a program to determine the distance he will have to fly for each offer (excluding the vertical distance at takeoff and landing)?
On the first line of the input is a positive integer, the number of test cases. Then for each test case:
• A line containing a three positive integers, n< 16 and h,o < 231, specifying the order of the Hilbert curve, and the house numbers of the offered house and the local city planning office.
For each test case:
• One line containing the distance Chris will have to fly to his work in meters, rounded to the nearest integer.
Exampl
etypedef struct{
int x;
int y;
}point;
point solve(long order,long num){
if(order == 1){
point a;
if(num == 1){
a.x = 1;
a.y = 1;
}
else if(num == 2){
a.x = 2;
a.y = 1;
}
else if(num == 3){
a.x = 2;
a.y = 2;
}
else if(num == 4){
a.x = 1;
a.y = 2;
}
return a;
}
long div = pow(2.0,2*(order-1));
FOR(i, 4){
if(num <= (i+1)*div){
if(i==0){
point tmp = solve(order-1, div - (num-i*div-1));
long double dx = (pow(2.0,order-1)+1)/2.0;
long double dy = dx;
long double cx = tmp.x - dx;
long double cy = dy - tmp.y;
long double fx = cy;
long double fy = -cx;
tmp.x = dx + fx;
tmp.y = dy - fy;
return tmp;
}
else if(i==1){
point tmp = solve(order-1, num-i*div);
tmp.x = tmp.x + pow(2.0,order-1);
return tmp;
}
else if(i==2){
point tmp = solve(order-1, num-i*div);
tmp.x = tmp.x + pow(2.0,order-1);
tmp.y = tmp.y + pow(2.0,order-1);
return tmp;
}
if(i==3){
point tmp = solve(order-1, div - (num-i*div-1));
long double dx = (pow(2.0,order-1)+1)/2.0;
long double dy = dx;
long double cx = tmp.x - dx;
long double cy = dy - tmp.y;
long double fx = -cy;
long double fy = cx;
tmp.x = dx + fx;
tmp.y = dy - fy;
tmp.y = tmp.y + pow(2.0,order-1);
return tmp;
}
}
}
}
int main(){
int testCase;
cin >> testCase;
FOR(testi, testCase){
long o,a,b;
cin >> o >> a >> b;
point pa = solve(o,a);
point pb = solve(o,b);
long double tmp = (pa.x-pb.x)*(pa.x-pb.x)+(pa.y-pb.y)*(pa.y-pb.y);
long result = (long)(sqrt(tmp)*10+0.5);
cout << result << endl;
}
return 0;
}
여러분의 커뮤니케이션을 기다리고 있습니다.

아이디어의 전쟁의 현장이었던 2010년 대한민국 매쉬업 경진대회에 다녀왔습니다. 이번 대회는 지난 2월 6일(토요일), 삼성동 코엑스 컨퍼런스룸 401에서 열렸습니다. 이번.....
2010년에 100가지가 넘는 안드로이드 폰 출시가 될 예정입니다. Mobile World Congress keynote에서 Google CEO Eric Schmidt의 연설.....
서울, 안양, 부산, 대구, 광주, 대전에 안드로이드 폰 체험 할 수 있는 곳이 있네요. 저는 코엑스 메가박스 입구에 있는 모토로라 체험 부스에서 우연히 모토로이를 만져보게 되.....
위 3D 갤러리는 http://www.fotoviewr.com/ 사이트의 Fotoviewr 입니다. Flex와 Papervision3D를 이용하여 위와같은 3D 갤러리를 구현해.....
Total : 115,868 Today : 224 Yesterday : 155