PAT-天梯赛习题集-L2-018-多项式A除以B

ACM模版

描述

题解

原本我以为这个是今年 CCCC 的初赛中最难得题了，毕竟 AC 率是最低的了，可是我发现这个其实真的不是最难的，最难的是那个多标签最短路，为什么这么说，因为我到现在那道题还是19分没有拿全，不知道 bug 出在哪儿了，而这道题，实际上不是代码难写，也不是思路复杂，而是我根本不知道什么是多项式除法，所以比赛时大眼瞪小眼瞪了半个多小时却没有看懂样例，一旦知道什么是多项式除法以及具体怎么个解法儿，那么一切都迎刃而解了，模糊的讲，多项式除法实际上和一般的除法的解法儿一样，不断除一直到无法除了之后剩下的就是余，具体的我也不知道怎么说，参考里提供一个多项式除法的视频链接，看懂了这个链接，一切都不是难题了，只消一个递归，一个多项式减法，就 GG 了。

好了，我继续去找那个多标签最短路问题的 bug 了，心塞啊~~~

代码

#include <iostream>
#include <cstdio>
#include <cmath>
#include <list>

using namespace std;

const double ESP = 0.04;

list<pair<int, double> > A, B, Q, R, T;

int E;
double C;
int NA, NB, NQ = 0, NR = 0;

void sub()
{
    T.clear();
    auto i = A.begin(), j = R.begin();
    for (; i != A.end(); i++)
    {
        if (j == R.end())
        {
            break;
        }
        for (; j != R.end(); j++)
        {
            if ((*i).first > (*j).first)
            {
                T.push_back(make_pair((*i).first, (*i).second));
                break;
            }
            else if ((*i).first < (*j).first)
            {
                T.push_back(make_pair((*j).first, -(*j).second));
                continue;
            }
            else
            {
                T.push_back(make_pair((*i).first, (*i).second - (*j).second));
                j++;
                break;
            }
        }
    }
    for (; i != A.end(); i++)
    {
        T.push_back(make_pair((*i).first, (*i).second));
    }
    for (; j != R.end(); j++)
    {
        T.push_back(make_pair((*j).first, -(*j).second));
    }
    swap(A, T);
}

void solve()
{
    if (A.front().first < E)
    {
        return ;
    }
    int tempE = A.front().first - E;
    double tempC = A.front().second / C;
    if (fabs(tempC) > ESP)
    {
        Q.push_back(make_pair(tempE, tempC));
        NQ++;
    }
    A.erase(A.begin());
    R.clear();
    for (auto i = B.begin(); i != B.end(); i++)
    {
        R.push_back(make_pair(tempE + (*i).first, tempC * (*i).second));
    }

    //  A = A - R
    sub();

    solve();
}

int main(int argc, const char * argv[])
{
//    freopen("/Users/zyj/Desktop/input.txt", "r", stdin);

    int e;
    double c;

    cin >> NA;
    for (int i = 0; i < NA; i++)
    {
        scanf("%d%lf", &e, &c);
        A.push_back(make_pair(e, c));
    }
    cin >> NB;
    for (int i = 0; i < NB; i++)
    {
        scanf("%d%lf", &e, &c);
        B.push_back(make_pair(e, c));
    }
    E = B.front().first;
    C = B.front().second;
    B.erase(B.begin());

    solve();

    if (NQ == 0)
    {
        cout << "0 0 0.0\n";
    }
    else
    {
        cout << NQ << ' ';
        printf("%d %.1f", Q.front().first, Q.front().second);
        auto i = Q.begin();
        i++;
        for (; i != Q.end(); i++)
        {
            printf(" %d %.1f", (*i).first, (*i).second);
        }
        puts("");
    }

    R.clear();
    for (auto i = A.begin(); i != A.end(); i++)
    {
        if (fabs((*i).second) > ESP)
        {
            R.push_back(make_pair((*i).first, (*i).second));
            NR++;
        }
    }

    if (NR == 0)
    {
        cout << "0 0 0.0\n";
    }
    else
    {
        cout << NR << ' ';
        printf("%d %.1f", R.front().first, R.front().second);
        auto i = R.begin();
        i++;
        for (; i != R.end(); i++)
        {
            printf(" %d %.1f", (*i).first, (*i).second);
        }
        puts("");
    }

    return 0;
}

参考

《多项式除法》在线视频