Flight Boarding Optimization

题目链接

• 题意：
  n个数字，范围是1-s，需要将1-s分成k段，使得每个数字必然只属于一段。分别计算每段的正序对，求和即为答案，现在需要让这个和最小，求最小值
• 分析：
  既然有k段这个限制，那么采用DP的时候必然需要将k作为一个状态：所以状态表示为dp[i][j]，以i为分段的最后一个数字，一共j段的答案。因为状态转移一定是dp[i][j] -> dp[k][j + 1]，所以将第二维优化。
  状态转移的时候，需要求一下在[i + 1, j]区间内的正序对个数，这个可以预处理解决：对于一个正序对x - y，在二维坐标中将（x，y）值加一。最后求得时候，其实就是要x和y均在给定的区间内的正序对个数，就是在一个正方形内的点的个数，预处理出每个点到原点的矩形中点的个数即可。
• tips：
  这题目本来是在湖南大学OJ上做的，不加离散化会超时。。。表示无奈。找到的原OJ是正常的，不加离散化更快
  

加离散化（最差情况会增加常数）原OJ 202ms：

const int MAXN = 1010;

int S[MAXN][MAXN];
int ipt[MAXN];
int dp[MAXN];

inline int query(int l, int r)
{
    return S[r][r] - S[l - 1][r] - S[r][l - 1] + S[l - 1][l - 1];
}

int xx[MAXN], tot = 0, to[MAXN];
int main()
{
    tot = 1;
    int n, Max, k;
    RIII(n, Max, k);
    REP(i, n)
    {
        RI(ipt[i]);
        xx[tot++] = ipt[i];
    }
    sort(xx, xx + tot);
    int idx = 1;
    FF(i, 1, tot)
        if (xx[i] != xx[i - 1])
            xx[idx++] = xx[i];
    tot = idx;
    REP(i, tot)
        to[xx[i]] = i;
    Max = tot;
    REP(i, n)
        ipt[i] = to[ipt[i]];

    for (int i = 0; i <= Max; i++)
        dp[i] = INF;
    dp[0] = 0;

    REP(i, n) FF(j, i + 1, n)
        if (ipt[i] < ipt[j])
            S[ipt[i]][ipt[j]]++;
    FE(i, 1, Max) FE(j, 1, Max)
        S[i][j] += S[i - 1][j] + S[i][j - 1] - S[i - 1][j - 1];

    REP(kase, k)
    {
        FED(i, Max, 0)
        {
            for (int j = Max; j >= i + 1; j--)
            {
                int nxt = dp[i] + query(i + 1, j);
                if (nxt < dp[j])
                    dp[j] = nxt;
            }
        }
    }
    WI(dp[Max]);
    return 0;
}

不加离散化（湖大OJ超时，原OJ 187ms）：

const int MAXN = 1010;

int S[MAXN][MAXN];
int ipt[MAXN];
int dp[MAXN];

inline int query(int l, int r)
{
   return S[r][r] - S[l - 1][r] - S[r][l - 1] + S[l - 1][l - 1];
}

int main()
{
   int n, Max, k;
   scanf("%d%d%d", &n, &Max, &k);
   REP(i, n)
       scanf("%d", &ipt[i]);
   for (int i = 0; i <= Max; i++)
       dp[i] = INF;
   dp[0] = 0;
   REP(i, n) FF(j, i + 1, n)
       if (ipt[i] < ipt[j])
           S[ipt[i]][ipt[j]]++;
   FE(i, 1, Max) FE(j, 1, Max)
       S[i][j] += S[i - 1][j] + S[i][j - 1] - S[i - 1][j - 1];

   REP(kase, k)
   {
       FED(i, Max, 0)
       {
           for (int j = Max; j >= i + 1; j--)
           {
               int nxt = dp[i] + query(i + 1, j);
               if (nxt < dp[j]) dp[j] = nxt;
           }
       }
   }
   WI(dp[Max]);
   return 0;
}