P1591 阶乘数码 （高精阶乘）

题目描述

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

输入格式

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

输出格式

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

输入输出样例

输入 #1

2
5 2
7 0

输出 #1

1
2

原题：https://www.luogu.com.cn/problem/P1591

 

高精阶乘，算完后统计下出现了多少个a即可。先放AC代码：

#include <iostream>
#include <cmath>
#include <string>
#include <cstring>
#include <algorithm>

using namespace std;


typedef long long LL;

//手动memset整个数组为-1
void bitset(char *s, int n)
{
    for (int i = 0; i < n; ++i)
        s[i] = -1;
}

//将-1作为数组结束的标志来复制数组
void arrcpy(char *s, char *s2)
{
    for (int i = 0; s2[i] != -1; ++i)
        s[i] = s2[i];
}

//以-1为数组结束的标志来计算数组内的有效长度
int len(char *s)
{
    int res = 0;
    while (s[res++] != -1)
        ;
    return res - 1;
}

//输出数组，用于debug
void show(char *s)
{
    for (int i = 0; s[i] != -1; ++i)
        cout << char(s[i] + '0');
    cout << endl;
}

LL fa(int n, int a)
{

    char f1[10000];
    char f2[4];
    char res[10000];

    bitset(f1, 10000);
    bitset(f2, 4);
    bitset(res, 10000);

    f1[0] = 1;//初始为1
    int l1;//l1 为 f1的有效长度，l2为f2的有效长度
    for (int i = 1; i <= n; ++i)
    {
        bitset(f2, 4);
        int tmp = i;
        int l2 = 0;
        while (tmp)//构造每次的乘数
        {
            f2[l2++] = tmp % 10;
            tmp /= 10;
        }
        l1 = len(f1);//由于每轮都会将res拷贝到f1里，所以需要统计一次长度
        int m;
        for (int j = 0; j < l1; ++j)
        {
            m = 0;
            int k;
            for (k = 0; k < l2; ++k)
            {
                m += f1[j] * f2[k];
                if (res[j + k] > 0)
                    m += res[j + k];
                res[j + k] = m % 10;
                m /= 10;
            }
            if (m)
                res[j + k] = (res[j + k] > 0 ? res[j + k] : 0) + m;
        }
        bitset(f1, 10000);
        arrcpy(f1, res);
        bitset(res, 10000);
        //show(f1);
    }
    
    //统计a出现的次数
    int c = 0;
    for (int i = 0; f1[i] >= 0; ++i)
        if (f1[i] == a)
            ++c;

    return c;
}

int main()
{
    //freopen("1.txt", "r", stdin);

    int t;
    cin >> t;
    int n, a;
    while (t--)
    {
        cin >> n >> a;
        cout << fa(n, a) << endl;
    }

    return 0;
}

我这里使用三个char数组来存放两个乘数和积，每计算完一次就把res数组拷贝到f1数组里继续下一轮计算，要注意的是高精乘法必须得有一个数组来放计算结果，一共需要三个数组，加法则可以不用，只需要两个数组，另外，由于题目没要求输出之类的，所以计算过程中所有数字都是倒序的，也不用翻转了。