算法复习——分治法重点问题

一、众数问题

问题描述：
　　给定含有n 个元素的多重集合S，每个元素在S 中出现的次数称为该元素的重数。多重集S 中重数最大的元素称为众数。
　　例如，S={1，2，2，2，3，5}。多重集S 的众数是2，其重数为3。
编程任务：
　　对于给定的由n 个自然数组成的多重集S，编程计算S 的众数及其重数。
数据输入：
　　输入数据由文件名为input.txt 的文本文件提供。文件的第1 行多重集S 中元素个数n；接下来的n 行中，每行有一个自然数。
结果输出:
　　程序运行结束时，将计算结果输出到文件output.txt 中。输出文件有2 行，第1 行给出众数，第2 行是重数。
输入文件示例          输出文件示例
input.txt             output.txt
6                      2
1                      3
2
2
2
3
5

分析：
暴力算法：只需再创建一个数组以记录原数组中每一个数出现的个数。遍历原数组，若某个数重复出现，则新数组++。最后找出新数组中最大的数和其对应的原数组中的数即可。
分治算法：
1、将原数组排序，找出中位数假定为众数，并且确定中位数的个数，记录左右两边的下标。这样，把数组分为了左右两部分。
2、如果左边部分的个数大于中位数的个数，则遍历左边，右边同理。
3、若左边剩余元素还大于重数，则有可能会存在更大的众数重数，递归寻找。
4、若右边剩余元素还大于重数，则有可能会存在更大的众数重数，递归寻找。

代码如下

#include <bits/stdc++.h>
using namespace std;
int maxsum=0;
int maxzhong=0;
void fun(int a[],int l,int r)
{
    int mid=(r+l)/2;
    int i,j;
    i=mid-1;
    j=mid+1;
    int sum=1;
    while(a[mid]==a[i]&&i>=l)
    {
        sum++;
        i--;
    }
    while(a[mid]==a[j]&&j<=r)
    {
        sum++;
        j++;
    }
    if(sum==maxsum&&a[mid]<maxzhong)
    {
        maxsum=sum;
        maxzhong=a[mid];
    }
    if(sum>maxsum)
    {
        maxsum=sum;
        maxzhong=a[mid];
    }
    if(i-l+1>maxsum)
    {
        fun(a,l,i);
    }
    if(r-j+1>maxsum)
    {
        fun(a,j,r);
    }
}
int main()
{
    int n;
    cin>>n;
    int a[1000];
    for(int i=0;i<n;i++)
        cin>>a[i];
    sort(a,a+n);
    fun(a,0,n-1);
    cout<<maxzhong<<endl;
    cout<<maxsum;
    return 0;
}

二、邮局选址问题

问题描述：
　　在一个按照东西和南北方向划分成规整街区的城市里，n 个居民点散乱地分布在不同的街区中。用x 坐标表示东西向，用y 坐标表示南北向。各居民点的位置可以由坐标(x,y) 表示。街区中任意2 点(x1,y1) 和(x2,y2) 之间的距离可以用数值|x1-x2|+|y1-y2| 度量。
居民们希望在城市中选择建立邮局的最佳位置，使n 个居民点到邮局的距离总和最小。
编程任务：
　　给定n 个居民点的位置,编程计算n 个居民点到邮局的距离总和的最小值。
数据输入：
　　由文件input.txt 提供输入数据。文件的第1 行是居民点数n，1<=n<=10000。接下来n 行是居民点的位置，每行2 个整数x 和y，-10000<=x，y<=10000。
结果输出:
　　程序运行结束时，将计算结果输出到文件output.txt 中。文件的第1 行中的数是n 个居民点到邮局的距离总和的最小值。
输入文件示例               输出文件示例
input.txt                  output.txt
5                          10
1 2
2 2
1 3
3 -2
3 3

分析：
分治法：
1、将输入的坐标（x,y）分别排序（快排）。
2、分别取两个数组的中位数，数组的每一个元素减去中位数的绝对值相加即可得到最短距离。
代码如下

