DP:Miking Time(POJ 3616)

　　　　　　　　　　　　　　2015-09-21

　　　　　　　　　　　　　　


　　　　　　　　　　　　　　　奶牛挤奶

　　题目大意就是这只Bessie的牛产奶很勤奋，某农民有一个时刻表，在N时间内分成M个时间段，每个时间段Bessie会一直产奶，然后有一定的效益，并且Bessie产奶后要休息两个小时。

　　这是一道很简答的DP，因为区间是不重复的，所以我们只要快排一下就好了，然后从第一个时间段到最后一个时间段

　　　　状态转移方程：

　　　　　　dp[i][j]=dp[i-1][j] j<i;

　　　　　　dp[i][i]=MAX(dp[i][i],dp[i][j]+list[i].eff); j<i

　　　　很简单是吧，我这道题做了两个小时？？？？

　　　　为什么？？？因为我快排写错了！


　　　　代码一开始交的时候是1001*1001矩阵，时间是0ms，但是内存要用到7000+，用滚动数组时间到16ms，但是内存变成了132，还是值得的

　　

#include <stdio.h>
#include <stdlib.h>
#include <string.h>
#define MAX(a,b) (a)>(b)?(a):(b)
#define CUTOFF 20

typedef int Position;
typedef struct input
{
int start_hour;
int end_hour;
int eff;
}LIST;

LIST list[1000];
long long dp1[1001];
long long dp2[1001];

void Quick_Sort(Position, Position);
void Swap(Position, Position);
void Insertion_Sort(Position, Position);
int Median_Of_Three(Position, Position, Position);
void Search(const int, const int, const int);

int main(void)
{
int interval, R, N, i;

while (~scanf("%d%d%d", &N, &interval, &R))
{
for (i = 1; i <= interval; i++)
scanf("%d%d%d", &list[i].start_hour, &list[i].end_hour, &list[i].eff);
list[0].start_hour = INT_MIN; list[0].end_hour = -R;
Quick_Sort(1, interval);
Search(N, interval, R);
}
return 0;
}

int Median_Of_Three(Position left, Position middle, Position right)
{
if (list[left].start_hour > list[middle].start_hour)
Swap(left, middle);
if (list[left].start_hour > list[right].start_hour)
Swap(left, right);
if (list[middle].start_hour > list[right].start_hour)
Swap(middle, right);
Swap(middle, right);//隐藏枢纽元
return list[right].start_hour;
}

void Swap(Position x, Position y)
{
list[x].start_hour ^= list[y].start_hour;
list[y].start_hour ^= list[x].start_hour;
list[x].start_hour ^= list[y].start_hour;

list[x].end_hour ^= list[y].end_hour;
list[y].end_hour ^= list[x].end_hour;
list[x].end_hour ^= list[y].end_hour;

list[x].eff ^= list[y].eff;
list[y].eff ^= list[x].eff;
list[x].eff ^= list[y].eff;
}

void Insertion_Sort(Position left, Position right)
{
Position i, j;
int tmp_s, tmp_e, tmp_eff;
for (i = left + 1; i <= right; i++)
{
tmp_s = list[i].start_hour;
tmp_e = list[i].end_hour;
tmp_eff = list[i].eff;
for (j = i; j > left && list[j - 1].start_hour > tmp_s; j--)
{
list[j].start_hour = list[j - 1].start_hour;
list[j].end_hour = list[j - 1].end_hour;
list[j].eff = list[j - 1].eff;
}
list[j].start_hour = tmp_s;
list[j].end_hour = tmp_e;
list[j].eff = tmp_eff;
}
}

void Quick_Sort(Position left, Position right)
{
Position mid = (left + right) / 2, i, j;
int pivot;

if (right - left > CUTOFF)
{
pivot = Median_Of_Three(left, mid, right);
i = left; j = right;
while (1)
{
while (list[++i].start_hour < pivot);
while (list[--j].start_hour > pivot);
if (i < j)
Swap(i, j);
else break;
}
Swap(i, right);
Quick_Sort(left, i - 1);
Quick_Sort(i + 1, right);
}
else Insertion_Sort(left, right);
}

void Search(const int N, const int interval, const int R)
{
int i, j;
long long ans = -1;
long long *now = dp2, *prev = dp1, *tmp = NULL;

for (i = 1; i <= interval; i++)
{
for (j = i - 1; j >= 0; j--)
{
now[j] = prev[j];
if (list[i].start_hour - list[j].end_hour >= R)
now[i] = MAX(now[i], now[j] + list[i].eff);
ans = MAX(ans, now[i]);
}
tmp = now; now = prev; prev = tmp;
}
printf("%lld\n", ans);
}