Codeforces Round 1292 简要题解

A. NEKO’s Maze Game

略

B. Aroma’s Search

略

C. Xenon’s Attack on the Gangs

略

D. Chaotic V.

vp的时候脑子抽了，不知道写了啥，最后也没过。
令$k=\max (k_i)$，考虑将$1!\sim k!$分解质因数，这里可以从小到大计算，每次在上次基础上分解最后一个数即可。将一个数$x$表示为每个质因子出现次数的集合$(c_1,c_2,...,c_m)$，那么$x$的父亲即为将第一个不为$0$的$c_i$减去$1$后所得的集合，也可以方便计算两个集合间的距离。
注意到一定存在一个最优的$P$是将某个$x!$对应集合某个前缀变为$0$所得到的，证明可以考虑对于某个$P$，找到最接近（公共后缀最长）的$x!$对应集合，那么显然需要满足最后一个不同的位$p$之前都是$0$，并且若位$p$不为$0$，那么显然可以选择变成一个与$P$从$p+1$位往后都相同的$y!$的对应集合的$p$往后的后缀，且答案不会变劣。
这样就可以暴力枚举每个$x!$的后缀作为可能的$P$计算答案了，实现的时候考虑预处理出了阶乘之间的最长公共后缀（只需要枚举相邻两个暴力计算，然后取中间的$\min$即可）以及一些前缀和，即可轻松从后往前计算答案。
时间复杂度为$\mathcal{O}(k^2+n)$。

#include <bits/stdc++.h>
#define FR first
#define SE second
#define inf 0x3f3f3f3f3f3f3f3fLL

using namespace std;

typedef long long ll;
typedef pair<int,int> pr;

const int M=5000;

int prime[5005],tot;
bool check[5005];

int num[5005][1005],val[5005],trans[5005];
ll sum1[1005][5005],sum2[1005][5005],sum4[1005][5005],sum5[1005][5005],sum3[5005];

void pre() {
   
  for(int i=2;i<=M;i++) {
   
  	if (!check[i]) prime[++tot]=i;
  	for(int j=1;j<=tot&&i*prime[j]<=M;j++) {
   
  		check[i*prime[j]]=1;
  		if (i%prime[j]==0) break;
	  }
  }
  for(int i=1;i<=M;i++) {
   
  	memcpy(num[i],num[i-1],sizeof(num[i]));
  	int x=i;
  	trans[i]=1;
  	for(int j=1;j<=tot&&prime[j]<=i;j++) {
   
  		while (x%prime[j]==0) {
   
  			num[i][j]++;
  			x/=prime[j];
  			trans[i]=j+1;
		  }
	  }
  }
  for(int i=1;i<=tot;i++)
    for(int j=1;j<=M;j++) sum1[i][j]=sum1[i][j-1]+num[j][i];
  for(int i=1;i<=tot;i++)
    for(int j=1;j<=M;j++) sum2[i][j]=sum2[i-1][j]+sum1[i][j];
  for(int i=1;i<=tot;i++)
    for(int j=1;j<=M;j++) sum4[i][j]=sum4[i][j-1]+(ll)num[j][i]*val[j];
  for(int i=1;i<=tot;i++)
    for(int j=1;j<=M;j++) sum5[i][j]=sum5[i-1][j]+sum4[i][j];
  for(int i=1;i<=M;i++) sum3[i]=sum3[i-1]+val[i];
}

vector <int> vt[5005];

ll solve(int x,int n) {
   
  for(int i=1;i<=tot+2;i++) vt[i].clear();
  int maxn=1;
  for(int i=x-1;i>0;i--) {
   
  	maxn=max(maxn,trans[i+1]);
  	vt[maxn].push_back(i);
  }
  maxn=1;
  for(int i=x+1;i<=M;i++) {
   
  	maxn=max(maxn,trans[i]);
  	vt[maxn].push_back(i);
  }
  ll ans=inf,s1=0,s2=0;
  int lx=1,rx=M;
  for(int i=tot;i>0;i--) {
   
  	s2+=s1*num[x][i+1];
  	for(int j=0;j<vt[i+2].size();j++) {
   
  		int u=vt[i+2][j];
  		if (u<x) lx=max(lx,u+1); else rx=min(rx,u-1);
  		s1+=val[u];
		s2+=(sum2[i][u]-sum2[i][u-1])*val[u]+(ll)abs(num[u][i+1]-num[x][i+1])*val[u]; 
	  }
	ll s=s2+sum5[i-1][rx]-sum5[i-1][lx-1];
	ans=min(ans,s+sum4[i][rx]-sum4[i][lx-1]);
	if (i==1) {
   
		int v=num[x][i];
		s+=s1*v;
		s+=(sum3[x]-sum3[lx-1])*v-(sum4[i][x]-sum4[i][lx-1]);
		s+=(sum4[i][rx]-sum4[i][x])-(sum3[rx]-sum3[x])*v;
		ans=min(ans,s);
	}
  }
  return ans;
}

int main() {
   
  int n;
  scanf("%d",&n);
  for(int i=1;i<=n;i++) {
   
  	int x;
  	scanf("%d",&x);
  	if (!x) x=1;
  	val[x]++;
  }
  pre();
  ll ans=inf;
  for(int i=1;i<=M;i++)
    ans=min(ans,solve(i,n));
  printf("%lld\n",ans