实验二—分治和递归

h105l

已于 2024-01-10 18:16:23 修改

阅读量939

点赞数 15

文章标签：算法

于 2024-01-10 18:15:32 首次发布

本文链接：https://blog.csdn.net/m0_73201425/article/details/135510049

版权

文章介绍了如何使用分治和递归算法解决众数问题，涉及快速排序的应用，以及如何计算给定集合的子集划分和优化问题，如石油管道布局、邮局选址和整数因子分解。

摘要由CSDN通过智能技术生成

实验二—分治和递归

众数问题

问题描述：

给定含有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

算法思想：

因为输入是一个数组，所以可以考虑采用分治的算法解决这类简单问题

但是怎么来给每个元素计数呢？我们采用快速排序的思想将数组分为两部分，左边的<=基准数，右边的>=基准数。在分区的过程中便可以记录基准数的频次

递归的对左右两边进行相同的操作，得到左右两边的众数及频次

比较左右两边的众数和频次，返回题目要求值

要用到分治算法，如果不排序就进行分治的话，肯定会漏掉一些频数，所以快速排序的算法就派上用场了，使得分治所得结果是满足我们需求的

代码块

//快速排序划区
int pkey = data[low];
int l=low,h=high;
while(l<h)
{
  while(l<h && data[h]>=pkey){
    if(data[h]==pkey)  p[1]++;
    h--;
  }
  data[l]=data[h];
  while(l<h && data[l]<=pkey){
    if(data[l]==pkey)  p[1]++;
    l++;
  }
  data[h]=data[l];
}
data[l]=pkey;

//分治思想
//注意：函数中的数组和数组长度是不需要改变的，只需要改变起始位置即可
int ploc = l;
int *p2 = GetMode(data,data_length,low,ploc-1);
int *p3 = GetMode(data,data_length,ploc+1,high);

综合代码如下

#include <iostream>
using namespace std;

int* GetMode(int* data, int data_length, int low, int high) {
    int* p = new int[2]; // 创建一个数组来存储众数和频次
    p[0] = data[low]; // 设定初始值为当前数据中的第一个数
    p[1] = 1; // 设定初始频次为 1

    if (low < high) { // 当下界小于上界时，继续递归划分
        int pkey = data[low]; // 选择第一个数作为基准数
        int l = low, h = high;

        // 使用基准数进行分区
        while (l < h) {
            while (l < h && data[h] >= pkey) {
                if (data[h] == pkey) p[1]++; // 记录基准数的频次
                h--;
            }
            data[l] = data[h];

            while (l < h && data[l] <= pkey) {
                if (data[l] == pkey) p[1]++; // 记录基准数的频次
                l++;
            }
            data[h] = data[l];
        }
        data[l] = pkey;

        int ploc = l; // 基准数的位置
        int* p2 = GetMode(data, data_length, low, ploc - 1); // 左部分的众数
        int* p3 = GetMode(data, data_length, ploc + 1, high); // 右部分的众数

        if (p2[1] > p3[1])
        {p2 = p2;}
        else
        {p2 = p3;}
        
        if (p[1] > p2[1])
        { p = p;}
        else
        { p = p2;}
    }
    return p;
}

int main() {
    int n;
    cin >> n;
    int* arr = new int[n];
    for (int i = 0; i < n; i++) {
        cin >> arr[i];
    }
    int* p = GetMode(arr, n, 0, n - 1);

    cout << p[0] << "\n" << p[1] << "\n";
    delete[] arr;
    delete[] p;

    return 0;
}

集合划分问题

问题描述：

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 个子集组成；集合{{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-1个元素中划分m-1个子集：

countPartitions(n-1,m-1);

将一个元素加入到已有的子集中

那么剩下的n-1个元素需要划分成m个子集：

countPartitions(n-1,m);

考虑子问题的情况：

if n= =m || m= =1 :那么只有一种划分方式，返回1即可。
if m>n:无法划分集合，返回0。

#include <iostream>
using namespace std;

int funct(int n,int m)
{
    if(n == m||m == 1)
    {
        return 1;
    }
    else
    {
        return funct(n-1,m-1)+m*funct(n-1,m);
    }
}

int main()
{
    int n,m;
    cin>>n>>m;
    cout<<funct(n,m);
    return 0;
}

输油管道问题

问题描述：

某石油公司计划建造一条由东向西的主输油管道。该管道要穿过一个有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

算法思想

最重要的就是选出纵坐标中的中位数，那就是主管道的纵坐标。

然而选出中位数又可以使用分治的思想（对Select和Partion的结合使用）
Partion函数—分区操作
//和快速排序类似
int Partion(vector<int> &a,int l,int h){
  int x = a[l];
  while(l<h){
    while(l<h && a[h]>=x) h--;
    swap(a[l],a[h]);
    while(l<h && a[l]<=x) l++;
    swap(a[l],a[h]);
  }
  a[l]=x;
  return l;
}
Select函数—选择第K小的元素
int Select(vector<int> &a,int l,int h,int k){
  if l==h return a[l];
  int i=Partion(a,l,h);
  int j=i-l+1;
  if(k<=j) return Select(a,l,i,k);
  else return Select(a,i+1,h,k-j);
}

#include <iostream>
#include <algorithm>
#include <vector>
using namespace std;

void Swap(int &a, int &b) {
    int t = a;
    a = b;
    b = t;
}

int Partion(vector<int> &a, int l, int h) {
    int x = a[l];
    while (l < h) {
        while (l < h && a[h] >= x) --h;
        Swap(a[l], a[h]);
        while (l < h && a[l] <= x) ++l;
        Swap(a[l], a[h]);
    }
    a[l] = x;
    return l;
}

int Select(vector<int> &a, int l, int h, int k) {
    if (l == h) return a[l];
    int i = Partion(a, l, h);
    int j = i - l + 1;
    if (k <= j) return Select(a, l, i, k);
    else return Select(a, i + 1, h, k - j);
}

int main() {
    const int N = 10000;
    vector<int> a(N), b(N);
    int n, mid, sum = 0;
    cin >> n;
    for (int i = 0; i < n; i++)
        cin >> b[i] >> a[i];

    mid = Select(a, 0, n - 1, n / 2 + 1);
    for (int i = 0; i < n; i++)
        sum += abs(a[i] - mid);
    
    cout << sum;
    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

算法思想

和上一个题类似，我们只需要分别求出横坐标们和纵坐标们的中位数，如何与各个坐标做差求和即可解决问题

特别提出的是，我们这里求中位数时利用了已有的sort算法直接求出第n/2个数对应的值

#include<iostream>
#include<algorithm>
using namespace std;
int main()
{
	int n;
	int a[100],b[100];
	cin>>n;
	for(int i=0;i<n;i++)
	{
		cin>>a[i]>>b[i];
	}
	sort(a,a+n);
	sort(b,b+n);
	int sum=0,num=0;
	for(int i=0;i<n;i++)
	{
		sum+=abs(a[n/2]-a[i]);
		num+=abs(b[n/2]-b[i]);
	}
	cout<<sum+num<<endl;


}

整数因子分解问题

问题描述：

大于1 的正整数n 可以分解为：n=x1x2…*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

很明显本题需要用到递归的思想来解决问题，过程如下图所示：

#include <iostream>
using namespace std;
 
int total;
 
void solve(int n) {
    if (n == 1)
        total++;
    else {
        for (int i = 2; i <= n; i++) {
            if (n % i == 0)
                solve(n / i);
        }
    }
}
 
int main() {
    int n;
    cin>>n;
    solve(n);
    cout << total << endl;
    
    return 0;
}