#include<stdio.h>
#include<stdlib.h>
#include<math.h>
 
typedef struct tag_point_t
{
    double x, y;
}point_t;
 
#define POINT_NUM 30
 
double calc_distance(const point_t *p, const point_t *q)
{
	double dx, dy;
 
	dx=p->x-q->x;
	dy=p->y-q->y;
	return sqrt(dx*dx+dy*dy);
}
 
int main(void)
{
	point_t point[POINT_NUM];
	unsigned long connection[POINT_NUM], A_connection;
	int i, j, is_changed, A, B;
 
	for(i=0;i<POINT_NUM;i++)
	{
		point[i].x=rand()/(RAND_MAX+1.0)*10.0;
		point[i].y=rand()/(RAND_MAX+1.0)*10.0;
		connection[i]=0;
	}
 
	A=rand()%POINT_NUM;
	do
	{
		B=rand()%POINT_NUM;
	}while(A==B);
 
	for(i=0;i<POINT_NUM;i++)
	{
		for(j=i;j<POINT_NUM;j++)
		{
			if(calc_distance(&point[i], &point[j])<=10.0)
			{
				connection[i]|=(1<<j);
				connection[j]|=(1<<i);
			}
		}
	}
 
	A_connection=connection[A];
	do
	{
		is_changed=0;
		for(i=0;i<POINT_NUM;i++)
		{
			if(((A_connection>>i)&1) && ((A_connection&connection[i])!=connection[i]))
			{
				A_connection|=connection[i];
				is_changed=1;
			}
		}
	}while(is_changed);
 
	for(i=0;i<POINT_NUM;i++)
	{
		printf("%2d : (%f,%f)\n", i, point[i].x, point[i].y);
	}
	printf("A=%d  B=%d\n", A, B);
 
	printf("移動%s\n", ((A_connection>>B)&1)?"可能":"不可能");
 
	return 0;
}
 

I2luY2x1ZGU8c3RkaW8uaD4KI2luY2x1ZGU8c3RkbGliLmg+CiNpbmNsdWRlPG1hdGguaD4KCnR5cGVkZWYgc3RydWN0IHRhZ19wb2ludF90CnsKICAgIGRvdWJsZSB4LCB5Owp9cG9pbnRfdDsKCiNkZWZpbmUgUE9JTlRfTlVNIDMwCgpkb3VibGUgY2FsY19kaXN0YW5jZShjb25zdCBwb2ludF90ICpwLCBjb25zdCBwb2ludF90ICpxKQp7Cglkb3VibGUgZHgsIGR5OwoKCWR4PXAtPngtcS0+eDsKCWR5PXAtPnktcS0+eTsKCXJldHVybiBzcXJ0KGR4KmR4K2R5KmR5KTsKfQoKaW50IG1haW4odm9pZCkKewoJcG9pbnRfdCBwb2ludFtQT0lOVF9OVU1dOwoJdW5zaWduZWQgbG9uZyBjb25uZWN0aW9uW1BPSU5UX05VTV0sIEFfY29ubmVjdGlvbjsKCWludCBpLCBqLCBpc19jaGFuZ2VkLCBBLCBCOwoKCWZvcihpPTA7aTxQT0lOVF9OVU07aSsrKQoJewoJCXBvaW50W2ldLng9cmFuZCgpLyhSQU5EX01BWCsxLjApKjEwLjA7CgkJcG9pbnRbaV0ueT1yYW5kKCkvKFJBTkRfTUFYKzEuMCkqMTAuMDsKCQljb25uZWN0aW9uW2ldPTA7Cgl9CgoJQT1yYW5kKCklUE9JTlRfTlVNOwoJZG8KCXsKCQlCPXJhbmQoKSVQT0lOVF9OVU07Cgl9d2hpbGUoQT09Qik7CgoJZm9yKGk9MDtpPFBPSU5UX05VTTtpKyspCgl7CgkJZm9yKGo9aTtqPFBPSU5UX05VTTtqKyspCgkJewoJCQlpZihjYWxjX2Rpc3RhbmNlKCZwb2ludFtpXSwgJnBvaW50W2pdKTw9MTAuMCkKCQkJewoJCQkJY29ubmVjdGlvbltpXXw9KDE8PGopOwoJCQkJY29ubmVjdGlvbltqXXw9KDE8PGkpOwoJCQl9CgkJfQoJfQoKCUFfY29ubmVjdGlvbj1jb25uZWN0aW9uW0FdOwoJZG8KCXsKCQlpc19jaGFuZ2VkPTA7CgkJZm9yKGk9MDtpPFBPSU5UX05VTTtpKyspCgkJewoJCQlpZigoKEFfY29ubmVjdGlvbj4+aSkmMSkgJiYgKChBX2Nvbm5lY3Rpb24mY29ubmVjdGlvbltpXSkhPWNvbm5lY3Rpb25baV0pKQoJCQl7CgkJCQlBX2Nvbm5lY3Rpb258PWNvbm5lY3Rpb25baV07CgkJCQlpc19jaGFuZ2VkPTE7CgkJCX0KCQl9Cgl9d2hpbGUoaXNfY2hhbmdlZCk7CgoJZm9yKGk9MDtpPFBPSU5UX05VTTtpKyspCgl7CgkJcHJpbnRmKCIlMmQgOiAoJWYsJWYpXG4iLCBpLCBwb2ludFtpXS54LCBwb2ludFtpXS55KTsKCX0KCXByaW50ZigiQT0lZCAgQj0lZFxuIiwgQSwgQik7CgoJcHJpbnRmKCLnp7vli5Ulc1xuIiwgKChBX2Nvbm5lY3Rpb24+PkIpJjEpPyLlj6/og70iOiLkuI3lj6/og70iKTsKCglyZXR1cm4gMDsKfQo=