2023百度之星第三场新的阶乘题解（数论/分解质因数）

不要麻酱

已于 2023-09-28 20:59:57 修改

阅读量451

点赞数

分类专栏：算法竞赛文章标签： c++ 算法

于 2023-09-27 22:25:10 首次发布

本文链接：https://blog.csdn.net/weixin_62701031/article/details/133364721

版权

算法竞赛专栏收录该内容

17 篇文章 0 订阅

订阅专栏

题目

思路

用线性筛把每个合数的分解成最小质数prims[j] 乘以一个常数x（x可能是质数可能是合数）的形式，然后从后往前操作，把每个合数的x再按照最小质数prims[j] 乘以一个常数x的方法分解，这样就能得到质数的次数。

代码

#include<bits/stdc++.h>
using namespace std;
#define int long long
const int mod = 1e9 + 7, N = 1e7 + 10;
int n, st[N], is_primes[N], ans[N];
pair<int, int> pre[N];
vector<int>primes;
void get()
{
    for(int i = 2; i <= n; i ++)
    {
        if(!st[i]) primes.push_back(i), is_primes[i] = 1, pre[i] = {i, 1};
        for(int j = 0; primes[j] * i <= n; j ++)
        {
            int t = primes[j] * i;
            pre[t] = {primes[j], i};
            st[t] = 1;
            if(i % primes[j] == 0) break;
        }
    }
    for(int i = n; i >= 1; i --)
    {
        if(is_primes[i]) continue;
        auto [x, y] = pre[i];
        ans[x] += ans[i], ans[y] += ans[i];
    }
}
signed main()
{
	ios::sync_with_stdio(0);
	cin.tie(0); cout.tie(0);
    cin >> n;
    for(int i = 2; i <= n; i ++) ans[i] = n + 1 - i;
    get();
    cout << "f(" << n << ")=";
    int cnt = 0; 
    for(int i = 2; i <= n; i ++)
    {
        if(is_primes[i])
        {
            if(cnt) cout << "*";
            if(ans[i] == 1) cout << i;
            else cout << i << "^" << ans[i];
            cnt ++;
        }

    }
}