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 ^ bint 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;}