这道题就是个简单的深搜题,不过需要剪枝。
首先addition chains在求幂的时候用到。比如a^15=a*a^14=a*(a^7)^2=a*(a^6*a)^2=...
所以至少有如下一种分解方法:如果n是奇数,分解(n-1);如果n是偶数,分解(n/2)。
这种分解方法可以作为初始上界。之后每次求得一个比这个值小的解时候,更新上界。
其次,避免过多的重复计算。新求出一个数之后,只需要判断这个数和它之前的数的组合(包括它自己),而不需考虑前面的两两组合(因为之前已经求过了)。
代码如下:
#include <iostream>
using namespace std;
void additionchain(int intermids[], int len, int n, int &limit, int res[])
{
if(len >= limit)
return;
for(int i = len - 1; i >= 0; --i)
{
intermids[len] = intermids[len - 1] + intermids[i];
if(intermids[len] == n)
{
if(len + 1 <= limit)
{
limit = len + 1;
for(int i = 0; i <= len; ++i)
{
res[i] = intermids[i];
}
}
return;
}
if(intermids[len] < n)
{
additionchain(intermids, len + 1, n, limit, res);
}
}
}
int main()
{
int n;
int intermids[100];
int res[100];
int len = 0;
while(true)
{
cin>>n;
if(n == 0)
break;
else if(n == 1)
cout<<"1"<<endl;
else if(n == 2)
cout<<"1 2"<<endl;
else
{
intermids[0] = 1;
intermids[1] = 2;
len = 2;
int limit = 0;
int tmp = n;
while(tmp)
{
if(tmp % 2)
{
tmp--;
}
else
{
tmp >>= 1;
}
limit++;
}
additionchain(intermids, len, n, limit, res);
for(int i = 0; i < limit; ++i)
{
cout<<res[i];
if(i != limit - 1)
cout<<" ";
else
cout<<endl;
}
}
}
return 0;
}