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;
}

原文地址:https://www.cnblogs.com/yimmortal/p/10160839.html