bzoj 4742 [Usaco2016 Dec]Team Building

http://www.elijahqi.win/archives/3896
Description

Every year, Farmer John brings his NN cows to compete for “best in show” at the state fair. His arch
-rival, Farmer Paul, brings his MM cows to compete as well (1≤N≤1000,1≤M≤1000).Each of the N+MN+
M cows at the event receive an individual integer score. However, the final competition this year wi
ll be determined based on teams of KK cows (1≤K≤10), as follows: Farmer John and Farmer Paul both
select teams of KK of their respective cows to compete. The cows on these two teams are then paired
off: the highest-scoring cow on FJ’s team is paired with the highest-scoring cow on FP’s team, the s
econd-highest-scoring cow on FJ’s team is paired with the second-highest-scoring cow on FP’s team, a
nd so on. FJ wins if in each of these pairs, his cow has the higher score.Please help FJ count the n
umber of different ways he and FP can choose their teams such that FJ will win the contest. That is,
each distinct pair (set of KK cows for FJ, set of KK cows for FP) where FJ wins should be counted.
Print your answer modulo 1,000,000,009.
每年农夫约翰都会带着他的N只牛去集会上参加“你是最棒哒“的比赛。他的对手农夫保罗也带了M只牛去参加比赛
(1 ≤ N ≤ 1000, 1 ≤ M ≤ 1000)。每只牛都有自己的分数。两人会选择K只牛组成队伍(1 ≤ K ≤ 10)，两队
牛在按分数大小排序后一一配对，并且约翰打败保罗当且仅当对于每一对牛，约翰的牛分数都比保罗的高。请帮助
约翰计算约翰打败保罗的方案数 mod 1000000009。两种方案不同，当且仅当约翰或保罗选择的牛的集合与另一种
方案不同。

Input

The first line of input contains N, M, and K. The value of K will be no larger than N or M.
The next line contains the N scores of FJ’s cows.
The final line contains the M scores of FP’s cows.

Output

Print the number of ways FJ and FP can pick teams such that FJ wins, modulo 1,000,000,009.

Sample Input

10 10 3
1 2 2 6 6 7 8 9 14 17
1 3 8 10 10 16 16 18 19 19
Sample Output

382
HINT

Source

Platinum 鸣谢Firstlast提供译文

设dp[i][j]k]前i头fj的牛 j头fp的牛 配对k次的方案数

考虑dp[i+1][j][k]->dp[i+1][j+1][k]

dp[i][j+1][k]->dp[i+1][j+1][k]会算重复 所以减去即可

#include<bits/stdc++.h>
using namespace std;
inline char gc(){
    static char now[1<<16],*S,*T;
    if (T==S){T=(S=now)+fread(now,1,1<<16,stdin);if (T==S) return EOF;}
    return *S++;
}
inline int read(){
    int x=0,f=1;char ch=gc();
    while(!isdigit(ch)) {if (ch=='-') f=-1;ch=gc();}
    while(isdigit(ch)) x=x*10+ch-'0',ch=gc();
    return x*f;
}
const int N=1100;
const int mod=1e9+9;
inline void inc(int &x,int v){x=x+v>=mod?x+v-mod:x+v;}
inline void dec(int &x,int v){x=x-v<0?x-v+mod:x-v;}
int dp[N][N][15],n,m,K,a[N],b[N];
int main(){
    freopen("bzoj4742.in","r",stdin);
    n=read();m=read();K=read();
    for (int i=1;i<=n;++i) a[i]=read();
    for (int i=1;i<=m;++i) b[i]=read();
    sort(a+1,a+n+1);sort(b+1,b+m+1);dp[0][0][0]=1;
    for (int i=0;i<=n;++i){
        for (int j=0;j<=m;++j){
            for (int k=0;k<=K;++k){
                if (!dp[i][j][k]) continue;
                inc(dp[i+1][j][k],dp[i][j][k]);inc(dp[i][j+1][k],dp[i][j][k]);
                dec(dp[i+1][j+1][k],dp[i][j][k]);
                if (a[i+1]>b[j+1]) inc(dp[i+1][j+1][k+1],dp[i][j][k]);
            }
        }
    }printf("%d\n",dp[n][m][K]);
    return 0;
}