ICode9

精准搜索请尝试: 精确搜索
首页 > 其他分享> 文章详细

CF1512G(数学,筛法,约数和)

2022-07-04 20:34:53  阅读:97  来源: 互联网

标签:约数 return 筛法 int CF1512G long include define


CF1512G(数学,筛法,约数和)

题意

求约数和为 \(c\) 的最小数字 \(n\)

思路

板子题,原题因为数据范围太小甚至能暴力(

因为约数和可以写成 \(\prod _{i = 1}^n(1+p_i+p_i^2+...+p_i^{q_i})\) 的形式,所以往下搜就行了。

中间要特判存在一个大质因子(\(p_i > \sqrt {n}\) )的情况。

被询问的答案拿map记一下不然会超时。

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<vector>
#include<set>
#include<queue>
#include<map>
#include<stack>
#include<string>
#include<random>
#include<iomanip>
#define yes puts("yes");
#define inf 0x3f3f3f3f
#define ll long long
#define linf 0x3f3f3f3f3f3f3f3f
#define ull unsigned long long
#define endl '\n'
#define int long long
#define rep(i,a,n) for(int i = a;i <= n;i++)
#define per(i,n,a) for(int i = n;i >= a;i--)
using namespace std;
mt19937 mrand(random_device{}());
int rnd(int x) { return mrand() % x;}
typedef pair<int,int> PII;
const int MAXN =10 + 2e5 ,mod=1e9 + 7;
bool st[MAXN];
vector<int> p;
void init() {
    const int N = 1e4 + 1;
    for(int i = 2;i <= N;i ++) {
        if(!st[i]) p.push_back(i);
        for(int j = 0;p[j] < N / i;j ++) {
            st[p[j] * i] = 1;
            if(i % p[j] == 0) break;
        }
    }
}
vector<int> ans;
bool isp(int x) {
    rep(i,2,x / i) if(x % i == 0) return false;
    return true;
}
void dfs(int id,int res,int s) {
    if(s == 1) {
        ans.push_back(res);
        return;
    }
    
    if(s - 1 >= p[id] && isp(s - 1)) ans.push_back(res * (s - 1));
    for(int i = id;p[i] <= s / p[i] && i < p.size();i ++) {
        for(int r = p[i],psum = 1 + p[i];psum <= s;r *= p[i], psum += r) {
            if(s % psum == 0) dfs(i + 1,res * r,s / psum);
        }
    }
}
map<int,int> mp;
void solve()
{    
    int n; cin >> n;
    if(mp.count(n)) {
        cout << mp[n] << endl;
        return;
    }
    ans.clear();
    dfs(0,1,n);
    if(ans.size()) {
        mp[n] = *min_element(ans.begin(),ans.end());
        cout << mp[n] << endl;
    }else {
        cout << -1 << endl;
        mp[n] = -1;
    }
}
signed main()
{
    ios::sync_with_stdio(false);cin.tie(0);cout.tie(0);
    init();
    int T;cin>>T;
    while(T--)
        solve();

    return 0;
}

标签:约数,return,筛法,int,CF1512G,long,include,define
来源: https://www.cnblogs.com/Mxrush/p/16444296.html

本站声明: 1. iCode9 技术分享网(下文简称本站)提供的所有内容,仅供技术学习、探讨和分享;
2. 关于本站的所有留言、评论、转载及引用,纯属内容发起人的个人观点,与本站观点和立场无关;
3. 关于本站的所有言论和文字,纯属内容发起人的个人观点,与本站观点和立场无关;
4. 本站文章均是网友提供,不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属;如您发现该文章侵犯了您的权益,可联系我们第一时间进行删除;
5. 本站为非盈利性的个人网站,所有内容不会用来进行牟利,也不会利用任何形式的广告来间接获益,纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。

专注分享技术,共同学习,共同进步。侵权联系[81616952@qq.com]

Copyright (C)ICode9.com, All Rights Reserved.

ICode9版权所有