可以说这道题目,涉及了很多知识点,是一道很好的题目。
首先这道题涉及了计算几何的知识,求圆心坐标的横距离进而求出圆心坐标;其次,这道题目数据不大,很容易想到dfs;最后这道题涉及到了全排列!理论上这道题还可以剪枝,步骤也不是很麻烦,但是不剪枝也是可以过的!
下面说一下我最开始的想法,这个想法很麻烦,最后wa了,简单说一下,就是从右向左判断,一次和圆相切,算出坐标,然后假定这个坐标,然后继续向左判断看看有没有和它相交的圆。话说这个想法很麻烦,然后自后求长度的时候,也求错了,其实找到最左边的坐标和最右边的坐标,相减就是长度,而我是把最右边的坐标作为长度,坐标超界的部分虽然做了处理,但是处理得也是有问题,我感觉!
说一下后来AC的的做法,这个做法是老大给我讲的,而且操作简单,代码写出来也很精简。首先,在某个位置,枚举圆,我们是一定要把这个圆放在这的,只是要确定圆心坐标,保证这个圆至少和一个圆相切,并且和其他的圆都不会相交。那么,我么已经放好的圆的左边的圆开始,让当前圆一次与这些圆相切,并求出和每个圆相切的坐标,然后找到最右边的x的值,就是当前可以放的位置,当然y值就是当前圆的半径。其次,我们要论证一下,为什么这么做是对的,首先我们要保证以当前这样的排列,使得相同排列的情况长度最小,那么就要都与左边的圆相切,这个没什么说的;那么为什么要选最大的x坐标,试想,求出的所有x坐标,都是与某个圆相切的出来的,那么比如一共有8个圆,放到第3个为位置的时候,求出与第一个圆相切的横坐标是1,于第二个圆相切的横坐标是2,那么相切代表它刚刚好不相交,如果我们选择了1的这个横坐标,把圆放在这个位置,那么这个圆必然会和第二个圆相交,因此,一定要选择求出的所有的坐标的最右边的。还有一个常犯的错误就是只求和前一个圆相切的坐标,这样的话,如果,倒数第二个圆十分大的话,不能保证当前圆和和倒数第二个不相交,一次类推,倒数第三个圆也有可能很大很大.....
下面列出一些圆的位置的特殊情况:
第一附图是说明可能会于前面的任何圆相交,第二幅图说明只判断和最后一个圆相切是不对的。
还有要注意的就是求长度
没有剪枝的代码:
#include <cstdio>
#include <cstring>
#include <cmath>
#include <algorithm>
#include <iostream>
using namespace std;
#define eps 1e-9
const int N = 10;
const double INF = 100000000.0;
int T, n;
int id[N];
bool vis[N];
double r[N], x[N], y[N], ans;
void dfs( int now )
{
if (now == n) {
double lx = x[0], rx = x[0], tmp;
for ( int i = 0; i < n; ++i ) {
//cout << x[i] << " ";
lx = min( x[i]-r[id[i]], lx );
rx = max( x[i]+r[id[i]], rx );
}
//cout << endl << lx << ' ' << rx << endl;
tmp = rx - lx;
if ( tmp + eps < ans ) ans = tmp;
}
for ( int i = 0; i < n; ++i )
if ( !vis[i] ) {
double mx = x[0], my = r[i];
for ( int j = 0; j < now; ++j ) {
double rj = r[id[j]], xx;
xx = sqrt((rj+r[i])*(rj+r[i])-(rj-r[i])*(rj-r[i])) + x[j];
mx = max ( xx, mx );
}
x[now] = mx, y[now] = r[i], id[now] = i;
vis[i] = true;
dfs( now+1 );
vis[i] = false;
}
}
int main()
{
scanf("%d", &T);
while ( T-- ) {
scanf("%d", &n);
for ( int i = 0; i < n; ++i ) scanf("%lf", &r[i]);
ans = INF;
memset( vis, 0, sizeof(vis) );
for ( int i = 0; i < n; ++i ) {
int cnt = 0;
vis[i] = true;
id[cnt] = i;
x[cnt] = y[cnt] = r[i];
dfs( cnt+1 );
//cout << "...................."<< endl;
vis[i] = false;
}
printf("%.3lf\n", ans);
}
}
剪枝的代码(一旦发现当前的长度大于以确定的长度,就返回):
#include <cstdio>
#include <cstring>
#include <cmath>
#include <algorithm>
#include <iostream>
using namespace std;
#define eps 1e-9
const int N = 10;
const double INF = 100000000.0;
int T, n;
int id[N];
bool vis[N];
double r[N], x[N], y[N], ans;
void dfs( int now, double tmp )
{
if (now == n) {
//double lx = x[0], rx = x[0], tmp;
/*for ( int i = 0; i < n; ++i ) {
//cout << x[i] << " ";
lx = min( x[i]-r[id[i]], lx );
rx = max( x[i]+r[id[i]], rx );
}
//cout << endl << lx << ' ' << rx << endl;
tmp = rx - lx;*/
if ( tmp + eps < ans ) ans = tmp;
return;
}
if ( tmp > ans + eps ) return;
for ( int i = 0; i < n; ++i )
if ( !vis[i] ) {
double mx = x[0], my = r[i], lx = x[0], rx = x[0];
for ( int j = 0; j < now; ++j ) {
double rj = r[id[j]], xx;
xx = sqrt((rj+r[i])*(rj+r[i])-(rj-r[i])*(rj-r[i])) + x[j];
mx = max ( xx, mx );
lx = min( x[j]-r[id[j]], lx );
rx = max( x[j]+r[id[j]], rx );
}
lx = min( mx-r[i], lx );
rx = max( mx+r[i], rx );
x[now] = mx, y[now] = r[i], id[now] = i;
vis[i] = true;
dfs( now+1, rx-lx );
vis[i] = false;
}
}
int main()
{
scanf("%d", &T);
while ( T-- ) {
scanf("%d", &n);
for ( int i = 0; i < n; ++i ) scanf("%lf", &r[i]);
ans = INF;
memset( vis, 0, sizeof(vis) );
for ( int i = 0; i < n; ++i ) {
int cnt = 0;
vis[i] = true;
id[cnt] = i;
x[cnt] = y[cnt] = r[i];
dfs( cnt+1, r[i]*2 );
//cout << "...................."<< endl;
vis[i] = false;
}
printf("%.3lf\n", ans);
}
}