POJ1065+POJ3636

由离散数学内容可知：任何集合A上的恒等关系，集合的幂集P（A）上的包含关系，实数集上的小于等于关系，正整数集上的整除关系等都是偏序关系。

而1065和3636两题就是要利用集合的包含关系来解决。

POJ1065：有n根木棒，每一根都具有各自的长度和重量，现如今要用机器需要处理这些木棒，如果第i+1根木棒的重量w’和长度l’和第i根木棒的w和l满足l<=l’&& w<=w’，那么机器将不会耗费时间，否则机器需要耗费一个单位的时间，要求处理完所有木棒需要的最少时间（包括机器启动的时间）。

POJ3636：套娃娃，只有当 w1 < w2 && h1 < h2 时才能将娃娃(w1,h1)套在娃娃(w2,h2)里面，问最终可以产生几组娃娃。

分析：两题之间的不同在于1065要求满足<=关系，而3636要求满足<关系。这点的不同导致了二者之间的排序方式的不同，两题中第一个参数都是从小到大排列，但是第二个参数的排列就有所区别了。前者要求是包含＝，所以第二个参数也是从小到大排；而后者则不包含 ＝ 所以要从大到小排 。因为w升序，h降序可以保证w相等时，一定不会出现覆盖的情形。

POJ1065代码：
O(nlogn)算法

#include<iostream>
#include<algorithm>
using namespace std;
const int Max = 5005;
struct Stick {
    int l, w;
}stick[Max];

int stk[Max];
bool cmp(Stick a, Stick b)
{
    return (a.l == b.l&&a.w < b.w) || a.l < b.l;
}

int main()
{
    int t, m;
    cin >> t;
    while (t--)
    {
        cin >> m;
        for (int i = 1; i <= m; i++)
            cin >> stick[i].l >> stick[i].w;

        sort(stick + 1, stick + m + 1, cmp);

        int top = 0;
        stk[top] = 10000000;
        for (int i = 1; i <= m; i++)
        {
            //最长递减子序列
            if (stick[i].w < stk[top]) {
                stk[++top] = stick[i].w;
            }
            else {//在stack中找第一个小于等于n[i].l的数替换
                int left = 1, right = top;
                while (left < right) {
                    int mid = (left + right) >> 1;
                    if (stk[mid] <= stick[i].w) {
                        right = mid;//stack[mid] == n[i].w时，right=mid保证不会越过stack中与n[i].w相等的数
                    }
                    else {
                        left = mid + 1;//stack[mid]>n[i].w 跳过stack[mid]到比stack[mid]小的数
                    }
                }//其实最终跳出的条件是left==right
                stk[right] = stick[i].w;
            }
        }
        cout << top << endl;
    }
    return 0;
}

O(n^2)算法：

#include<iostream>
#include<algorithm>

using namespace std;

const int Max = 5005;

struct Stick {
    int l, w;
}stick[Max];

int stk[Max];
bool cmp(Stick a, Stick b)
{
    return (a.l == b.l&&a.w < b.w) || a.l < b.l;
}

int main()
{
    int t, m;
    cin >> t;
    while (t--)
    {
        cin >> m;
        for (int i = 0; i < m; i++)
            cin >> stick[i].l >> stick[i].w;

        sort(stick, stick + m, cmp);

        int num = 0;
        for (int i = 0; i < m; i++)
        {
            int flag = 0;
            for (int j = 0; j < num; j++) {
                if (stick[i].w >= stk[j])//找到stk中第一个小于或等于n[i].w的数将其替换成n[i].w
                {
                    stk[j] = stick[i].w;
                    flag = 1;
                    break;
                }
            }
            if (!flag) { //n[i].w比stk中任何数都小则将其添加到stk中
                stk[num++] = stick[i].w;
            }
        }
        cout << num << endl;
    }
    return 0;
}

O(n^2)DP:

#include <iostream>
#include <algorithm>

using namespace std;

const int N = 5005;
const int inf = 0x3f3f3f3f;

struct Stick
{
    int l, w;
    bool operator < (const Stick &a)const
    {
        return (l == a.l&&w < a.w) || l < a.l;
    }
}stick[N];

int n, dp[N];

int main()
{
    int t;
    stick[0].w = inf;
    cin >> t;
    while (t--)
    {
        cin >> n;
        memset(dp, 0, sizeof(dp));
        for (int i = 1; i <= n; i++) {
            cin >> stick[i].l >> stick[i].w;
        }

        sort(stick + 1, stick + n + 1);

        for (int i = 1; i <= n; i++)
        {
            for (int j = 0; j < i; j++)
            {
                if (stick[j].w > stick[i].w)
                    dp[i] = max(dp[i], dp[j] + 1);
            }
        }

        int ans = 0;
        for (int i = 1; i <= n; i++) {
            ans = max(ans, dp[i]);
        }
        cout << ans << endl;
    }
    return 0;
}

POJ3636代码：

暴力贪心O(n^2)：

#include<iostream>
#include<algorithm>

using namespace std;

const int Max = 20005;

struct Doll {
    int w, h;
    bool operator<(const Doll& a) const {
        return w < a.w || (w == a.w&&h > a.h);
    }
}doll[Max];
Doll d[Max];

int main()
{
    int n, m;
    cin >> n;
    while (n--)
    {
        int result = 0;
        cin >> m;
        for (int i = 0; i < m; i++) {
            cin >> doll[i].w >> doll[i].h;
        }

        sort(doll, doll + m);

        for (int i = 0; i<m; i++)
        {
            int flag = 0;
            for (int j = 0; j<result; j++)
                if (doll[i].w>d[j].w&&doll[i].h>d[j].h)
                {
                    d[j] = doll[i];
                    flag = 1;
                    break;
                }
            if (!flag) {
                d[result++] = doll[i];
            }
        }
        cout << result << endl;
    }
    return 0;
}

O(nlogn):

#include<iostream>
#include<algorithm>

using namespace std;

const int Max = 20005;
const int inf=10000000;
struct Doll {
    int w, h;
    bool operator<(const Doll& a) const {
        return w < a.w || (w == a.w&&h > a.h);
    }
}doll[Max];
int d[Max];

int main()
{
    int n, m;
    cin >> n;
    while (n--)
    {
        cin >> m;
        for (int i = 1; i <= m; i++) {
            cin >> doll[i].w >> doll[i].h;
        }

        sort(doll+1, doll + m+1);

        int top = 0;
        d[top] = inf;
        for (int i = 1; i <= m; i++)
        {
            //最长非升子序列
            if (d[top] >= doll[i].h) {
                d[++top] = doll[i].h;
            }
            else {
                int low = 0, high = top;
                while (low <= high)
                {
                    int mid = (low + high) >> 1;
                    if (d[mid] < doll[i].h) {
                        high = mid - 1;
                    }
                    else {
                        low = mid + 1;
                    }
                }
                d[low] = doll[i].h;
            }
        }
        cout << top << endl;
    }
    return 0;
}

---

链（chain）是一个偏序集S的全序子集（所谓全序是指任意两个元素可比较）
反链（antichain）是一个偏序集S的子集，其中任意两个元素不可比较

偏序集的两个定理：
定理1 令（X,≤）是一个有限偏序集，并令r是其最大链的大小。则X可以被划分成r个但不能再少的反链。

其对偶定理称为Dilworth定理：
定理2 令（X,≤）是一个有限偏序集，并令m是反链的最大的大小。则X可以被划分成m个但不能再少的链。即：链的最少划分数=反链的最长长度。