#include using namespace std; namespace arithmetic { template class Addition { public: template T operator+(const V& v) const { T res(static_cast(*this)); return res += static_cast(v); } }; template class Subtraction { public: template T operator-(const V& v) const { T res(static_cast(*this)); return res -= static_cast(v); } }; template class Multiplication { public: template T operator*(const V& v) const { T res(static_cast(*this)); return res *= static_cast(v); } }; template class Division { public: template T operator/(const V& v) const { T res(static_cast(*this)); return res /= static_cast(v); } }; template class Modulus { public: template T operator%(const V& v) const { T res(static_cast(*this)); return res %= static_cast(v); } }; } template class IndivisibleArithmetic : public arithmetic::Addition, public arithmetic::Subtraction, public arithmetic::Multiplication {}; template class Arithmetic : public IndivisibleArithmetic, public arithmetic::Division {}; template class Polynomial : public Arithmetic>, public arithmetic::Modulus> { private: vector val; void normalize() { while (val.size() > 1u && val.back() == 0) val.pop_back(); if (val.empty()) val.emplace_back(0); } public: Polynomial() { normalize(); } Polynomial(const vector& val) : val(val) { normalize(); } Polynomial operator+=(const Polynomial& p) { for (int i = 0; i < p.size(); ++i) { if (int(val.size()) == i) val.emplace_back(p[i]); else val[i] += p[i]; } normalize(); return *this; } Polynomial operator-=(const Polynomial& p) { for (int i = 0; i < p.size(); ++i) { if (int(val.size()) == i) val.emplace_back(-p[i]); else val[i] -= p[i]; } normalize(); return *this; } // TODO FFT Polynomial operator*=(const Polynomial& p) { Polynomial res; for (int i = 0; i < size(); ++i) { for (int j = 0; j < p.size(); ++j) { res[i + j] += val[i] * p[j]; } } *this = res; normalize(); return *this; } Polynomial operator/=(const Polynomial& p) { Polynomial res; for (int i = size() - p.size(); i >= 0; --i) { res[i] = val[p.size() + i - 1] / p.back(); for (int j = 0; j < p.size(); ++j) val[i + j] -= res[i] * p[j]; } *this = res; normalize(); return *this; } Polynomial operator%=(const Polynomial& p) { for (int i = size() - p.size(); i >= 0; --i) { T d = val[p.size() + i - 1] / p.back(); for (int j = 0; j < p.size(); ++j) val[i + j] -= d * p[j]; } normalize(); return *this; } T& operator[](int i) { if (i >= int(val.size())) val.resize(i + 1, 0); return val[i]; } const T& operator[](int i) const { return val[i]; } int size() const { return val.size(); } T& back() { return val.back(); } const T& back() const { return val.back(); } typename vector::iterator begin() { return val.begin(); } typename vector::iterator end() { return val.end(); } Polynomial identity() const { return Polynomial(1, 1); } }; template string to_string(const T& v) { string str; for (const auto& i : const_cast(v)) str += to_string(i) + " "; return str.substr(0, max(0, (int)str.size() - 1)); } int main() { int d; cin >> d; Polynomial poly; for (int i = 0; i <= d; ++i) cin >> poly[i]; poly %= Polynomial(vector({0, -1, 0, 1})); cout << poly.size() - 1 << endl; cout << to_string(poly) << endl; }