C 解法, 执行用时: 2ms, 内存消耗: 364KB, 提交时间: 2019-08-04

  
int main()
{
    long long  mod =1000000007;
  
    int cnk[51][51];
    int s;
    scanf("%d",&s);
    int x[3];
    int i,j,t;
    for(i=0;i<3;i++)
    {
        scanf("%d",&x[i]);
          
    }
  
    for (i = 1, cnk[0][0] = 1; i<51; i++) //动态规划求c(n,k)
    {
        cnk[i][0] = 1;
        for (j = 1; j <51; j++)
        {
            cnk[i][j] = (cnk[i - 1][j - 1] + cnk[i - 1][j]) % mod; //没有取余的话会溢出
        }
    }
  
    for(i=0;i<3-1;i++)//对x排序，从小到大,保证x2最大，x0最小
    {
        for(j=i+1;j<3;j++)
        {
            if(x[i]>x[j])
            {
                t=x[i];
                x[i]=x[j];
                x[j]=t;
            }
        }
    }
      
  
    long long cnt=0;
    int left;
    int tmp;
    left = s - x[2]; //给工作量最大的分配后，剩余工作量
  
  
    for(i = 0; i <= left; i++) //i是x1在left中分担的工作量
    {
        if(i <= x[1] && left - i <= x[0])//left - i是x0在left中分担的工作量
        {
            tmp = cnk[x[2]][x[1] - i] % mod; // 在x2中，有x[1] - i分配给x1
            tmp = (tmp * cnk[left][i]) % mod; // left中有i个分配给x1
            tmp = (tmp * cnk[s - (left - i)][x[0] - (left - i)]) % mod; // x0的随意挑选部分
            cnt = (cnt + tmp) % mod; //每个case求和
        }
           
    }
    cnt = (cnt * cnk[s][x[2]]) % mod; //求和后与最初的情况求积
      
    printf("%d\n",cnt);
       
    return 0;
}

C++ 解法, 执行用时: 2ms, 内存消耗: 376KB, 提交时间: 2020-08-11

#include <iostream>
#include <vector>
#include <math.h>
#include <algorithm>
using namespace std;

int main() {
    int s;
    vector<int> member(3);
    while (cin >> s >> member[0] >> member[1] >> member[2]) {
        sort(member.begin(), member.end());
        vector<vector<int> > C(s  + 1, vector<int> (s + 1));
        C[0][0] = 1;
        for (int i = 1; i <= s; i++) {
            C[i][0] = 1;
            for (int j = 1; j <= s; j++) {
                C[i][j] = (C[i - 1][j - 1] + C[i - 1][j]) % 1000000007;
            }
        }
        int left = s - member[2];
        long temp, ans = 0;
        for (int i = 0; i <= left; i++) {
            if (i <= member[1] && left - i <= member[0]) {
                temp = C[member[2]][member[1] - i] % 1000000007;
                temp = temp * C[left][i] % 1000000007;
                temp = temp * C[s - left + i][member[0] - left + i] % 1000000007;
                ans = (ans + temp) % 1000000007;
            }
        }
        ans = (ans * C[s][member[2]]) % 1000000007;
        cout << ans << endl;
    }
    return 0;
}