<算法导论>第八章4 水壶(Kettle)


8-4  水壶

   假设给定了n个红色的水壶和n个蓝色的水壶,它们的形状和尺寸都不相同。所有红色水壶中所盛水的量都不一样,蓝色水壶也是一样。此外,对于每个红色的水壶,都有一个对应的蓝色水壶,两者所盛的水量是一样的。反之亦然。

  你的任务是将所盛水量一样的红色水壶和蓝色水壶找出来。为了达到这一目的,可以执行如下操作:挑选出一对水壶,其中一个是红色的,另一个是蓝色的:将红色水壶中倒满水;再将水倒入到蓝色的水壶中。通过这个操作,可以判断出来这两只水壶的容量哪一个大,或者是一样大。假设这样的比较需要一个时间单位。你的目标是找出一个算法,它通过执行最少次数的比较,来确定分组和配对问题。记住不能直接比较两个红色的或两个蓝色的水壶。

  a)给出一个确定型的算法,它利用Θ(n^2)次比较来完成水壶的配对。

  b)证明:对于一个可以解决本问题的算法,必须执行的比较次数的下界为Ω(nlgn)。

  c)给出一个随机化的算法,其期望的比较次数为O(nlgn),并证明这个界是正确的。对你的算法来说,最坏情况下的比较次数是什么?

a)先用朴素的方法,拿出第1个红色水壶,然后从n个蓝色水壶中去配对,需要Θ(n)次比较。重复上述操作,知道所有的红色水壶都完成了配对。总的时间复杂度为Θ(n^2)

b)采用决策树的方法进行分析,总共可能的配对有n!种,每次比较的产生的分支有3个(大于、相等或小于),若决策树的高度为h,则Ω(nlgn)就是比较次数的下界。

c) 由于每个红色的水壶和蓝色的水壶都有一个配对,因此理论上我们只需要分别将其进行排序即可。但是实际问题中,是不容许两个各红色或蓝色的水壶进行比较的,必须颜色不一样的水壶才能进行比较。因此需要一定的技巧。

     整个排序的过程中采用快速排序的思想,现对最关键也是有所不同的PATITION部分,这个过程有点类似交叉验证,具体描述如下:

1. 从红色水壶序列中随机选择一个作为枢轴元素

2. 在蓝色水壶序列中查找和该红色水壶容量一样的那个作为枢轴元素

3. 利用红色水壶枢轴元素对蓝色水壶序列进行PATITION。

4. 利用蓝色水壶枢轴元素对红色水壶序列进行PATITION。

如此PATIION部分就完成了,期望的时间复杂度为Θ(n)。整个过程类似快排迭代进行即可完成排序。

整个时间的复杂度为:时间复杂度为O(nlgn)。最坏的情况下时间复杂度为O(n^2)。

下面是一个简单的代码实现:

struct Kettle
{
    int I,V;//水壶的初始下标,容量
    Kettle():I(0),V(0){}
    Kettle(int _i,int _v):I(_i),V(_v){}
};

/*随机产生[L,R]上的随机数*/
int Rand(int L,int H)
{
    return rand()%(H-L+1) + L;
}
istream& operator>>(istream& in,Kettle& k)
{
    in>>k.I>>k.V;
    return in;
}
void swap(Kettle& lh,Kettle& rh)//重载swap函数
{
    swap(lh.I,rh.I);
    swap(lh.V,rh.V);
}

/*
A代表红色水壶
B代表蓝色水壶
L,R表示当前处理的数据在数组中的上下界
step 1:从A中随机选择一个壶r
step 2:从蓝色水壶中找到和r容量一样的蓝色水壶b
step 3:将r作为蓝色水壶Partition的分隔元素,b作为红色水壶Partition的分隔元素
*/
int Partition(Kettle* A,Kettle *B,int L,int R)
{
    int I = Rand(L,R),J;
    int bKey = A[I].V;//从红色中随机选取的分隔元素,做蓝色的分隔元素

    for(J=L; J<=R; ++J)//找B中和rKey一样的元素
        if(B[J].V==bKey)
        {
            cout<<"A["<<A[I].I<<"] = B["<<B[J].I<<"] = "<<bKey<<endl;
            break;
        }
    int rKey = B[J].V; //做红色的分隔元素

    int j,ir=L-1;
    swap(A[I],A[R]);
    for(j=L; j<R; ++j)//对红色水壶进行划分
    {
        if(A[j].V<=rKey)
        {
            ++ir;
            swap(A[ir],A[j]);
        }
    }
    swap(A[ir+1],A[R]);

    int ib=L-1;
    swap(B[J],B[R]);
    for(j=L; j<R; ++j)//对蓝色水壶进行划分
    {
        if(B[j].V<=bKey)
        {
            ++ib;
            swap(B[ib],B[j]);
        }
    }
    swap(B[ib+1],B[R]);
    return ir+1;//这里肯定 ib==ir
}

void Match(Kettle *A,Kettle *B,int L,int R)//整个过程和快排基本一致
{
    if(L < R)
    {
        int M = Partition(A,B,L,R);
        Match(A,B,L,M-1);
        Match(A,B,M+1,R);
    }
    else if(L==R) //单个的元素直接配对
        cout<<"A["<<A[L].I<<"] = B["<<B[L].I<<"] = "<<A[L].V<<endl;
}

void TestPartition()
{
    srand(11);
    int a[] = {2, 3, 7, 5, 1, 9, 6, 8, 4};
    int b[] = {4, 7, 9, 5, 1, 3, 8, 2, 6};
    Kettle A[9],B[9];
    for(int i=0;i<9;++i)
    {
        A[i].I = i; A[i].V = a[i];
        B[i].I = i; B[i].V = b[i];
    }
    Match(A,B,0,8);
}

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值