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);
}