uoj22 外星人（dp）

题目大意：
给定一个$n$个数的序列$a$，给定一个$x$，其中$a$数组可以进行顺序的调换，每一个$a_i$都能使$x=x \mod a_i$,
求最后经过一系列计算后的$y$，满足$abs(x-y)$尽可能小，并求出方案数

QwQ 哇，一看到这个题。说实话，没什么好的思路。

也就发现了几个性质：

1.最后的$y$一定小于最小的$a_i$

2.如果存在一个$a_i<a_j$，且$i<j$ 那么$a_j$就没有任何作用了，对答案没有任何一点影响

那我们不妨将整个数组从大到小排序

先考虑第一问:
我们定义$f[i][j]$表示，考虑到第$i$个数，当前的值为$j$是否可行，首先我们令$f[0][x]=1$，然后对于当前的$i$，我们可以选择用它 ，也可以选择不用（换句话说，就是放一个比它更小的在前面，就可以实现不使用它了）但是后者需要满足$i\ !=n$ 然后分别对应转移即可

这里有部分分的代码！

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cmath>

using namespace std;

const int maxn = 1010;
const int maxx = 5010;

int f[maxn][maxx];
int a[maxn];
int n,x;

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

bool cmp(int a,int b)
{
    return a>b;
}

int main()
{
    n=read(),x=read();
    for (int i=1;i<=n;i++) a[i]=read();
    sort(a+1,a+1+n,cmp);
    f[0][x]=1;
    for (int i=1;i<=n;i++)
    {
        for (int j=0;j<=x;j++) f[i][j%a[i]]=max(f[i][j%a[i]],f[i-1][j]);
        if (i!=n) for (int j=0;j<=x;j++) f[i][j]=max(f[i][j],f[i-1][j]);
    }
    for (int i=x;i>=0;i--) if (f[n][i]) {
        cout<<i<<endl<<0<<endl;
        return 0;
    }
    return 0;   
}

那么加上第二问呢，该怎么解决呢。

看了一些排列组合的题解，不过并不知道怎么做呀。倒是有一种更好理解的方法QwQ

我们令$g[i][j]$表示处理第$i$个数，当前值是$j$的方案数

如果我们使用这个点$g[i][j \mod a_i ]+=g[i-1][j]$（说明他待在当前的位置，且后面比他小的位置，都在他后面

如果不用$g[i][j]=g[i-1][j]*(n-i)$ （表示他可以和他之后的任意一个比他小的数换位置，都不会使用这个点）（或者理解为他有$n-i$个空隙可以插进去

直接上代码

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
#include<cmath>

using namespace std;

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

const int maxn = 1010;
const long long mod =  998244353;

int f[maxn][5010];
long long g[maxn][5010];
int n,x;
int a[maxn];

bool cmp (int a,int b)
{
    return a>b;
}

int main()
{
  scanf("%d%d",&n,&x);
  for (int i=1;i<=n;i++) a[i]=read();
  f[0][x]=1;
  g[0][x]=1;
  sort(a+1,a+1+n,cmp);
  for (int i=1;i<=n;i++)
  {
    for (int j=0;j<=x;j++) f[i][j%a[i]]=max(f[i][j%a[i]],f[i-1][j]),g[i][j%a[i]]=(g[i][j%a[i]]+g[i-1][j])%mod;
    if (i!=n) for (int j=0;j<=x;j++) f[i][j]=max(f[i][j],f[i-1][j]),g[i][j]=(g[i][j]+g[i-1][j]%mod*(long long
      )(n-i)%mod)%mod; 
  }  
  for(int i=a[n];i>=0;i--)
  {
     if (f[n][i])
     {
        cout<<i<<endl;
        cout<<g[n][i]<<endl;
        return 0;
     }
  }
  return 0;
}