HDU 4574 Bombs (枚举+搜索)

转载请注明出处，谢谢http://blog.csdn.net/ACM_cxlove?viewmode=contents    by---cxlove

基本上想到了解法 ，但是一直在想[a,b]中乘积在[l,r]的组合有多少种怎么预处理。。。。

觉得当b - a + 1很大的时候，个数很少，大概可以承受n * m 的DP预处理。。。

结果结果就是个爆搜嘛。。。。不过加了些优化 。。。

取的数中最小是2，2 ^ 30 > 1e9 。所以个数不多

首先枚举长度，然后爆搜得到个数，确定长度之后，依旧是按位确定。

枚举每一位的数，继续爆搜判断是否满足第k大。。。。

有一些优化 ：

1、当个数为1的时候，就不用枚举了，分析下上下界，O(1)可以得到结果

2、考虑剩下的c - 1个数都取上界，然后利用区间的下界可以得到当前数的下界，考虑剩下的c - 1个数都取下界，然后 利用区间的上界可以得到当前数的上界，这个剪枝很重要吧。

#include <iostream>
#include <cstdio>
#include <cmath>
#include <cstring>
#include <queue>
#include <vector>
#define mp(a , b) make_pair(a , b)
using namespace std;
typedef long long LL;
const LL inf = 1000000000LL;
int a , b , l , r , k;
// a ^ b
int pow (int a , int b) {
    int ret = 1;
    for (int i = 0 ; i < b ; i ++) {
        ret = ret * a;
        if (ret > inf || ret <= 0)
            return inf + 1;
    }
    return ret;
}
int gao (int c , int a , int b , int l , int r , int k) {
    if (c == 0) return l <= r;
    if (c == 1) {
        if (l > r) return 0;
        return max (0 , min(r , b) - max (l , a) + 1);
    }
    //上下界
    int down = max(a , l / pow(b , c - 1));
    int up = min (b , r / pow(a , c - 1));
    int ret = 0;
    for (int i = down ; i <= up ; i ++) {
        ret += gao (c - 1 , i , b , (l + i - 1) / i , r / i , k - ret);
        if (ret > k) return ret;  //剪枝
    }
    return ret;
}
int ret , ans[50];
void fuck (int c , int a , int b , int l , int r , int k) {
    if (c == 0) return;
    if (c == 1) {
        int num = k + max(l , a) - 1;
        ret = ret * num;
        ans[c] = num;
        return ;
    }
    // 上下界
    int down = max(a , l / pow(b , c - 1));
    int up = min (b , r / pow(a , c - 1));
    for (int i = down ; i <= up ; i ++) {
        int cnt = gao (c - 1 , i , b , (l + i - 1) / i , r / i , k);
        if (k > cnt) {
            k -= cnt;
            continue;
        }
        ret = ret * i;
        ans[c] = i;
        fuck (c - 1 , i , b , (l + i - 1) / i , r / i , k);
        return ;
    }
}
int main () {
    #ifndef ONLINE_JUDGE
        freopen ("input.txt" , "r" , stdin);
    #endif
    int t , cas = 0;
    scanf ("%d" , &t);
    while (t --) {
        bool ok = false;
        scanf ("%d %d %d %d %d" , &a , &b , &l , &r , &k);
        printf ("Case #%d: " , ++ cas);
        //枚举长度
        for (int i = 1 ; i <= 30 ; i ++) {
            int cnt = gao (i , a , b , l , r , k);
            if (k > cnt) {
                k -= cnt;
                continue;
            }
            ret = 1;
            fuck (i , a , b , l , r , k);
            printf ("%d\n" , ret);
            for (int j = i ; j >= 1 ; j --) {
                printf ("%d%c" , ans[j] , j == 1 ? '\n' : ' ');
            }
            ok = true;
            break;
        }
        if (!ok) puts("-1");
    }
    return 0;
}