zoukankan      html  css  js  c++  java
  • hdu 5428 The Factor(数学)

    Problem Description
    There is a sequence of n positive integers. Fancycoder is addicted to learn their product, but this product may be extremely huge! However, it is lucky that FancyCoder only needs to find out one factor of this huge product: the smallest factor that contains more than 2 factors(including itself; i.e. 4 has 3 factors so that it is a qualified factor). You need to find it out and print it. As we know, there may be none of such factors; in this occasion, please print -1 instead. 
     
    Input
    The first line contains one integer T (1T15), which represents the number of testcases. 
    For each testcase, there are two lines:
    1. The first line contains one integer denoting the value of n (1n100).
    2. The second line contains n integers a1,,an (1a1,,an2×109), which denote these n positive integers. 
     
    Output
    Print T answers in T lines.
     
    Sample Input
    2 3 1 2 3 5 6 6 6 6 6
     
    Sample Output
    6 4
     
    Source
     

    对于每一个数字,它有用的部分其实只有它的所有质因子(包括相等的)。求出所有数的所有质因子中最小的两个,相乘就是答案。如果所有数字的质因子个数不到两个,那么就是无解。时间复杂度

    O(n*sqrt(a))O(nsqrt(a))

     1 #include<iostream>
     2 #include<cstdio>
     3 #include<cstring>
     4 #include<set>
     5 #include<cmath>
     6 #include<algorithm>
     7 using namespace std;
     8 #define ll long long 
     9 multiset<ll> s;
    10 multiset<ll>::iterator it,it1;
    11 int main()
    12 {
    13     int t;
    14     scanf("%d",&t);
    15     while(t--){
    16         s.clear();
    17         ll n;
    18         scanf("%I64d",&n);
    19         for(ll i=0;i<n;i++){
    20             ll m;
    21             scanf("%I64d",&m);
    22             for(ll j=2;j*j<=m;j++){
    23                 if(m%j==0){
    24                     while(m%j==0){
    25                         s.insert(j);
    26                         m/=j;
    27                     }
    28                 }
    29             }
    30             if(m>1) s.insert(m);
    31         }
    32         ll size=s.size();
    33         
    34         if(size<2) printf("-1
    ");
    35         else{
    36             it=s.begin();
    37             ll ans1=(*it);
    38             it++;
    39             ll ans2=(*it);
    40             printf("%I64d
    ",ans1*ans2);
    41         } 
    42     }
    43     return 0;
    44 }
    View Code
  • 相关阅读:
    SQL EXPLAIN优化详解
    2019数据智能算法大赛赛后复盘
    K-D树详解
    点云配准相关
    pandas速查手册(中文版)
    Windows CMD命令大全
    获取时间戳
    SQL 字符串转成临时表
    delphi cxgrid明细新增第三行的报错【Key violation】
    cxgrid列的Properties(cxPopupEditPopup)的关闭方法
  • 原文地址:https://www.cnblogs.com/UniqueColor/p/4784042.html
Copyright © 2011-2022 走看看