topcoder srm 520 div1

problem1 link

设$f[i][j][k]$表示考虑了前$i$道题，剩下时间为$j$，剩下技能为$k$的最大得分．

从小到大计算二元组$(j,k)$的话，在存储上可以省略掉$i$这一维．

problem2 link

首先，不同的提交状态有８种．预计算每一种提交状态的每一个分值的种数，设为$g[mask][score]$．对于$g[mask][score]$的计算，可以枚举$mask$包含的一个题目，设为$i$,其分值为$p_{i}$，那么第$i$道题得分如果为$x$，那么其他$mask^{'}=mask$^$2^{i}$题目得分为$score-x$，所以$g[mask][score]=\sum_{x=1}^{p_{i}}g[mask^{'}][score-x]$

最后从前向后考虑每个人的得分．设$f[i][score]$表示前$i$个人的得分降序排列且第$i$个人得分为$score$的方案数．设第$i+1$个人的得分为$y$,由于第$i$个人的得分一定大于$y$，所以$f[i+1][y]=\sum_{x=y+1}^{MaxScore_{i}}f[i][x]$

由于计算数组$g$和$f$都是连续求和，可以通过预处理前缀和或者后缀和来加速计算

problem3 link

将所有人按照是否攻击别人和是否被攻击成功分为四类，即YY,YN,NY,NN．设这四类的个数分别为$A,B,C,D,A+B+C+D=n$．所有的攻击一共有$T=A+B$个．

明显如果$A+C>B+C$，那么答案为０．下面认为$C \le B$

现在就是要构造一个$T*2$的二维表，某一行的两列的两个值$r_{1},r_{2}$，表示有一条攻击是$r_{1}->r_{2}$

对于第一类中的每个元素$e$，它一定在第一列第二列同时出现过，并且可能在第二列出现过多次．现在只考虑出现在第二列且被攻击成功的那一行，设其为$y$,设出现的第一列的行为$x$，那么有$x<y$．

现在来计算把$A$个第一类元素放到这个表中的方案数．可以用动态规划来计算．设$f[T][A]$表示把$A$个元素放入到一个$T*2$表中的方案数．现在先假设所有的$A$个元素是一样的，即没有先后顺序．

考虑第一行的第一列，有两种情况：

(1)没有放置$A$中的元素，那么情况为$f[T-1][A]$

(2)放置了一个，那么这个元素出现在第二列的行数在后面的$T-1$行中，且这$T-1$行的第二列中已经放置了$A-1$个，所以还有$(T-1)-(A-1)=T-A$个位置，即答案为$(T-A)*f[T-1][A-1]$

所以$f[T][A]=f[T-1][A]+(T-A)*f[T-1][A-1]$

现在对于第二列放置的每个第一类的元素，这一行所代表的攻击是成功的．也就是说对应的第一列也已经有了攻击者（可能是第一类或者第二类的元素）．$f[T][A]$已经保证了对于第一类元素的合理性．

现在现在第一列对应的$T-A=B$个元素（这些元素可能是第二类也可能是第一类但还没有攻击对手）还没有攻击的对象，因为这$B$个攻击中有$C$个需要是成功的，所以有$C_{B}^{C}$个选择．现在第一列还有$B-C$个还没有攻击对象，并且这些攻击都是失败的，所以只要不是它自己就行，所以每个有$n-1$种．

最后所有的元素都是不同的，所以要乘以它们的阶乘．所以答案为$f[A+B][A]*C_{B}^{C}*(n-1)^{B-C}*A!*B!*C!$

code for problem1

#include <iostream>
#include <vector>
#include <algorithm>
using namespace std;

class SRMCodingPhase {
  public:
    int countScore(vector <int> points, vector <int> skills, int luck) {
      vector<vector<int>> f(76, vector<int>(luck + 1, -1));
      f[75][luck] = 0;
      const int b[] = {2, 4, 8};
      for (int i = 0; i < 3; ++ i) {
        const int point = points[i];
        const int skill = skills[i];
        for (int j = 0; j <= 75; ++ j) {
          for (int k = 0; k <= luck; ++ k) {
            if (f[j][k] == -1) {
              continue;
            }
            for (int t = 0; t <= k && t < skill; ++ t) {
              const int last_point = point - b[i] * (skill - t);
              if (last_point <= 0) {
                continue;
              }
              if (skill - t > j) {
                continue;
              }
              f[j - (skill - t)][k - t] = max(f[j - (skill - t)][k - t], last_point + f[j][k]);
            }
          }
        }
      }
      int result = 0;
      for (int j = 0; j <= 75; ++ j) {
        for (int k = 0; k <= luck; ++ k) {
          result = max(result, f[j][k]);
        }
      }
      return result;
    }
};

　　

