HDU4516 威威猫系列故事——因式分解 多项式分解

题意：给定一个多项式，对其进行因式分解。

解法：由于多项式每一项系数绝对值不超过1000，由于最后解的形式为(x-a)(x-b)(x-c)(x-d)(x-e)其中a*b*c*d*e一定是最后的常数项系数，因此a, b, c, d, e的取值范围都在[-1000, 1000]内，因此枚举所有的根，剩下的就是重根的时候该怎么办？一个解决办法就是对原多项式进行求导，如果一个根t是f(x)的K重根的话，那么t一定是f(x)'的K-1重根。该题的字符串处理我没写好，后面调了很久。还有就是由于有5次方存在，因此代入时使用long long计算。

代码如下：

#include <cstdlib>
#include <cstring>
#include <cstdio>
#include <algorithm>
#include <cctype>
#include <vector>
using namespace std;

char str[500];
int seq[10];

long long _pow(int a, int b) {
    long long ret = 1;
    for (int i = 0; i < b; ++i) {
        ret *= a;
    }
    return ret;
}

int jiechen[10] = {1, 1, 2, 6, 24, 120};

void qiudao(int *rec, int k) {
    for (int i = k; i <= 5; ++i) {
        rec[i-k] = jiechen[i] / jiechen[i-k] * seq[i];
    }
}

bool judge(int rec[], int x) {
    long long sum = 0;
    for (int i = 0; i <= 5; ++i) {
        sum += 1LL * rec[i] * _pow(x, i);
    }
    return sum == 0;
}

void gao(char ts[]) {
    int len = strlen(ts);
    int p = -1, a, b;
    for (int i = 0; i < len; ++i) {
        if (ts[i] == 'x') {
            if (isdigit(ts[i-1])) { 
                ts[i] = '\0';
            } else {
                ts[i] = '1';
            }
            p = -2;
        } else if (ts[i] == '^') {
            ts[i] = '\0';
            p = i+1;    
        }
    }
    a = atoi(ts);
    if (!a && p != -1) a = 1;
    if (p == -1) {
        b = 0;
    } else if (p == -2) {
        b = 1;
    }else {
        b = atoi(ts+p);
    }
    seq[b] += a;
}

void solve() {
    vector<int>v;
    int cnt = 0;
    memset(seq, 0, sizeof (seq));
    char ts[50], *p;
    p = strtok(str, "+");
    while (p) {
        strcpy(ts, p);
        gao(ts);
        p = strtok(NULL, "+");
    }
    for (int i = 5; i >= 0; --i) {
        if (seq[i] != 0) {
            cnt = i;
            break;
        }
    }
    for (int i = -1000; i <= 1000; ++i) { 
        for (int j = 0; j < cnt; ++j) {
            int rec[10]    = {0};
            qiudao(rec, j);
            if (judge(rec, i)) {
                v.push_back(i);
            } else {
                break;
            }
        }
    }
    //x^4-x^2
    //x^4-7x^3+18x^2-20x+8
    //x^3-13x^2+55x-75
    //x^2+5x^2-6x^2+x^2+2x-20x+30x-10x+8-7
    //x^5-10x^4+39x^3-74x^2+68x-24
    //以上都是能够分解的式子 
    if (v.size() != cnt || seq[cnt] != 1 || cnt == 0) {
        printf("-1\n");
    } else {
        sort(v.begin(), v.end());
        for (int i = v.size()-1; i >= 0; --i) {
            if (v[i] < 0) {
                printf("(x+%d)", -v[i]);
            } else if (v[i] == 0) {
                printf("x");    
            } else {
                printf("(x-%d)", v[i]);    
            }
        }
        puts("");
    }
}

int main() {
    int T, ca = 0;
    scanf("%d", &T);
    while (T--) {
        scanf("%s", str);
        int len = strlen(str);
        for (int i = 0; i < len; ++i) {
            if (str[i] == '-') {
                for (int j = len-1; j >= i; --j) {
                    str[j+1] = str[j];
                }
                str[i] = '+';
                len += 1;
                ++i;
                str[len] = '\0';
            }
            
        }
        printf("Case #%d: ", ++ca);
        solve();
    }
    return 0;
}