Description
给出一个 n×n n × n 矩阵,问从左上角到右下角的所有路径中(每次只能往下或往右),经过的数字之积末尾连续的 0 0 的数量最少是多少
Input
第一行一整数表示矩阵行列数,之后输入一个 n×n n × n 矩阵 Aij A i j (2≤n≤1000,0≤Aij≤109) ( 2 ≤ n ≤ 1000 , 0 ≤ A i j ≤ 10 9 )
Output
输出路径经过数字之积末尾连续的 0 0 的最少数量和这条路径
Sample Input
3
1 2 3
4 5 6
7 8 9
Sample Output
0
DDRR
Solution
经过数字之积末尾连续的数量即为积中 2 2 和作为因子出现次数的较小值,由于最优方案对应的路径,其 2 2 的数量或的数量必然是所有路径中最少的,故分别求出 2 2 最少的路径和最少的路径,取一个较小值即可,注意由于矩阵中可能有 0 0 存在,故如果上述方法求出的答案大于,那么不如选择走有 0 0 的路径,这样答案是
求 2 2 作为因子出现最少的路径,首先求出每个元素中因子的数量 a[i][j] a [ i ] [ j ] ,设 dp[i][j] d p [ i ] [ j ] 为从 (1,1) ( 1 , 1 ) 到 (i,j) ( i , j ) 的路径中 2 2 的最少数量,那么有转移,在转移过程中记录路径即可,求 5 5 <script type="math/tex" id="MathJax-Element-27">5</script>作为因子出现最少的路径同理可得
Code
#include<cstdio>
#include<iostream>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<vector>
#include<queue>
#include<map>
#include<set>
#include<ctime>
using namespace std;
typedef long long ll;
typedef pair<int,int>P;
const int INF=0x3f3f3f3f,maxn=1005;
int n,dp[maxn][maxn],a[maxn][maxn],b[maxn][maxn];
char pa[maxn][maxn],pb[maxn][maxn],path[maxn];
int deal(int x,int y)
{
if(x==0)return 0;
int ans=0;
while(x%y==0)x/=y,ans++;
return ans;
}
int main()
{
scanf("%d",&n);
int flag=0,x0,y0;
for(int i=1;i<=n;i++)
for(int j=1;j<=n;j++)
{
int temp;
scanf("%d",&temp);
if(temp==0)flag=1,x0=i,y0=j;
a[i][j]=deal(temp,2),b[i][j]=deal(temp,5);
}
for(int i=1;i<=n;i++)dp[1][i]=dp[1][i-1]+a[1][i],dp[i][1]=dp[i-1][1]+a[i][1],pa[1][i]='R',pa[i][1]='D';
for(int i=2;i<=n;i++)
for(int j=2;j<=n;j++)
{
if(dp[i-1][j]>dp[i][j-1])pa[i][j]='R';
else pa[i][j]='D';
dp[i][j]=min(dp[i-1][j],dp[i][j-1])+a[i][j];
}
int ans=dp[n][n];
for(int i=1;i<=n;i++)dp[1][i]=dp[1][i-1]+b[1][i],dp[i][1]=dp[i-1][1]+b[i][1],pb[1][i]='R',pb[i][1]='D';
for(int i=2;i<=n;i++)
for(int j=2;j<=n;j++)
{
if(dp[i-1][j]>dp[i][j-1])pb[i][j]='R';
else pb[i][j]='D';
dp[i][j]=min(dp[i-1][j],dp[i][j-1])+b[i][j];
}
if(flag&&min(ans,dp[n][n])>1)
{
printf("1\n");
for(int i=1;i<x0;i++)printf("D");
for(int i=1;i<n;i++)printf("R");
for(int i=x0;i<n;i++)printf("D");
printf("\n");
}
else if(ans<=dp[n][n])
{
printf("%d\n",ans);
int x=n,y=n;
for(int i=0;i<2*n-2;i++)
{
path[2*n-3-i]=pa[x][y];
if(pa[x][y]=='R')y--;
else x--;
}
path[2*n-2]='\0';
printf("%s\n",path);
}
else
{
printf("%d\n",dp[n][n]);
int x=n,y=n;
for(int i=0;i<2*n-2;i++)
{
path[2*n-3-i]=pb[x][y];
if(pb[x][y]=='R')y--;
else x--;
}
path[2*n-2]='\0';
printf("%s\n",path);
}
return 0;
}