首页 > 其他 > 详细

HDU 5778 abs 数学

时间:2017-03-04 21:59:28      阅读:229      评论:0      收藏:0      [点我收藏+]

http://acm.hdu.edu.cn/showproblem.php?pid=5778

这题的意思就是找离x最近的一个数y,且y是一个完全平方数,还是所有质因子都只能出现两次的完全平方数

 

一开始的思路是直接枚举这个差值,然后去两边找,val - res和val + res找,然后超时了。

其实也很正常,因为两个完全平方数的间隔实在太大了。中间有很多的数字,

那么把思路换一下,

题意是找a^2   <=  x  <=  b^2这样子的东西,那么可以同时开方,a <= sqrt(x) <= b

这样找到的结果是一样的,然后就相当于找一个没有质因子出现过两次的数字了。

技术分享
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <cmath>
#include <algorithm>
#include <assert.h>
#define IOS ios::sync_with_stdio(false)
using namespace std;
#define inf (0x3f3f3f3f)
typedef long long int LL;


#include <iostream>
#include <sstream>
#include <vector>
#include <set>
#include <map>
#include <queue>
#include <string>
#include <bitset>
const int maxn = 1e6 + 20;
int prime[maxn], total;
bool check[maxn];
void initprime() {
    for (int i = 2; i <= maxn - 20; ++i) {
        if (!check[i]) {
            prime[++total] = i;
        }
        for (int j = 1; j <= total; ++j) {
            if (i * prime[j] > maxn - 20) break;
            check[i * prime[j]] = 1;
            if (i % prime[j] == 0) break;
        }
    }
}
bool isok(LL val) {
    if (val < 2) return false;
    bool flag = false;
    for (int i = 1; i <= total; ++i) {
        if (val < prime[i]) break;
        if (val % prime[i] == 0) {
            val /= prime[i];
            if (val % prime[i] == 0) return false;
            flag = true;
//            while (val % prime[i] == 0) {
//                val /= prime[i];
//                has++;
//                if (has == 3) return false;
//                flag = true;
//            }
        }
    }
//    if (val != 1) return false;
    return true;
}
void work() { // big prime * big prime  TLE   && 2^6 not ok
    LL val;
//    scanf("%I64d", &val);
    cin >> val;
    LL tpos = val;
    val = (LL)sqrt(val * 1.0);
    LL mid = 2e18;
    if (isok(val)) {
        mid = val * val;
    }
    LL big = 2e18, small = 2e18;
    for (int res = 1; ; ++res) {
        if (isok(val + res)) {
            big = (val + res) * (val + res);
            break;
        }
    }
    for (int res = 1; ; ++res) {
        if (val - res < 2) break;
        if (isok(val - res)) {
            small = (val - res) * (val -res);
            break;
        }
    }
    LL ans = min(abs(tpos - big), abs(tpos - small));
    ans = min(ans, abs(tpos - mid));
//    cout << big << " " << mid << " " << small << endl;
    cout << ans << endl;
}
int main() {
#ifdef local
    freopen("data.txt", "r", stdin);
//    freopen("data.txt", "w", stdout);
#endif
//    IOS;
    initprime();
//    cout << isok(33) << endl;
    int t;
    scanf("%d", &t);
    while (t--) work();
    return 0;
}
View Code

 

HDU 5778 abs 数学

原文:http://www.cnblogs.com/liuweimingcprogram/p/6502669.html

(0)
(0)
   
举报
评论 一句话评论(0
关于我们 - 联系我们 - 留言反馈 - 联系我们:wmxa8@hotmail.com
© 2014 bubuko.com 版权所有
打开技术之扣,分享程序人生!