code for problem2

#include <iostream>
#include <vector>
#include <algorithm>
#include <string>
using namespace std;

#define mod 1000000007
#define MAX_SCORE 200000

int g[8][MAX_SCORE + 1];
int f[20][MAX_SCORE + 1];

class SRMIntermissionPhase {
  public:
    int countWays(vector<int> points, vector<string> description) {
      Initialize(points);

      const int mask0 = GetMask(description[0]);
      for (int i = 0; i <= MAX_SCORE; ++ i) {
        f[0][i] = GetRangeSum(mask0, i, i);
      }
      CalculateSuffixSum(0);
      const int person_number = (int)description.size();
      for (int i = 1; i < person_number; ++ i) {
        const int mask = GetMask(description[i]);
        for (int j = 0; j < MAX_SCORE; ++ j) {
          f[i][j] = (long long)GetRangeSum(mask, j, j) * f[i - 1][j + 1] % mod;
        }
        CalculateSuffixSum(i);
      }
      return f[person_number - 1][0];
    }
  private:

    void CalculateSuffixSum(int idx) {
      for (int i = MAX_SCORE - 1; i >= 0; -- i) {
        Add(f[idx][i], f[idx][i + 1]);
      }
    }
    void CalculatePrefixSum(int idx) {
      for (int i = 1; i <= MAX_SCORE; ++ i) {
        Add(g[idx][i], g[idx][i - 1]);
      }
    }

    int GetMask(const std::string& s) {
      int mask = 0;
      for (int i = 0; i < 3; ++ i) {
        if (s[i] == 'Y') {
          mask |= 1 << i;
        }
      }
      return mask;
    }

    void Initialize(const vector<int>& points) {
      g[0][0] = 1;
      CalculatePrefixSum(0);
      for (int i = 1; i < 8; ++ i) {
        int low_bit = 0;
        for (int j = 0; j < 3; ++ j) {
          if ((i & (1 << j)) != 0) {
            low_bit = j;
            break;
          }
        }
        const int other = i ^ (1 << low_bit);
        const int point = points[low_bit];
        for (int k = 1; k <= MAX_SCORE; ++ k) {
          g[i][k] = GetRangeSum(other, k - point, k - 1);
        }
        CalculatePrefixSum(i);
      }
    }

    int GetRangeSum(int mask, int left, int right) {
      if (left <= 0) {
        return g[mask][right];
      }
      int result = g[mask][right] - g[mask][left - 1];
      if (result < 0) {
        result += mod;
      }
      return result;
    }
    void Add(int &x, int y) {
      x += y;
      if (x >= mod) {
        x -= mod;
      }
    }
};

code for problem3

#include <cstring>
#include <iostream>
#include <map>
#include <sstream>
#include <stdio.h>
#include <string>
#include <vector>
using namespace std;

#define mod 1000000007
#define MAX_SIZE 2501

int f[MAX_SIZE][MAX_SIZE];

class SRMChallengePhase {
  public:
  int countWays(vector<string> codersAttempted, vector<string> codersChallenged)
  {
    const std::string s1 = Concat(codersAttempted);
    const std::string s2 = Concat(codersChallenged);
    const int n = (int)s1.size();
    int A = 0, B = 0, C = 0;
    for (int i = 0; i < n; ++i) {
      if (s1[i] == 'Y' && s2[i] == 'Y') {
        ++A;
      } else if (s1[i] == 'Y' && s2[i] == 'N') {
        ++B;
      } else if (s1[i] == 'N' && s2[i] == 'Y') {
        ++C;
      }
    }
    if (B < C) {
      return 0;
    }
    f[0][0] = 1;
    for (int i = 1; i <= A + B; ++i) {
      for (int j = 0; j <= A; ++j) {
        f[i][j] = f[i - 1][j];
        if (j > 0) {
          Add(f[i][j], (long long)(i - j) * f[i - 1][j - 1] % mod);
        }
      }
    }
    int result = f[A + B][A];
    for (int i = 1; i <= B - C; ++i) {
      result = (long long)result * (n - 1) % mod;
    }
    for (int i = 2; i <= A; ++i) {
      result = (long long)result * i % mod;
    }
    for (int i = B; i >= 2; --i) {
      result = (long long)result * i % mod;
      if (i >= B - C + 1) {
        result = (long long)result * i % mod;
      }
    }
    return result;
  }

  private:
  void Add(int& x, int y)
  {
    if (y < 0) {
      y += mod;
    }
    x += y;
    if (x >= mod) {
      x -= mod;
    }
  }
  std::string Concat(const vector<string>& str)
  {
    std::string s = "";
    for (auto x : str) {
      s += x;
    }
    return s;
  }
};