洛谷P1591题解

一、问题引出

阶乘数码

题目描述

求 $n!$ 中某个数码出现的次数。

输入格式

第一行为 $t(t\le 10)$，表示数据组数。接下来 $t$ 行，每行一个正整数 $n(n\le 1000)$ 和数码 $a$。

输出格式

对于每组数据，输出一个整数，表示 $n!$ 中 $a$ 出现的次数。

样例 #1

样例输入 #1

2
5 2
7 0

样例输出 #1

1
2

二、求解思路

这是一道很经典的高精度*单精度的问题，我们应记住这种题目的模板，在之前的高精度乘以高精度的题目当中，我们使用了字符串去保存每次计算的中间结果，但是在这里，我们使用数组去存储每次计算所得的结果，因为之前的使用字符串的方法会超时，虽然结果是正确的，在这里我们也附上采用字符串方式的结果：

#include <iostream>
using namespace std;
string add(string s,string t)
{
    int len1=s.length();
    int len2=t.length();
    if(len1<len2)
    {
        for(int i=0;i<(len2-len1);i++)
        {
            s="0"+s;
        }
    }
    else
    {
        for(int i=0;i<(len1-len2);i++)
        {
            t="0"+t;
        }
    }
    int len=s.length();
    int r=0;
    string result;
    for(int i=len-1;i>=0;i--)
    {
        int sum=(int(s[i]-'0')+int(t[i]-'0')+r)%10;
        r=(int(s[i]-'0')+int(t[i]-'0')+r)/10;
        result=char(sum+'0')+result;
    }
    if(r!=0)
        result=char(r+'0')+result;
    return result;
}
string mul(string s,string t)
{
    if (s=="0" || t=="0")
    {
        return "0";
    }  
    int len1=s.length();
    int len2=t.length();
    string result;
    for(int i=len2-1;i>=0;i--)
    {
        string midResult;
        //获取乘数的值方便待会进行运算
        int x=t[i]-'0';
        for(int j=0;j<x;j++)
        {
            midResult=add(midResult,s);
        }
        //移位操作
        for (int j = 0; j < len2-1-i; j++)
        {
            midResult+="0";
        }
        result=add(result,midResult);
    }
    return result;
}
string Factorial(int n)
{
    string result="1";
    for (int i = 1; i <= n; i++)
    {
        result=mul(result,to_string(i));
    }
    return result;
}
int main()
{
    int t,n,a,sum;
    cin>>t;
    for (int i = 0; i < t; i++)
    {
        sum=0;
        cin>>n>>a;
        string s=Factorial(n);
        // cout<<s;
        for (int j = 0; j < s.length(); j++)
        {
            if (s[j]-'0'==a)
            {
                sum++;
            }  
        }
        cout<<sum<<endl;
    }
}

三、正确思路

#include <iostream>
using namespace std;
/**
 * @brief 
 * 数组a用于保存阶乘的结果，因为n的值最多也就是1000，
 * 最终结果的位数不会超过50000位，所以我们选择将数组的大小设置为50000
 */
int a[5000];
int main()
{
    int t;
    cin>>t;
    while (t--)
    {
        int n,m,sum=0;
        cin>>n>>m;
        /**
         * @brief 
         * 1.记得每轮把数组a清空
         * 2.数组a的第1位要设为1，不然没法去乘
         */
        for (int i = 1; i <= 11; i++)
        {
            a[i]=0;
        }
        a[1]=1;
        /**
         * @brief 
         * 我们求阶乘的大体思路就是：
         * 1.对每次的数组a，我们都用2~n之间的数去乘以它，并保存乘得的中间结果，
         *   这样n-1轮过后a中存放的数就是最终的结果（前面放低位）
         * 2.每次对a乘以i的时候，需要考虑进位等事项
         * 3.用jw去表示进位
         * 4.用p表示每次数组a保存的中间结果的长度，初始为1
         * 5.用变量j去记录
         */
        int jw=0,p=1,j;
        for (int i = 2; i <= n; i++)
        {
            for (j = 1; j <=p; j++)
            {
                a[j]=a[j]*i+jw;
                jw=a[j]/10;
                a[j]%=10;
            }
            /**
             * @brief 
             * 如果还有进位，就往前放
             * 记得j每次都得++，因为j记录的是中间结果的长度
             * 记得更新p
             */
            while (jw>0)
            {
                a[j]=jw%10;
                jw/=10;
                j++;
            }
            p=j-1;
        }
        for (int i = 1; i <= p; i++)
        {
            if (a[i]==m)
            {
                sum++;
            }
        }
        cout<<sum<<endl;
    }
}