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

hdu 1215(因子和)

时间:2016-05-23 21:10:17      阅读:122      评论:0      收藏:0      [点我收藏+]

标签:

七夕节

Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 65536/32768 K (Java/Others)
Total Submission(s): 40119    Accepted Submission(s): 12613


Problem Description
七夕节那天,月老来到数字王国,他在城门上贴了一张告示,并且和数字王国的人们说:"你们想知道你们的另一半是谁吗?那就按照告示上的方法去找吧!"
人们纷纷来到告示前,都想知道谁才是自己的另一半.告示如下:

技术分享


数字N的因子就是所有比N小又能被N整除的所有正整数,如12的因子有1,2,3,4,6.
你想知道你的另一半吗?
 

 

Input
输入数据的第一行是一个数字T(1<=T<=500000),它表明测试数据的组数.然后是T组测试数据,每组测试数据只有一个数字N(1<=N<=500000).
 

 

Output
对于每组测试数据,请输出一个代表输入数据N的另一半的编号.
 

 

Sample Input
3 2 10 20
 

 

Sample Output
1 8 22
n的因子和为s(n) 令g(p, e) = (p^(e+1) - 1) / (p-1),则s(n) = g(p1, e1) * g(p2, e2) * ... * g(pk, ek)
要分类讨论一下。
#include <stdio.h>
#include <algorithm>
#include <string.h>
#include <math.h>
#include <queue>
using namespace std;
typedef long long LL;
const int N = 10000;
struct Node{
    int p,num;
}node[N];
LL g(int p,int e){ ///g(p, e) = (p^(e+1) - 1) / (p-1)
    LL sum = 1;
    for(int i=0;i<e+1;i++){
        sum*=p;
    }
    sum--;
    return sum/(p-1);
}
int main()
{
    int tcase;
    scanf("%d",&tcase);
    while(tcase--){
        int n;
        scanf("%d",&n);
        if(n==1){
            printf("0\n");
            continue;
        }
        int m = n;
        int id = 0;
        for(int i=2;i*i<=n;i++){
            if(n%i==0){
                node[id].p = i;
                node[id].num = 0;
                while(n%i==0) {
                    n/=i;
                    node[id].num++;
                }
                id++;
            }
        }
        if(n>1&&m!=n) {node[id].p =n,node[id++].num=1;}
        LL sum = 1;
        for(int i=0;i<id;i++){
            sum*=g(node[i].p,node[i].num);
        }
        if(id==0) printf("%d\n",sum);
        else printf("%d\n",sum-m);
    }
    return 0;
}

 

hdu 1215(因子和)

标签:

原文地址:http://www.cnblogs.com/liyinggang/p/5521284.html

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