【UOJ22】. [UR #1]外星人

题目链接

题解

http://vfleaking.blog.uoj.ac/blog/33
先从大到小排序，$f[i][j]$表示当前处理到$a_i$，现在的结果是j可不可能。再用$g[i][j]$记录有几种得到当前情况的可能。
那么对于每个$a_i$，要么就把它放在当前位置，现在生效，要么把它放在后面的$n-i$个位置，使它永不生效，因为如果你先模了一个小于$a_i$的数，再模$a_i$结果是不会变的。那么转移就很简单了。

#include<bits/stdc++.h>
using namespace std;

inline int read(){
    int x = 0, f = 1; char c = getchar();
    while(!isdigit(c)) { if(c == '-') f = -1; c = getchar(); }
    while(isdigit(c)) { x = x * 10 + c - '0'; c = getchar(); }
    return x * f;
}

typedef long long ll;
const int N = 1000 + 10, X = 5000 + 10, mo = 998244353;
int n, x, a[N];
bool f[N][X];
ll g[N][X];

inline bool cmp(int x, int y){
    return x > y;
}

void work(){
    f[0][x] = 1; g[0][x] = 1;
    for(int i = 1; i <= n; i++){
        for(int j = 0; j <= x; j++){
            f[i][j%a[i]] |= f[i-1][j];
            (g[i][j%a[i]] += g[i-1][j]) %= mo;
        }
        if(i < n){
            for(int j = 0; j <= x; j++){
                f[i][j] |= f[i-1][j];
                (g[i][j] += g[i-1][j] * (n - i)) %= mo;
            }
        }
    }
    for(int i = x; i >= 0; i--)
        if(f[n][i]){
            printf("%d\n%lld\n", i, g[n][i]);
            break;
        }
}

void init(){
    n = read(), x = read();
    for(int i = 1; i <= n; i++) a[i] = read();
    sort(a + 1, a + n + 1, cmp);
}

int main(){
    init();
    work();
    return 0;
}