码迷,mamicode.com
首页 > 其他好文 > 详细

洛谷 [P1436] 反素数

时间:2018-03-07 20:17:41      阅读:188      评论:0      收藏:0      [点我收藏+]

标签:素数   最小   oid   include   cin   sim   约数个数   prim   cout   

算术基本定理的应用

算术基本定理:

一个正整数 \(N\) 能唯一分解成如下形式
\[N=p_1 ^ {c_1}p_2^{c_2}\cdots P_m ^ {c_m}\]
其中 \(c_i\) 都是正整数, \(P_i\)都是质数

推论:

\(N\) 的正约数集合可以写成如下形式:
\[ \{p_1^{b_1}P_2^{b_2}\cdots P_m^{b_m}\}\], 其中\(0\leq b_i \leq c_i\)

\(N\) 的正约数个数可表示为:
\[ (c_1 + 1) * (c_2 + 1) * \cdots *(c_m + 1) = \prod_{i=1}^m\quad(c_i + 1)\]

\(N\) 的所有正约数的和可表示为:
\[(1 + p_1 + p_1 ^2 + \cdots p_1 ^ {c_1})* \cdots *(1 + p_m + p_m ^2 + \cdots + p_m ^{c_m}) = \prod_{i =1}^1(\sum_{j = 1}^{c_1}(p_j)^i)\]

首先,对于本题来说,反素数一定是$ 1 \sim N $中约数个数最多的数中最小的一个

其次\(1 \sim N\)的的所有数中,拥有的不相同的质因子不会超过 10 个
每个质因子的幂次不会超过 30

最后我们发现,答案 x 的质因子是连续的若干个质数,且质数单调递减,
如果不是这样的话,我们一定可以通过交换质因子次数的方法,得到一个约数一样,但更小的答案

注意:本题要开 long long

#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <cmath>
#define ll long long
using namespace std;
ll n, ans, cnt, prime[15] = {0, 2, 3, 5, 7, 11, 13, 17, 19, 23, 29, 31};
void dfs(int k, ll num, ll c, int cur) {
    //cout<<k<<endl;
    if(k == 11 || !cur) {
        if(cnt < c) {ans = num;cnt = c;}
        if(cnt == c) ans = min(ans, num);
        return;
    }
    ll ttt = 1;
    for(int i = 0 ; i <= cur ; i++) {
        if(i) ttt *= prime[k];
        if(num * ttt > n) return;
        dfs(k + 1, num * ttt, c * (i + 1), i);
    }
}
int main() {
    cin>>n;
    if(n == 1) {printf("%lld\n", n);return 0;}
    int k = 0;
    for( ; (1ll<<k) < n && k <= 30; k++) ;
    dfs(1, 1, 1, k);
    cout<<ans<<endl;
    return 0;
}

洛谷 [P1436] 反素数

标签:素数   最小   oid   include   cin   sim   约数个数   prim   cout   

原文地址:https://www.cnblogs.com/Mr-WolframsMgcBox/p/8524682.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!