浙江机电职业技术学院第八届新生亮相赛（同步赛）D.一种很新的阶乘【线性筛+一分为二法】

原题链接：https://ac.nowcoder.com/acm/contest/73799/D

时间限制：C/C++ 2秒，其他语言4秒
空间限制：C/C++ 262144K，其他语言524288K
64bit IO Format: %lld

题目描述

夜雷很喜欢特别大的数字，曾经一度算阶乘算了好几页纸

现在夜雷定义了一种很新的阶乘  f(x) = x * (x-1)^2 * (x-2)^3 *………* 2^(x-1) * 1^x

夜雷算麻了，小小的脑子装不下那么多数，他希望你能帮他算一算。

为了省下一点纸，他希望你能帮他算出该函数的质因数分解形式

输入描述:

一行，一个整数 x（1<x<=1e7），表示运算式的输入。

输出描述:

一个字符串，表示f(n)的质因子分解形式，要求按照质因子从小到大排列

示例1

输入

5

输出

f(5)=2^8*3^3*5

解题思路：

首先根据题目表达式只需要对1~n中的每一个数进行质因数分解即可，对于指数只要对分解出来的所有质因数乘上相应的值即可，1不需要分解，所以只需要考虑2~n即可，分解质因数时间为O(根号n)，对于1~n中的每个数都分解质因数，时间为O(根号1+根号2+根号3+...+根号n)，这个时间不是很好计算，我写个代码计算一下大概=21081852648，大概是2e10，这个时间复杂度非常高了，肯定过不了，我们考虑看有没有什么优化方式。

首先我们知道对于任意一个非质数k肯定都可以分解为俩个数a和b的乘积，我们让a为质数，那么b可能为质数，可能不是质数，但是a和b都比k小，例如k=a*b，cnt[a]+=k,cnt[b]+=k,cnt[k]=0，我们从大到小枚举n~2,那么就能将所有的大的非质数分解为小的质数，这样每个数只需要枚举一次，时间复杂度为O(n),n=1e7，这个时间是肯定可以过的，对于质数就不需要分解了，根据上面分析我们还需要先拿到1~n中的所有质数，我们线性筛拿到所有质数即可。

总结：这个题目出的挺好的，不是很难而且运用到了一些很好的数学优化技巧。

时间复杂度：首先线性筛时间复杂度为O(n)，然后枚举每一个数时间复杂度为O(n)，最终时间复杂度为O(n)。

空间复杂度：空间复杂度为O(n)。

cpp代码如下：

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

using namespace std;
typedef long long LL;
const int N=1e7+10;

LL cnt[N],vis[N];
int primes[N], number;
bool st[N];

void get_primes(int n)
{
    for (int i = 2; i <= n; i ++ )
    {
        if (!st[i]) primes[number ++ ] = i,vis[i]=i;
        for (int j = 0; primes[j] <= n / i; j ++ )
        {
            vis[primes[j]*i]=primes[j];
            st[primes[j] * i] = true;
            if (i % primes[j] == 0) break;
        }
    }
}
int main()
{
    int x;
    cin>>x;
    get_primes(x);  //线性筛
    //枚举每一个数进行一分为二
    for(int i=1;i<x;i++)
    {
        int k=x-i+1;
        cnt[k]+=i;
        if(vis[k]==k)continue;  //vis[k]==k表示k为质数
        int x1=vis[k],x2=k/vis[k];
        cnt[x1]+=cnt[k],cnt[x2]+=cnt[k];
        cnt[k]=0;
    }
    //按照要求的输出格式输出答案
    printf("f(%d)=",x);
    bool ok=false;
    for(int i=2;i<=x;i++)
        if(cnt[i]!=0)
        {
            if(ok)printf("*");
            printf("%d",i);
            if(cnt[i]==1){
                continue;
            }
            printf("^%lld",cnt[i]);
            ok=true;
        }
    return 0;
}