UVa 10012 How Big is It? (计算几何+DFS)

可以说这道题目，涉及了很多知识点，是一道很好的题目。

首先这道题涉及了计算几何的知识，求圆心坐标的横距离进而求出圆心坐标；其次，这道题目数据不大，很容易想到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);
    }
}