poj 1945 Power Hungry Cows 启发式搜索

这个题目用了启发式搜索来做，hash判重，优先队列储存节点。启发式函数就是大的数字x*2^k>n，即为k。

x为节点中较大值，y为节点中较小值

然后后面就是一些剪枝了，当x>2*n不能得出最优解。当y==0，不能得出最优解。当x==y，不能得出最优解。

还有一个比较关键的剪枝，gcd(x,y)不能整除n时，得不出结果。

#include

#include

#include

#include

#include

using namespace std;

const int maxn=160000+10;

int n;

int d[22222];

int time;

struct node

{

    int x,y,h,cost;

    bool operator <(const node &xx) const

    {

        return(cost>xx.cost);

    }

};

priority_queue que;

int hash[maxn],lon;

struct

{

    int x,y,cost,next;

}data[11111111];

void insert(struct node *p)

{

    inttmp=p->x+p->y;

    data[++lon].x=p->x;

    data[lon].y=p->y;

    data[lon].cost=p->h;

    data[lon].next=hash[tmp];

    hash[tmp]=lon;

}

int gcd(int a,int b)

{

    int tmp;

    while(a%b)

    {

        tmp=a%b;

        a=b;

        b=tmp;

    }

    return(b);

}

int check(int x,int y,int h)

{

    if(x<=0)return(0);

    if(y!=0)

    if(n%gcd(x,y)!=0)

    {

        return(0);

    }

    if(x==n||y==n)

    {

        printf("%d\n",h);

        exit(0);

    }

    if(x

    swap(x,y);

    if(x>n&&y>n) return(0);

    if(x>n*2)return(0);

    if(x>n&&y==0)return(0);

    if(x==y)return(0);

    for(intk=hash[x+y];k!=-1;k=data[k].next)

    {

        if(data[k].x==x&&data[k].y==y)

        {

            if(data[k].cost>h)

            {

                data[k].cost=h;

                return(1);

            }

            else

            return(0);

        }

    }

    return(1);

}

int cal(int x,int y)

{

    if(x==n)

    return (0);

    if(x

    return(d[x]);

    if(x>n)

    {

        int k=1;

        int tmp=x-n;

        while(y

        {

            k++;

            y<<=1;

        }

        return(k);

    }

}

void work(int x,int y,int h)

{

    if(x

    swap(x,y);

    struct node p;

    p.x=x;

    p.y=y;

    p.h=h;

    p.cost=h+cal(x,y);

    insert(&p);

    que.push(p);

}

int bfs()

{

    while(!que.empty())

    que.pop();

    memset(hash,-1,sizeof(hash));

    lon=0;

    struct node p;

    p.x=1;

    p.y=0;

    p.h=0;

    p.cost=p.h+cal(1,0);

    insert(&p);

    que.push(p);

    while(!que.empty())

    {

        intx=que.top().x,y=que.top().y,h=que.top().h+1,cost=que.top().cost;

        que.pop();

        if(x==n||y==n)

        return(cost);

        if(check(x+y,y,h))

        work(x+y,y,h);

        if(check(x,x+y,h))

        work(x,x+y,h);

        if(check(x+x,y,h))

        work(x+x,y,h);

        if(check(x,x+x,h))

        work(x,x+x,h);

        if(check(y+y,y,h))

        work(y+y,y,h);

        if(check(y+y,x,h))

        work(y+y,x,h);

        if(check(x-y,y,h))

        work(x-y,y,h);

        if(check(x,x-y,h))

        work(x,x-y,h);

    }

}

int main()

{

    scanf("%d",&n);

    for(inti=n,j,t=0;i>=1;i=j)

    {

        j=i/2,t++;

        for(int k=i-1;k>j;k--)

        d[k]=t+1;

        d[j]=t;

    }

    bfs();

    return 0;

}