HDU4394 Digital Square【分支限界法BnB】

Digital Square

Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Others)
Total Submission(s): 2062    Accepted Submission(s): 830


Problem Description
Given an integer N,you should come up with the minimumnonnegative integer M.M meets the follow condition: M2%10x=N (x=0,1,2,3....)
 

Input
The first line has an integer T( T< = 1000), the number of test cases.
For each case, each line contains one integer N(0<= N <=109), indicating the given number.
 

Output
For each case output the answer if it exists, otherwise print “None”.
 

Sample Input
3 3 21 25
 

Sample Output
None 11 5
 

Source


问题链接HDU4394 Digital Square

题意简述

  输入测试用例数量t,输入t个正整数n,求最小的m,满足m^2%10^x=n,其中k=0,1,2,...。

问题分析

  x是不定的,用暴力法试探m是不现实的。有关模除%的问题,可以用从低位开始逐步试探的方法,即先试探个位,然后十位、百位、千位等等。已知n的情况下,m的低位是有限的。例如n=21,那么m的个位只能是1或9,其他的数字是无法满足给定条件的。解决本问题需要明确以下几点:

  1.对于m,从最低位开始试探,m的每位的数字可以是0-9,这时匹配n的最低位;

  2.试探完1位再试探2位、3位和4位等等,分别匹配n的最低1位、2位、3位和4位等等

  3.m起始值:0,1,2,3,......;

  4.显式剪枝条件:低位不匹配;

  5.LC函数(优先函数):为了找到最小的满足条件的m,低位匹配的m中,值最小的节点优先展开。

程序说明

  节点node的成员变量说明如下:

  curr:当前低位匹配的m值;

  len:当前已经匹配的位数;

  digit:当前位正在试探的数字(0-9)。


题记

  分支限界法有三种策略,分别是FIFO、LIFO和LC(least cost)。BFS属于分支限界法的一种,通常采用FIFO策略,采用LIFO策略的情况比较少见,因为多数情况下这两种策略效果几乎相同。分支限界法采用LC策略时,通常用BFS+优先队列来实现。



AC的C++语言程序如下:

/* HDU4394 Digital Square */

#include <iostream>
#include <queue>

using namespace std;

typedef unsigned long long ULL;

const int MAXN = 18;
const int MAXDIGIT = 9;

struct node {
    ULL curr;
    int len, digit;
    node(){}
    node(ULL c, int l, int n):curr(c), len(l), digit(n){}
    bool operator<(const node n) const {
        return curr > n.curr;
    }
};

node start;

int n;
ULL fact[18];
int nlen;
ULL ans;

void bfs()
{
    nlen = 0;
    ULL temp = n;
    while(temp) {
        nlen++;
        temp /= 10;
    }

    ans = 0;
    priority_queue<node> q;
    q.push(node(0, 0, 0));

    while(!q.empty()) {
        node top = q.top();
        q.pop();

        if(top.len == nlen) {
            ans = top.curr;
            break;
        } else if(top.digit != MAXDIGIT)
            q.push(node(top.curr, top.len, top.digit+1));

        node v;
        v.curr = top.curr + fact[top.len] * top.digit;
        v.len = top.len + 1;
        v.digit = 0;

        if(v.curr * v.curr % fact[v.len] == n % fact[v.len])
            q.push(v);
    }
}

int main()
{
    int t;

    fact[0] = 1;
    for(int i=1; i<MAXN; i++)
        fact[i] = fact[i-1] * 10;

    cin >> t;
    while(t--) {
        cin >> n;

        bfs();

        if(ans)
            printf("%llu\n", ans);
        else
            printf("None\n");
    }

    return 0;
}


阅读更多
想对作者说点什么?

博主推荐

换一批

没有更多推荐了,返回首页