活动安排问题 动态规划、贪心算法C语言实现

问题描述

  设有n个活动的集合E={1,2,…,n}，其中每个活动都要求使用同一资源，如演讲会场等，而在同一时间内只有一个活动能使用这一资源。每个活动i都有一个要求使用该资源的起始时间si和一个结束时间fi,且si <fi 。如果选择了活动i，则它在半开时间区间[si, fi)内占用资源。
  若区间[si, fi)与区间[sj, fj)不相交，则称活动i与活动j是相容的。也就是说，当si≥fj或sj≥fi时，活动i与活动j相容。活动安排问题就是要在所给的活动集合中选出最大的相容活动子集合。

思路

  其实就是一个不相交的子集合，妥妥的一个经典贪心问题。做法其实挺多的，个人感觉精髓还是寻找贪心的点：你可以贪心时长，越短越好；也可以贪心结束时间或者开始时间之类的。我们这里采取贪心结束时间和时长，就是越早结束越好，如果一起结束，那就选时间最短的。
  步骤如下：

• Step1 以结束时间为主要关键字，开始时间为次要关键字排序
• Step2 直接把第一个定为起始活动
• Step3 遍历并判断与已经在队列中的数值是否相容

  经典贪心在此是可以得到最优解的，只需证明：

1.  对于排序后的序列，一定存在一种最优解包括第一个活动（最早结束，时间最短）；
2.  说明活动选择问题具有优化子结构；
3.  算法按照Greedy 选择性计算解；

  这不是我们实现的重点，在此不赘述。

  做到这里，我们发现，最后的贪心过程跟01背包处的的dp[i]=max(dp[i], dp[i-a[j]]+1)有点像，因此联系到动态规划的思想。当然我认为这并不是严格意义上的动态规划，起码并不符合《背包九讲》的说法，但是问题不大，仅仅是一个尝试而已，写着玩的。

贪心实现

#include<bits/stdc++.h>
using namespace std;

struct activity{
	int start, end;
};

bool cmp(const activity &a,const activity &b)
{
	if(a.end!=b.end)
		return a.end<b.end;
	else
		return a.start>b.start;
}

int main(){
	int n;
	scanf("%d", &n);
	int satrt[n], end[n], dp[n];
	struct activity a[n];
	for(int i=0;i<n;i++)
		scanf("%d",&a[i].start);
	for(int i=0;i<n;i++)
		scanf("%d",&a[i].end);
	sort(a,a+n,cmp);
	for (int i = 0;i<n;i++){
		printf("%d\t%d\n",a[i].start, a[i].end);
	}
	int now_end = a[0].end;
	int len = 1;
	for(int i=1;i<n;i++){
		if (a[i].start>now_end){
			now_end = a[i].end;
			len++;
		}
		printf("%d\t", now_end);
	}
	printf("%d", len);
	return 0;
} 

//测试数据：
//11
//1 3 0 5 3 5 12 2 8 8 6
//4 5 6 7 7 9 14 13 12 11 10 

动态规划实现

#include<bits/stdc++.h>
using namespace std;

struct activity{
	int start, end;
};

bool cmp(const activity &a,const activity &b)
{
	if(a.end!=b.end)
		return a.end<b.end;
	else
		return a.start>b.start;
}


int main(){
	int n;
	scanf("%d", &n);
	int satrt[n], end[n], dp[n];
	struct activity a[n];
	for(int i=0;i<n;i++){
		scanf("%d",&a[i].start);
		scanf("%d",&a[i].end);
	}
	sort(a,a+n,cmp);
	for (int i = 0;i<n;i++){
		printf("%d\t%d\n",a[i].start, a[i].end);
	} 
	memset(dp,0,sizeof(dp));
	dp[0]=1;
	for(int i=1;i<n;i++){
		for(int k=0;k<i;k++)
			dp[i]=max(dp[k],dp[i]);
		for(int j=0;j<i;j++)
			if(a[j].end<=a[i].start)
				dp[i]=max(dp[i], dp[j]+1);
	}
	for(int i=0;i<n;i++)
		dp[n-1]=max(dp[n-1],dp[i]);
	printf("%d",dp[n-1]);
	return 0;
} 

//测试数据：
//11
//1 3 0 5 3 5 12 2 8 8 6
//4 5 6 7 7 9 14 13 12 11 10