HDU4447 Yuanfang, What Do You Think?

2012 Asia JinHua Regional Contest
Problem F. Yuanfang, What Do You Think?
官方解题报告：
手工或写暴力程序计算较小的n，可观察到(x^n-1)的分解式中含有(x^ni-1)的分解式。
其中ni是n的约数（n本身除外）。
除了(x^ni-1)的分解式外，(x^n-1)还含有一项自己独有的分解式，记为P(n)。
于是我们就得到了(x^n-1)的分解：(x^n-1)=P(n1)P(n2)P(n3)...P(n)
其中，n1,n2,n3,...是n的约数（n本身除外）。
从n=1算起，P(1)=x-1。
对于所有的n，P(n)均可通过前面的结果算得：P(n)=(x^n-1)/P(n1)/P(n2)/P(n3)/...
最后对P(n1)P(n2)P(n3)...P(n)排序即可。
时间复杂度：O(n^2*logn)。

[cpp]
#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;
const int MAXN = 1111;

inline void printX(int x, int v, bool show)
{
    if(0 == v)
    {
        return;
    }
    if(show)
    {
        if(v > 0)
        {
            printf("+");
        }
        else
        {
            printf("-");
        }
    }
    if(0 == x)
    {
        printf("1");
    }
    else
    {
        if(v != 1 && v != -1)
        {
            printf("%d", abs(v));
        }
        printf("x");
        if(x > 1)
        {
            printf("^%d", x);
        }
    }
}

struct Polynomial
{
    int n;
    int value[MAXN];

    void output() const
    {
        printX(n, value[n], false);
        for(int i=n-1;i>=0;--i)
        {
            printX(i, value[i], true);
        }
    }
} intrinsic[MAXN], sets[MAXN];
int setNumber;
bool visit[MAXN];

bool operator < (const Polynomial &a, const Polynomial &b)
{
    if(a.n == b.n)
    {
        for(int i=a.n;i>=0;--i)
        {
            if(abs(a.value[i]) == abs(b.value[i]))
            {
                if(a.value[i] != b.value[i])
                {
                    return a.value[i] < b.value[i];
                }
            }
            else
            {
                return abs(a.value[i]) < abs(b.value[i]);
            }
        }
    }
    return a.n < b.n;
}

Polynomial operator / (Polynomial a, const Polynomial &b)
{
    Polynomial c;
    c.n = a.n - b.n;
    for(int i=c.n;i>=0;--i)
    {
        if(a.n - b.n < i)
        {
            c.value[i] = 0;
            continue;
        }
        c.value[i] = a.value[a.n] / b.value[b.n];
        int d = a.n - b.n;
        a.n -= b.n;
        for(int j=b.n;j>=0;--j)
        {
            a.value[j + d] -= b.value[j] * c.value[i];
            if(a.value[j + d])
            {
                a.n = max(a.n, j + d);
            }
        }
    }
    return c;
}

void makeSet(int n)
{
    setNumber = 0;
    for(int i=1;i<=n;++i)
    {
        if(0 == n % i)
        {
            if(!visit[i])
            {
                visit[i] = true;
                intrinsic[i].n = i;
                intrinsic[i].value[i] = 1;
                intrinsic[i].value[0] = -1;
                for(int j=1;j<i;++j)
          &

补充：软件开发 , C++ ,