【USACO】奶牛求幂

题目

【问题描述】

　　约翰的奶牛门想要快速计算整数的P次幂，他需要你的帮助。因为计算极大数的幂，所以他们同一时间仅能使用2个存储器，每个存储器可记录某个结果值。

　　第一件工作是初始化存储器内的值：一个为底数x，另一个为1。奶牛可以相乘法或相除法2个存储器中的值，并把结果存在其中某个存储器内，但所有存储器的结果必须是整数，例如他们想计算x^31，一种计算方法是：
　　
　　　
　　

因此，x^31可以通过6次计算得出。给出要计算幂次P，要求求出最少需要几次计算。

【输入格式】

　　仅一行一个正整数P。

【输出格式】

　　仅一个整数：最少计算次数。

【输入样例】

31

【输出样例】

6

【数据范围】

对于30%的数据，P<=10；
对于60%的数据，P<=500；
对于全部的数据，P<=20000；

思路

模型提炼：
有a=1,b=0，经过最小次数加减得到p
这种感觉不到明显算法的题目就先考虑暴搜，因为这个步数不定，就考虑DFSID和BFS，BFS的状态开不下，且无法剪枝，决定用DFSID

剪枝
1、log求出最大情况可否达到p，然后我不想用log就倍增了一下做出来了
2、P%gcd(a,b)!=0则无解

剪枝2证明：
因为a+-b=（a0+b0）*gcd(a,b)=k*gcd(a,b),k为整数，所以p必须能被gcd(a,b) 整除

代码

#include<cmath>
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm> 
using namespace std;
int k,p,qkpow[40];
void ready()
{
    qkpow[0]=1;
    for(int i=1;i<=20;i++)
        qkpow[i]=qkpow[i-1]*2;
}
int SGP(int x,int i)
{
    if(i>16)
        return p;
    return x*qkpow[i];
}
int gcd(int a,int b)
{
    return b==0?a:gcd(b,a%b);
}
bool DFSID(int i,int a,int b)
{
    if(i>k)
    {
        if(a==p||b==p)
            return 1;
        return 0;
    }
    if(a!=0&&b!=0&&p%gcd(a,b)!=0)
        return 0;
    int A=max(a,b);
    if(SGP(A,k-i+1)<p)return 0;
    if(DFSID(i+1,a,a+a))return 1;
    if(DFSID(i+1,a+a,b))return 1;
    if(DFSID(i+1,a,a+b))return 1;
    if(DFSID(i+1,a+b,b))return 1;
    if(DFSID(i+1,b+b,b))return 1;
    if(DFSID(i+1,a,b+b))return 1;
    if(DFSID(i+1,a,abs(a-b)))return 1;
    if(DFSID(i+1,abs(a-b),b))return 1;
    return 0;
}
int main()
{
    //freopen("in.txt","r",stdin);
    scanf("%d",&p);
    ready();
    for(k=0;;k++)
    {
        if(DFSID(1,1,0))
        {
            printf("%d\n",k);
            break;
        }
    }
    return 0;
}

这道题比较重要的地方就是gcd的剪枝，感觉很巧妙