#include <bits/stdc++.h>
using namespace std;
void quicksort(int a[],int l,int r)
{
    int temp=a[l];
    int t;
    if(l>r) return;
    int i=l,j=r;
    while(i!=j)
    {
        while(a[j]>=temp&&i<j) j--;
        while(a[i]<=temp&&i<j) i++;
        if(i<j)
        {
            t=a[i];
            a[i]=a[j];
            a[j]=t;
        }
    }
    a[l]=a[i];
    a[i]=temp;
    quicksort(a,l,i-1);
    quicksort(a,i+1,r);
    return ;
}
void getdistance(int x[],int y[],int n)
{
    int Distance=0;
    quicksort(x,0,n-1);
    quicksort(y,0,n-1);
    int dx=x[n/2];
    int dy=y[n/2];
    for(int i=0;i<n;i++)
        Distance += abs(x[i]-dx)+abs(y[i]-dy);
        cout<<Distance<<endl;
}
int main()
{
    int n;
    cin>>n;
    int x[100],y[100];
    for(int i=0;i<n;i++)
        cin>>x[i]>>y[i];
    getdistance(x,y,n);
}

三、输油管道问题

问题描述：
　　某石油公司计划建造一条由东向西的主输油管道。该管道要穿过一个有n 口油井的油田。从每口油井都要有一条输油管道沿最短路经(或南或北)与主管道相连。如果给定n 口油井的位置,即它们的x 坐标（东西向）和y 坐标（南北向）,应如何确定主管道的最优位置, 即使各油井到主管道之间的输油管道长度总和最小的位置?证明可在线性时间内确定主管道的最优位置。
编程任务：
　　给定n 口油井的位置,编程计算各油井到主管道之间的输油管道最小长度总和。
数据输入：
　　由文件input.txt 提供输入数据。文件的第1 行是油井数n，1<=n10000。接下来n 行是油井的位置，每行2 个整数x 和y，-10000<=x，y<=10000 。
结果输出:
　　程序运行结束时，将计算结果输出到文件output.txt 中。文件的第1 行中的数是油井到主管道之间的输油管道最小长度总和。
输入文件示例       输出文件示例
input.txt           output.txt
5                       6
1 2
2 2
1 3
3 -2
3 3

分析：
注意此题只需考虑y。
代码如下

#include <bits/stdc++.h>
using namespace std;
void quicksort(int a[],int l,int r)
{
    int temp=a[l];
    int t;
    if(l>r) return;
    int i=l,j=r;
    while(i!=j)
    {
        while(a[j]>=temp&&i<j) j--;
        while(a[i]<=temp&&i<j) i++;
        if(i<j)
        {
            t=a[i];
            a[i]=a[j];
            a[j]=t;
        }
    }
    a[l]=a[i];
    a[i]=temp;
    quicksort(a,l,i-1);
    quicksort(a,i+1,r);
    return ;
}
void getdistance(int y[],int n)
{
    int Distance=0;
    quicksort(y,0,n-1);
    int dy=y[n/2];
    for(int i=0;i<n;i++)
        Distance += abs(y[i]-dy);
        cout<<Distance<<endl;
}
int main()
{
    int n;
    cin>>n;
    int x[100],y[100];
    for(int i=0;i<n;i++)
        cin>>x[i]>>y[i];
    getdistance(y,n);
}

四、集合划分问题

