hdu 1028 Ignatius and the Princess III

转载链接

加了点注释方便自己看

#include <cstring>
#include <iostream>
using namespace std;

int cur[123];//存中间结果
int pol[123][123];//表示多项式，里面存的是项的系数。pol[i]=第i个多项式，pol[i][j]=第i个多项式里，指数为j的项的系数
int ans[123] = {1};//存最后的多项式乘积结果，也就是可能的组合，都初始化为1

int main()
{
    // 初始化各个多项式
    for (int i = 1; i <= 120; ++i)
    {
        for (int j = 0; j <= 120; j += i)
        {//注意j+的是i，比如当i等于2时，只有偶次项会等于1 => X^0 + 0 + X^2 + 0 + X^4
            pol[i][j] = 1;
        }
    }

    // ans依次乘上所有多项式，ans暂时表示上一次的2个多项式的乘积结果，都乘完了以后，存的就是最后结果
    // 比如1*2*3，ans先存1*2的值，再用ans乘3
    //              ans                                 pol[1]                                pol[2]...
    //（X^0 + X^1 + X^2 + X^3 + X^4 + X^5）*（X^0 + X^1 + X^2 + X^3 + X^4 + X^5）*（X^0 + 0 + X^2 + 0 + X^4 + 0）...
    for (int i = 1; i <= 120; ++i) 
    {
        //先用上一次的乘积结果ans里的 X^0乘第一个多项式，再用X^1乘第一个多项式，以此类推。结果存回ans里，这样循环乘i个多项式，就完成了多项式乘积
        for (int j = 0; j <= 120; ++j)
        {
            for (int k = 0; j + k <= 120; ++k)// j+K<=5防止越界,只要算到5次项就好，多了没必要
            {                                    
                //ans[j]=上一次乘积的结果里，指数为j的项的系数，pol[i][k] = 第i个多项式里，指数为k的项的系数
                //cur[j+k]=指数为(j+k)这个项的系数 = ans[j]与pol[i][k]这2个项的系数相乘。（j+k = 2个项相乘后的指数值）
                cur[j + k] += ans[j] * pol[i][k];
            }
        }

        memcpy(ans, cur, sizeof(cur));//存结果
        memset(cur, 0, sizeof(cur));//置零
    }

    int q;
    while (cin >> q)
        cout << ans[q] << endl;
    return 0;
}