集合划分问题 分治

思路来源 海子 : http://www.cnblogs.com/dolphin0520/archive/2011/07/12/2103917.html

问题描述:
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}}

问题一:给出n,m,  输出按照此划分集合方法,  子集由m个集合组成的划分方式有几种, 如: -----inout: 4  3 ------output: 6

问题二:给出n, 输出可划分非空集合的划分方式的总数, 如: -----input: 4 -----output: 15

 

以上两个问题可以同时解决

---- f(n,m)表示将n个元素划分成由m个子集构成的集合的个数

case 1--------m=1时, f(n, m) = 1

case 2--------n=m时, f(n, m) = 1

case 3--------  f(n, m) = m*f(n-1, m) + f(n-1,m-1)

问题二答案就是f(n, m)  问题一就是cnt=f(n, 1) + f(n, 2) + f(n, 3) + ... + f(n, n)

#include <iostream>
#include <stdio.h>
using namespace std;

long long f[100][100];

void init()
{
	// initial
	for(int i = 1;i <= 99;i ++)
		f[i][1] = 1;

	for(int i = 2;i <= 50;i ++)
	{
		for(int j = 2;j <= i;j ++)
		{
			if(i == j)
				f[i][j] = 1;
			else
				f[i][j] = f[i-1][j-1] + (j)*f[i-1][j];
		}
	}
}

int main ()
{
	int n, m;
	init();
	//freopen("in.txt", "r", stdin);
	
	// 问题一
	while (scanf("%d %d", &n, &m) != EOF)
	{
		printf("%lld\n", f[n][m]);
	}
	// 问题二
	while (scanf("%d", &n) != EOF)
	{
		int cnt = 0;
		for(int i = 1;i <= n;i ++)
			cnt += f[n][i];
		printf("%lld\n", cnt);
	}
}

 

  • 0
    点赞
  • 7
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
邮局选址问题是指在一定范围内选取若干个邮局,使得所有居民到最近的邮局的距离之和最小。分治算法可以用于解决这个问题。 具体思路如下: 1. 将区域划分为若干个子区域,每个子区域内只有一个邮局。 2. 计算每个居民到最近的邮局的距离之和。 3. 选取距离之和最小的子区域,将其继续划分为子区域,直到每个子区域内只有一个居民点。 4. 在每个子区域内选取一个居民点作为邮局,计算所有居民到最近的邮局的距离之和。 5. 比较所有方案,选择距离之和最小的方案。 下面是使用C语言实现邮局选址问题分治算法的示例代码: ```c #include <stdio.h> #include <stdlib.h> #include <math.h> #define MAX 100 // 最大居民点数 // 定义居民点结构体 struct point { double x; // x坐标 double y; // y坐标 }; // 计算两点间的距离 double distance(struct point p1, struct point p2) { return sqrt(pow(p1.x - p2.x, 2) + pow(p1.y - p2.y, 2)); } // 计算所有居民到最近的邮局的距离之和 double cost(struct point* p, int n, struct point* q, int m) { double sum = 0; for (int i = 0; i < n; i++) { double min = distance(p[i], q[0]); for (int j = 1; j < m; j++) { double d = distance(p[i], q[j]); if (d < min) { min = d; } } sum += min; } return sum; } // 分治算法求解邮局选址问题 double solve(struct point* p, int n, struct point* q, int m) { if (m == 1) { // 只有一个邮局,直接计算距离之和 return cost(p, n, q, m); } else { // 将区域划分为两个子区域,递归求解 int k = n / 2; double x = p[k].x; double sum1 = solve(p, k, q, m); double sum2 = solve(p + k, n - k, q, m); double min = sum1 + sum2; int j = 0; for (int i = 0; i < n; i++) { if (p[i].x < x) { // 属于左边子区域 double d1 = cost(p, i, q, m); double d2 = solve(p + i, n - i, q, m); double sum = d1 + d2; if (sum < min) { min = sum; } } else { // 属于右边子区域 double d1 = cost(p + i, n - i, q, m); double d2 = solve(p, i, q, m); double sum = d1 + d2; if (sum < min) { min = sum; } q[j++] = p[i]; // 将右边子区域的居民点加入右边的邮局集合中 } } return min; } } int main() { int n, m; struct point p[MAX], q[MAX]; scanf("%d%d", &n, &m); for (int i = 0; i < n; i++) { scanf("%lf%lf", &p[i].x, &p[i].y); } printf("%.2lf\n", solve(p, n, q, m)); return 0; } ``` 注意在使用分治算法求解时,每次递归都需要开辟新的邮局集合q,以避免对原有集合的影响。

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值