问题描述：
　　n 个元素的集合{1,2,., n }可以划分为若干个非空子集。例如，当n=4 时，集合{1，2， 3，4}可以划分为15 个不同的非空子集如下：
{{1}，{2}，{3}，{4}}，
{{1，2}，{3}，{4}}，
{{1，3}，{2}，{4}}，
{{1，4}，{2}，{3}}，
{{2，3}，{1}，{4}}，
{{2，4}，{1}，{3}}，
{{3，4}，{1}，{2}}，
{{1，2}，{3，4}}，
{{1，3}，{2，4}}，
{{1，4}，{2，3}}，
{{1，2，3}，{4}}，
{{1，2，4}，{3}}，
{{1，3，4}，{2}}，
{{2，3，4}，{1}}，
{{1，2，3，4}}
其中，集合{{1，2，3，4}} 由1 个子集组成；集合{{1，2}，{3，4}}，{{1，3}，{2， 4}}，{{1，4}，{2，3}}，{{1，2，3}，{4}}，{{1，2，4}，{3}}，{{1，3，4}，{2}}，{{2， 3，4}，{1}} 由2 个子集组成；集合{{1，2}，{3}，{4}}，{{1，3}，{2}，{4}}，{{1，4}， {2}，{3}}，{{2，3}，{1}，{4}}，{{2，4}，{1}，{3}}，{{3，4}，{1}，{2}} 由3 个子集组成；集合{{1}，{2}，{3}，{4}} 由4 个子集组成。
编程任务：
　　给定正整数n 和m，计算出n 个元素的集合{1,2,., n }可以划分为多少个不同的由m 个非空子集组成的集合。
数据输入：
　　由文件input.txt 提供输入数据。文件的第1 行是元素个数n 和非空子集数m。
结果输出:
　　程序运行结束时，将计算出的不同的由m个非空子集组成的集合数输出到文件output.txt 中。
输入文件示例            输出文件示例
input.txt               output.txt
4    3                      6

分析：
如果n==m或m=1，则答案为1。

4个元素分成3个非空子集组成的集合，分为两步：
1、3个元素分成2个非空子集（3）

{{1，2}，{3}}，
{{1，3}，{2}}，
{{2，3}，{1}}，

2、3个元素分为3个子集

{{1}，{2}，{3}}

3、向1、2步骤的子集中插入4这个元素

{{1，2}，{3}，{4}}，
{{1，3}，{2}，{4}}，
{{2，3}，{1}，{4}}，//个数与1重合

{{1，4}，{2}，{3}}，
{{2，4}，{1}，{3}}，
{{3，4}，{1}，{2}}，//个数与3*2重合

综上所述，其实F（4,3）=F（3,2）+3*F（4,2）
归纳为F(n,m)=F(n-1,m-1)+m*F(n,m-1)
代码如下

#include<bits/stdc++.h>
using namespace std;
int fun(int n,int m)
{
    if(m==n||m==1)
        return 1;
    else
        return fun(n-1,m-1)+m*fun(n-1,m);
}
int main()
{
    int n,m;
    cin>>n>>m;
    cout<<fun(n,m);
    return 0;
}

五、整数因子分解问题

问题描述：
 　大于1 的正整数n 可以分解为：n=x1*x2*…*xm。
例如，当n=12 时，共有8 种不同的分解式：
12=12；
12=6*2；
12=4*3；
12=3*4；
12=3*2*2；
12=2*6；
12=2*3*2；
12=2*2*3 。
编程任务：
　　对于给定的正整数n，编程计算n 共有多少种不同的分解式。
数据输入：
　　由文件input.txt 给出输入数据。第一行有1 个正整数n (1≤n≤2000000000)。
结果输出:
将计算出的不同的分解式数输出到文件output.txt 。
输入文件示例          输出文件示例
input.txt            output.txt
　12                      8

分析：
对数n从2到本身n依次取余，若对i取余数结果为0，则i为n的整数因子，对n/i继续递归取余，直至余数为1。

代码如下

#include<bits/stdc++.h>
using namespace std;
int sum=0;
void fun(int n)
{
    int i=2;
    if (n==1) sum++;
    else{
        while(i<=n)
        {
            if(n%i==0)
            {
                fun(n/i);
            }
            i++;
        }
    }
}
int main()
{
    int n;
    cin>>n;
    fun(n);
    cout<<sum;
}