UVa11137 - Ingenuous Cubrency（dp|递推）

题目链接

简介：
把一个数分成若干立方的和，求方案数

分析：
设计状态：f[i][j]表示考虑了前i个数之后，组成的数是j的方案数
因为n不超过10000 , 22^3>10000
所以最后的答案是f[21][n]

状态转移：
f[i][j+a*i^3]=sum{f[i-1][j]}

这样的转移是n^3的（枚举i，j，a）
虽然是时间承受范围之内，但是我们还有更优美的转移O(n^2)：

f[i][j]=f[i-1][j]+f[i][j-i^3]（j>=i*i*i）

f[i][j]=f[i-1][j] (j < i*i*i)

//这里写代码片
#include<cstdio>
#include<cstring>
#include<iostream>
#define ll long long

using namespace std;

ll f[22][10010];

void doit()
{
    int i,j;
    memset(f,0,sizeof(f));
    f[0][0]=1;
    for (i=1;i<=21;i++)
        for (j=0;j<=10000;j++)
        {
            if (j>=i*i*i) f[i][j]=f[i-1][j]+f[i][j-i*i*i];
            else f[i][j]=f[i-1][j];
        }           
}

int main()
{
    int n;
    doit();
    while (scanf("%d",&n)!=EOF)
        printf("%lld\n",f[21][n]);
    return 0;
}