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

版权声明:本文为博主原创文章,注明作者和出处即可转载。 https://blog.csdn.net/KikiDMW/article/details/70142714

题目链接
这里写图片描述

题解

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

#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;
}
阅读更多

没有更多推荐了,返回首页