[BZOJ 2790]POI2012 Distance

乍一看以为是神题，想了想，原来很暴力。。。

可以看到乘和除是没有顺序关系的，那么为了我们就统一先除再乘。

假设我们已经把原数除完了变成x，那么再在x上乘一些数使得得到的数在所有书中所需的步数最小，记为min1，次小的几为min2。

预处理大法~~~~

对每个数字枚举所有约数d，更新的最小值。

在计算答案，同样枚举约数，利用min1或min2计算即可。

代码

#include <cstdio>
#include <cstring>
#include <algorithm>

using namespace std;
const int INF=1e9;
int check[1000005],prime[1000000],cnt,ans,i,j,tot,x,n,minx;
int p[100005][10],k[100005][10],s[100005],t[100005];
int min1[1000005],min2[1000005],num1[1000005],num2[1000005];

void init(){
for (i=2;i<=1000000;i++){
if (check[i]==0) check[i]=i, prime[++tot]=i;
for (j=1;j<=tot;j++){
if (i*prime[j]>1000000) break;
check[i*prime[j]]=prime[j];
if (i%prime[j]==0) break;
}
}
}

void dfs(int x,int dig,int cnt){
if (x>t[i]){
if (min1[dig]>cnt){
min2[dig]=min1[dig]; num2[dig]=num1[dig];
min1[dig]=cnt; num1[dig]=i;
} else
if (min2[dig]>cnt){
min2[dig]=cnt;
num2[dig]=i;
}
return;
}
for (int j=0;j<=k[i][x];j++){
dfs(x+1,dig,cnt-j);
dig=dig*p[i][x];
}

}

void work(int x,int dig,int cnt){
if (x>t[i]){
if (i==num1[dig]){
if (min2[dig]+cnt<minx || (min2[dig]+cnt==minx &&  ans>num2[dig])){
ans=num2[dig];
minx=min2[dig]+cnt;
}
} else
{
if (min1[dig]+cnt<minx || (min1[dig]+cnt==minx &&  ans>num1[dig])){
ans=num1[dig];
minx=min1[dig]+cnt;
}
}
return;
}
for (int j=0;j<=k[i][x];j++){
work(x+1,dig,cnt-j);
dig=dig*p[i][x];
}

}

int main(){
//freopen("odl.in","r",stdin);
//freopen("odl.out","w",stdout);
init();
scanf("%d",&n);
memset(min1,127/2,sizeof(min1));
memset(min2,127/2,sizeof(min2));
for (i=1;i<=n;i++){
scanf("%d",&x);
s[i]=0; t[i]=0;
while (x>1){
p[i][ ++t[i] ] = check[x];
k[i][ t[i] ] = 0;
while (x%p[i][t[i]]==0)
x/=p[i][t[i]], k[i][ t[i] ]++, s[i]++;
}
dfs(1,1,s[i]);
}
for (i=1;i<=n;i++){
minx=INF; ans=INF;
work(1,1,s[i]);
//printf("%d %d\n",minx,ans);
printf("%d\n",ans);
}
return 0;
}