寒假训练——第四周（唯一分解定理又称算术基本定理）

最新推荐文章于 2024-07-19 22:36:19 发布

AC自动寄

最新推荐文章于 2024-07-19 22:36:19 发布

阅读量105

点赞数

分类专栏： vj集训补题文章标签：算法

本文链接：https://blog.csdn.net/m0_61409183/article/details/126592309

版权

vj集训补题专栏收录该内容

19 篇文章 0 订阅

订阅专栏

A - 樱花

思路：

约数个数 + 线性筛 + 阶乘求质因数的指数

具体实现：

变换公式：（思路差不多，盗图一下🤣）
如此我们仅需求 $n!)^{2}$ 的约数个数之和
先用线性筛求出所有质因数
在用阶乘求质因数的指数在 $O (l o g (n))$ 的时间复杂度求出答案

时间复杂度：

预计 $O(nlog(n))~但~O(n\sqrt{n})过不去$
实际：小于 $O (n l o g (n))$

代码如下：

#include <unordered_map>
#include <unordered_set>
#include <algorithm>
#include <iostream>
#include <cstring>
#include <string>
#include <bitset>
#include <vector>
#include <queue>
#include <cmath>
#include <map>
#include <set>

#define fast ios::sync_with_stdio(false), cin.tie(nullptr); cout.tie(nullptr)

#define mkpr make_pair
#define x first
#define y second
#define int long long

using namespace std;

typedef long long LL;
typedef pair<int, int> PII;

const int mod = 1e9 + 7;
const int INF = 0x3f3f3f3f;
const LL LL_INF = 0x3f3f3f3f3f3f3f3f;
const double eps = 1e-9;

const int N = 1e6 + 10, M = N * 2;
const int YB = 8, YM = 1e8;

const int dx[4] = {-1, 0, 1, 0}, dy[4] = {0, 1, 0, -1};

int T, cases;
int n, m, times;

int primes[N], cnt;
bool st[N];

void init(int n)
{
    for (int i = 2; i <= n; i ++ )
    {
        if(!st[i]) primes[cnt ++ ] = i;
        for (int j = 0; primes[j] <= n / i; j ++ )
        {
            st[i * primes[j]] = true;
            if(i % primes[j] == 0) break;
        }
    }
}

void solve()
{
    cin >> n;

    init(n);
    
    int res = 1;
    for (int i = 0; i < cnt; i ++ )
    {
        int p = primes[i];
        int s = 0;
        
        // 求 n! 的 质因数p 的指数的简单方法
        for (int j = n; j; j /= p) s += j / p; 
        
        res = res * (s* 2 + 1) % mod;
    }
    
    cout << res << endl;

    return;
}

signed main()
{
    T = 1;
    //fast;cin >> T;
    //scanf("%d", &T);

    //for (cases = 1; cases <= T; cases ++ )
    while(T -- )
        solve();

    return 0;
}