Codeforces Round #338 div 2 D. Multipliers 数论

题意:给出一个数 x 的所有质因子,求这个数所有因子的乘积%1e9+7。

思路:一个因子a必然有另一个因子 b 与它对应使得ab=x。设总因子数为 k ,所以答案就应该是xk2。这个 k 是怎么求的呢?因为
x=a1k1a2k2ankn
k=(k1+1)(k2+1)(kn+1)

处理 xk2 要非常注意。在这里被坑了好久。
PS:c++11 auto真是好用的说。

http://codeforces.com/contest/615/problem/D

/*********************************************
    Problem : Codeforces
    Author  : NMfloat
    InkTime (c) NM . All Rights Reserved .
********************************************/

#include <map>
#include <set>
#include <queue>
#include <stack>
#include <cmath>
#include <ctime>
#include <cstdio>
#include <cstring>
#include <cstdlib>
#include <iostream>
#include <algorithm>

#define rep(i,a,b)  for(int i = (a) ; i <= (b) ; i ++) //遍历
#define rrep(i,a,b) for(int i = (b) ; i >= (a) ; i --) //反向遍历
#define repS(it,p) for(auto it = p.begin() ; it != p.end() ; it ++) //遍历一个STL容器
#define repE(p,u) for(Edge * p = G[u].first ; p ; p = p -> next) //遍历u所连接的点
#define cls(a,x)   memset(a,x,sizeof(a)) 
#define eps 1e-8

using namespace std;

//const int MOD = 1e9+7;
const int INF = 0x3f3f3f3f;
const int MAXN = 1e5+5;
const int MAXE = 2e5+5;
long long MOD = 1e9+7;

typedef long long LL;
typedef unsigned long long ULL;

int T,n,m,k;

int fx[] = {0,1,-1,0,0};
int fy[] = {0,0,0,-1,1};

map<int,int>mp;

LL ans ;

void input() {
    int tmp;
    ans = 1;
    mp.clear();
    rep(i,1,n) {
        scanf("%d",&tmp);
        mp[tmp]++ ; ans = (ans * (LL)tmp) % MOD;
    }
}

LL quick_pow(LL base , LL mult) {
    LL ret = 1;
    while(mult) {
        if(mult & 1) ret = (ret * base) % MOD;
        base = ( base * base ) % MOD;
        mult >>= 1;
    }
    return ret;
}

void solve() {
    LL mult = 1 , ss = 1;
    int ok = 0 ;
    auto pos = mp.begin();
    repS(it,mp) {
        if(it->second & 1) { ok = 1; pos = it;}
        ss = (ss * quick_pow(it->first,it->second/2)) % MOD;
    }
    if(ok) {
        repS(it,mp) {
            if(pos == it) {
                ans = quick_pow(ans,(it->second+1)/2);
            }
            else {
                ans = quick_pow(ans,it->second+1);
            }
        }
    }
    else {
        ans = ss;
        repS(it,mp) {
            ans = quick_pow(ans,it->second+1);
        }
    }
    printf("%I64d\n",ans);

}

int main(void) {
   // freopen("a.in","r",stdin);
    while(~scanf("%d",&n)) {
        input();
        solve();
    }
    return 0;
}
1、资源项目源码均已通过严格测试验证,保证能够正常运行;、 2项目问题、技术讨论,可以给博主私信或留言,博主看到后会第一时间与您进行沟通; 3、本项目比较适合计算机领域相关的毕业设计课题、课程作业等使用,尤其对于人工智能、计算机科学与技术等相关专业,更为适合; 4、下载使用后,可先查看README.md或论文文件(如有),本项目仅用作交流学习参考,请切勿用于商业用途。 5、资源来自互联网采集,如有侵权,私聊博主删除。 6、可私信博主看论文后选择购买源代码。 1、资源项目源码均已通过严格测试验证,保证能够正常运行; 2、项目问题、技术讨论,可以给博主私信或留言,博主看到后会第一时间与您进行沟通; 3、本项目比较适合计算机领域相关的毕业设计课题、课程作业等使用,尤其对于人工智能、计算机科学与技术等相关专业,更为适合; 4、下载使用后,可先查看README.md或论文文件(如有),本项目仅用作交流学习参考,请切勿用于商业用途。 5、资源来自互联网采集,如有侵权,私聊博主删除。 6、可私信博主看论文后选择购买源代码。 1、资源项目源码均已通过严格测试验证,保证能够正常运行;、 2项目问题、技术讨论,可以给博主私信或留言,博主看到后会第一时间与您进行沟通; 3、本项目比较适合计算机领域相关的毕业设计课题、课程作业等使用,尤其对于人工智能、计算机科学与技术等相关专业,更为适合; 4、下载使用后,可先查看README.md或论文文件(如有),本项目仅用作交流学习参考,请切勿用于商业用途。 5、资源来自互联网采集,如有侵权,私聊博主删除。 6、可私信博主看论文后选择购买源代码。
评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值