DPトピックのCF 4D

質問意図されている：所与大きい矩形より最小矩形2. 1の幅と高スト厳密に単調増加、単調に増加するシーケンス：長方形のサイズには現在n番目の矩形の態様、満足配列要件の組成を決定しました。
その数のQ最長のシーケンス。
やや似LISので、我々はDPを選択するn個の印刷結果の利便性のためのアルゴリズムを乗。
（DP順序DP [i]は、次に[I] = maxのDP、iはサイズシーケンスの終わりまで上昇させる | DP [J] <DP [I] [I]、DP [J] +1）
、すなわち現在の缶DPそれはへの移籍の全てよりも少ないです。次いで、前項のレコード別のプレ配列を開きます。
<=次いで、従来の直接書き込み、第二とを比較する第1比較器によれば、これは一般的なサイズである（Waが、カスタムオペレータのうち、すなわち、w及びhは、<=、それが使用されるべきであることはできない|| ）
PS：この質問は、前半は8太行HHHHを見つけ書いた,,それは水問題のシミュレーションだと思いました

#include<iostream>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<cstdlib>
#include<climits>
#include<stack>
#include<vector>
#include<queue>
#include<set>
#include<map>
//#include<regex>
#include<cstdio>
#define up(i,a,b)  for(int i=a;i<b;i++)
#define dw(i,a,b)  for(int i=a;i>b;i--)
#define upd(i,a,b) for(int i=a;i<=b;i++)
#define dwd(i,a,b) for(int i=a;i>=b;i--)
//#define local
typedef long long ll;
const double esp = 1e-6;
const double pi = acos(-1.0);
const int INF = 0x3f3f3f3f;
const int inf = 1e9;
using namespace std;
int read()
{
    char ch = getchar(); int x = 0, f = 1;
    while (ch<'0' || ch>'9') { if (ch == '-')f = -1; ch = getchar(); }
    while (ch >= '0' && ch <= '9') { x = x * 10 + ch - '0'; ch = getchar(); }
    return x * f;
}
typedef pair<int, int> pir;
#define lson l,mid,root<<1
#define rson mid+1,r,root<<1|1
#define lrt root<<1
#define rrt root<<1|1
int n, w, h;
const int N = 5005;
int pre[N];
int dp[N];
struct node { int w, h, pos; bool operator<(const node a)const { return w == a.w ? h < a.h : w < a.w; }
bool operator<=(const node a)const { return w <= a.w||h <= a.h; }
}vec[N];
int main()
{
    n = read(), w = read(), h = read();
    int x, y;
    up(i, 0, n) { x = read(), y = read(); vec[i] = node{ x,y,i }; }
    sort(vec, vec + n);
    memset(pre, -1, sizeof(pre));
    up(i, 0, n)
    {
        if (vec[i] <= node{ w,h,0 })continue;
        dp[i] = 1;
        up(j, 0, i)
        {
            if (vec[i] <= vec[j])continue;
            if (vec[j] <= node{ w,h,0 })continue;
            if (dp[i] < dp[j] + 1)
            {
                dp[i] = dp[j] + 1;
                pre[i] = j;
            }
        }
    }
    int ans = -1; int pos =-1;
    up(i, 0, n)
    {
        if (ans < dp[i])ans = dp[i], pos = i;
    }
    if (ans == 0) { cout << "0" << endl; }
    else
    {
        vector<int>ansvec;
        for (int i = pos; ~i; i = pre[i])
        {
            ansvec.push_back(i);
        }
        reverse(ansvec.begin(), ansvec.end());
        cout << ansvec.size() << endl;
        for (auto k : ansvec)cout << vec[k].pos+1 << " ";
    }
    return 0;
}