Language: C++

#ifndef __INTMOD_H__
#define __INTMOD_H__

/* updated: 2020-06-24 */

#include <iostream>
#include <cassert>
#include <vector>
#include <cstdint>
#include <functional>
#include <unordered_map>

template <uint32_t Mod>
class IntMod {
	static_assert(Mod != 0, "Mod must not be 0");
	static_assert(Mod < 0x80000000, "Mod must be less than 0x80000000");

private:
	uint32_t value_m;

public:
	IntMod();
	IntMod(int32_t value);
	IntMod(uint32_t value);
	IntMod(int64_t value);
	IntMod(uint64_t value);

	const IntMod<Mod>& operator+() const;
	IntMod<Mod> operator-() const;
	IntMod<Mod>& operator++();
	IntMod<Mod>& operator--();
	IntMod<Mod> operator++(int dummy);
	IntMod<Mod> operator--(int dummy);
	IntMod<Mod>& operator+=(const IntMod<Mod>& right);
	IntMod<Mod>& operator-=(const IntMod<Mod>& right);
	IntMod<Mod>& operator*=(const IntMod<Mod>& right);
	IntMod<Mod>& operator/=(const IntMod<Mod>& right);
	IntMod<Mod> operator[](uint64_t exp) const;

	static IntMod<Mod> add(const IntMod<Mod>& a, const IntMod<Mod>& b);
	static IntMod<Mod> subtract(const IntMod<Mod>& a, const IntMod<Mod>& b);
	static IntMod<Mod> multiply(const IntMod<Mod>& a, const IntMod<Mod>& b);
	static IntMod<Mod> divide(const IntMod<Mod>& a, const IntMod<Mod>& b);
	static int compare(const IntMod<Mod>& a, const IntMod<Mod>& b);

	IntMod<Mod> inverse() const;
	IntMod<Mod> power(uint64_t exp) const;
	uint32_t get() const;

	static IntMod<Mod> fact(uint32_t num);
	static IntMod<Mod> invFact(uint32_t num);
	static IntMod<Mod> combi(uint32_t n, uint32_t r);
	static IntMod<Mod> homcombi(uint32_t n, uint32_t r);
	static IntMod<Mod> perm(uint32_t n, uint32_t r);
	static std::vector<IntMod<Mod>> calcInverseList(int max);
};

template <uint32_t Mod>
IntMod<Mod> operator+(const IntMod<Mod>& a, const IntMod<Mod>& b);
template <uint32_t Mod>
IntMod<Mod> operator-(const IntMod<Mod>& a, const IntMod<Mod>& b);
template <uint32_t Mod>
IntMod<Mod> operator*(const IntMod<Mod>& a, const IntMod<Mod>& b);
template <uint32_t Mod>
IntMod<Mod> operator/(const IntMod<Mod>& a, const IntMod<Mod>& b);

template <uint32_t Mod, class Int>
IntMod<Mod> operator+(const IntMod<Mod>& a, const Int& b);
template <uint32_t Mod, class Int>
IntMod<Mod> operator-(const IntMod<Mod>& a, const Int& b);
template <uint32_t Mod, class Int>
IntMod<Mod> operator*(const IntMod<Mod>& a, const Int& b);
template <uint32_t Mod, class Int>
IntMod<Mod> operator/(const IntMod<Mod>& a, const Int& b);

template <uint32_t Mod, class Int>
IntMod<Mod> operator+(const Int& a, const IntMod<Mod>& b);
template <uint32_t Mod, class Int>
IntMod<Mod> operator-(const Int& a, const IntMod<Mod>& b);
template <uint32_t Mod, class Int>
IntMod<Mod> operator*(const Int& a, const IntMod<Mod>& b);
template <uint32_t Mod, class Int>
IntMod<Mod> operator/(const Int& a, const IntMod<Mod>& b);

template <uint32_t Mod>
bool operator==(const IntMod<Mod>& a, const IntMod<Mod>& b);
template <uint32_t Mod>
bool operator!=(const IntMod<Mod>& a, const IntMod<Mod>& b);
template <uint32_t Mod>
bool operator<(const IntMod<Mod>& a, const IntMod<Mod>& b);
template <uint32_t Mod>
bool operator<=(const IntMod<Mod>& a, const IntMod<Mod>& b);
template <uint32_t Mod>
bool operator>(const IntMod<Mod>& a, const IntMod<Mod>& b);
template <uint32_t Mod>
bool operator>=(const IntMod<Mod>& a, const IntMod<Mod>& b);

template <uint32_t Mod, class Int>
bool operator==(const IntMod<Mod>& a, const Int& b);
template <uint32_t Mod, class Int>
bool operator!=(const IntMod<Mod>& a, const Int& b);

template <uint32_t Mod>
std::istream& operator>>(std::istream& ist, IntMod<Mod>& object);
template <uint32_t Mod>
std::ostream& operator<<(std::ostream& ost, const IntMod<Mod>& object);

using MInt = IntMod<1000000007>;
using MF = IntMod<998244353>;

// for unordered_set/map
namespace std {
	template <uint32_t Mod>
	struct hash<IntMod<Mod>> {
		size_t operator() (const IntMod<Mod>& object) const {
			return std::hash<uint32_t>()(object.get());
		}
	};
}

/*------------------------------------*/
/*-------   Implementations   --------*/
/*------------------------------------*/

template <uint32_t Mod>
IntMod<Mod>::IntMod() 
	: value_m(0) {
}

template <uint32_t Mod>
IntMod<Mod>::IntMod(int32_t value) {
	int32_t tmp = value % int32_t(Mod);
	value_m = tmp >= 0 ? tmp : Mod - uint32_t(-tmp);
}

template <uint32_t Mod>
IntMod<Mod>::IntMod(uint32_t value)
	: value_m(value % Mod) {
}

template <uint32_t Mod>
IntMod<Mod>::IntMod(int64_t value) {
	int32_t tmp = value % int64_t(Mod);
	value_m = tmp >= 0 ? tmp : Mod - uint32_t(-tmp);
}

template <uint32_t Mod>
IntMod<Mod>::IntMod(uint64_t value)
	: value_m(value % Mod) {
}

template <uint32_t Mod>
const IntMod<Mod>& IntMod<Mod>::operator+() const {
	return *this;
}

template <uint32_t Mod>
IntMod<Mod> IntMod<Mod>::operator-() const {
	IntMod<Mod> ret;
	if (value_m == 0) {
		ret.value_m = 0;
	} else {
		ret.value_m = Mod - value_m;
	}
	return ret;
}

template <uint32_t Mod>
IntMod<Mod>& IntMod<Mod>::operator++() {
	++value_m;
	if (value_m == Mod) {
		value_m = 0;
	}
	return *this;
}

template <uint32_t Mod>
IntMod<Mod>& IntMod<Mod>::operator--() {
	if (value_m == 0) {
		value_m = Mod;
	}
	--value_m;
	return *this;
}

template <uint32_t Mod>
IntMod<Mod> IntMod<Mod>::operator++(int dummy) {
	IntMod ret(*this);
	++*this;
	return ret;
}

template <uint32_t Mod>
IntMod<Mod> IntMod<Mod>::operator--(int dummy) {
	IntMod ret(*this);
	--*this;
	return ret;
}

template <uint32_t Mod>
IntMod<Mod>& IntMod<Mod>::operator+=(const IntMod<Mod>& right) {
	value_m += right.value_m;		// note: value_m < 0x80000000
	if (value_m >= Mod) {
		value_m -= Mod;
	}
	return *this;
}

template <uint32_t Mod>
IntMod<Mod>& IntMod<Mod>::operator-=(const IntMod<Mod>& right) {
	if (value_m < right.value_m) {
		value_m += Mod;				// note: value_m < 0x80000000
	}
	value_m -= right.value_m;
	return *this;
}

template <uint32_t Mod>
IntMod<Mod>& IntMod<Mod>::operator*=(const IntMod<Mod>& right) {
	value_m = (uint64_t(value_m) * right.value_m) % Mod;
	return *this;
}

template <uint32_t Mod>
IntMod<Mod>& IntMod<Mod>::operator/=(const IntMod<Mod>& right) {
	*this *= right.inverse();
	return *this;
}

// for power
template <uint32_t Mod>
IntMod<Mod> IntMod<Mod>::operator[](uint64_t exp) const {
	return power(exp);
}

template <uint32_t Mod>
IntMod<Mod> IntMod<Mod>::add(const IntMod<Mod>& a, const IntMod<Mod>& b) {
	IntMod<Mod> ret(a);
	ret += b;
	return ret;
}

template <uint32_t Mod>
IntMod<Mod> IntMod<Mod>::subtract(const IntMod<Mod>& a, const IntMod<Mod>& b) {
	IntMod<Mod> ret(a);
	ret -= b;
	return ret;
}

template <uint32_t Mod>
IntMod<Mod> IntMod<Mod>::multiply(const IntMod<Mod>& a, const IntMod<Mod>& b) {
	IntMod<Mod> ret(a);
	ret *= b;
	return ret;
}

template <uint32_t Mod>
IntMod<Mod> IntMod<Mod>::divide(const IntMod<Mod>& a, const IntMod<Mod>& b) {
	IntMod<Mod> ret(a);
	ret /= b;
	return ret;
}

template <uint32_t Mod>
int IntMod<Mod>::compare(const IntMod<Mod>& a, const IntMod<Mod>& b) {
	if (a.value_m < b.value_m) return -1;
	if (a.value_m > b.value_m) return 1;
	return 0;
}

/* Mod must be a prime. */
template <uint32_t Mod>
IntMod<Mod> IntMod<Mod>::inverse() const {
	assert(value_m != 0);
	return power(Mod - 2);
}

template <uint32_t Mod>
IntMod<Mod> IntMod<Mod>::power(uint64_t exp) const {
	IntMod<Mod> product(1);
	IntMod<Mod> factor(*this);
	while (exp > 0) {
		if (exp & 1) {
			product *= factor;
		}
		factor *= factor;
		exp >>= 1;
	}
	return product;
}

template <uint32_t Mod>
uint32_t IntMod<Mod>::get() const {
	return value_m;
}

template <uint32_t Mod>
IntMod<Mod> IntMod<Mod>::fact(uint32_t num) {
	assert(num <= 0x80000000);

	static std::vector<IntMod<Mod>> table(1, 1);
	if (table.size() > num) return table[num];

	uint32_t old_size = table.size();
	table.resize(num + 1);
	for (uint32_t i = old_size; i <= num; ++i) {
		table[i] = table[i - 1] * i;
	}
	return table[num];
}

template <uint32_t Mod>
IntMod<Mod> IntMod<Mod>::invFact(uint32_t num) {
	assert(num <= 0x80000000);

	static std::vector<IntMod<Mod>> table(1, 1);
	if (table.size() > num) return table[num];

	uint32_t old_size = table.size();
	table.resize(num + 1);
	table[num] = fact(num).inverse();
	for (uint32_t i = num; i >= old_size; --i) {
		table[i - 1] = table[i] * i;
	}
	return table[num];
}

template <uint32_t Mod>
IntMod<Mod> IntMod<Mod>::combi(uint32_t n, uint32_t r) {
	assert(n >= r);
	return fact(n) * invFact(n - r) * invFact(r);
}

template <uint32_t Mod>
IntMod<Mod> IntMod<Mod>::homcombi(uint32_t n, uint32_t r) {
	assert(n != 0 || r == 0);

	if (n == 0 && r == 0) return IntMod<Mod>(1);
	return combi(n + r - 1, r);
}

template <uint32_t Mod>
IntMod<Mod> IntMod<Mod>::perm(uint32_t n, uint32_t r) {
	assert(n >= r);
	return fact(n) * invFact(n - r);
}

template <uint32_t Mod>
std::vector<IntMod<Mod>> IntMod<Mod>::calcInverseList(int max) {
	assert(max < Mod);

	std::vector<IntMod<Mod>> table(max + 1);
	table[1] = 1;
	for (int i = 2; i <= max; ++i) {
		table[i] = (uint64_t(Mod - Mod / i) * table[Mod % i].get()) % Mod;
	}
	return table;
}

template <uint32_t Mod>
IntMod<Mod> operator+(const IntMod<Mod>& a, const IntMod<Mod>& b) {
	return IntMod<Mod>::add(a, b);
}

template <uint32_t Mod>
IntMod<Mod> operator-(const IntMod<Mod>& a, const IntMod<Mod>& b) {
	return IntMod<Mod>::subtract(a, b);
}

template <uint32_t Mod>
IntMod<Mod> operator*(const IntMod<Mod>& a, const IntMod<Mod>& b) {
	return IntMod<Mod>::multiply(a, b);
}

template <uint32_t Mod>
IntMod<Mod> operator/(const IntMod<Mod>& a, const IntMod<Mod>& b) {
	return IntMod<Mod>::divide(a, b);
}

template <uint32_t Mod, class Int>
IntMod<Mod> operator+(const IntMod<Mod>& a, const Int& b) {
	return a + IntMod<Mod>(b);
}

template <uint32_t Mod, class Int>
IntMod<Mod> operator-(const IntMod<Mod>& a, const Int& b) {
	return a - IntMod<Mod>(b);
}

template <uint32_t Mod, class Int>
IntMod<Mod> operator*(const IntMod<Mod>& a, const Int& b) {
	return a * IntMod<Mod>(b);
}

template <uint32_t Mod, class Int>
IntMod<Mod> operator/(const IntMod<Mod>& a, const Int& b) {
	return a / IntMod<Mod>(b);
}

template <uint32_t Mod, class Int>
IntMod<Mod> operator+(const Int& a, const IntMod<Mod>& b) {
	return IntMod<Mod>(a) + b;
}

template <uint32_t Mod, class Int>
IntMod<Mod> operator-(const Int& a, const IntMod<Mod>& b) {
	return IntMod<Mod>(a) - b;
}

template <uint32_t Mod, class Int>
IntMod<Mod> operator*(const Int& a, const IntMod<Mod>& b) {
	return IntMod<Mod>(a) * b;
}

template <uint32_t Mod, class Int>
IntMod<Mod> operator/(const Int& a, const IntMod<Mod>& b) {
	return IntMod<Mod>(a) / b;
}

template <uint32_t Mod>
bool operator==(const IntMod<Mod>& a, const IntMod<Mod>& b) {
	return IntMod<Mod>::compare(a, b) == 0;
}

template <uint32_t Mod>
bool operator!=(const IntMod<Mod>& a, const IntMod<Mod>& b) {
	return IntMod<Mod>::compare(a, b) != 0;
}

template <uint32_t Mod>
bool operator<(const IntMod<Mod>& a, const IntMod<Mod>& b) {
	return IntMod<Mod>::compare(a, b) < 0;
}

template <uint32_t Mod>
bool operator<=(const IntMod<Mod>& a, const IntMod<Mod>& b) {
	return IntMod<Mod>::compare(a, b) <= 0;
}

template <uint32_t Mod>
bool operator>(const IntMod<Mod>& a, const IntMod<Mod>& b) {
	return IntMod<Mod>::compare(a, b) > 0;
}

template <uint32_t Mod>
bool operator>=(const IntMod<Mod>& a, const IntMod<Mod>& b) {
	return IntMod<Mod>::compare(a, b) >= 0;
}

template <uint32_t Mod, class Int>
bool operator==(const IntMod<Mod>& a, const Int& b) {
	return a == IntMod<Mod>(b);
}

template <uint32_t Mod, class Int>
bool operator!=(const IntMod<Mod>& a, const Int& b) {
	return a != IntMod<Mod>(b);
}

template <uint32_t Mod>
std::istream& operator>>(std::istream& ist, IntMod<Mod>& object) {
	uint64_t value;
	ist >> value;
	object = IntMod<Mod>(value);
	return ist;
}

template <uint32_t Mod>
std::ostream& operator<<(std::ostream& ost, const IntMod<Mod>& object) {
	ost << object.get();
	return ost;
}

#endif
#ifndef __RUNTIME_MOD_H__
#define __RUNTIME_MOD_H__

/* updated: 2020-06-24 */

// required: IntMod

#ifdef ONLY_MY_ENVIR
#include "IntMod.h"
#endif

#include <iostream>
#include <cassert>
#include <vector>
#include <cstdint>
#include <functional>
#include <unordered_map>

/* 実行時に mod が与えられるときのための IntMod, 内容は IntMod と同じ */
/* mod の変更は禁止: fact() などが動作しないため */
template <>
class IntMod<0> {
	uint32_t value_m;

	static uint32_t mod_m;
public:
	IntMod();
	IntMod(int32_t value);
	IntMod(uint32_t value);
	IntMod(int64_t value);
	IntMod(uint64_t value);

	const IntMod<0>& operator+() const;
	IntMod<0> operator-() const;
	IntMod<0>& operator++();
	IntMod<0>& operator--();
	IntMod<0> operator++(int dummy);
	IntMod<0> operator--(int dummy);
	IntMod<0>& operator+=(const IntMod<0>& right);
	IntMod<0>& operator-=(const IntMod<0>& right);
	IntMod<0>& operator*=(const IntMod<0>& right);
	IntMod<0>& operator/=(const IntMod<0>& right);
	IntMod<0> operator[](uint64_t exp) const;

	static IntMod<0> add(const IntMod<0>& a, const IntMod<0>& b);
	static IntMod<0> subtract(const IntMod<0>& a, const IntMod<0>& b);
	static IntMod<0> multiply(const IntMod<0>& a, const IntMod<0>& b);
	static IntMod<0> divide(const IntMod<0>& a, const IntMod<0>& b);
	static int compare(const IntMod<0>& a, const IntMod<0>& b);

	static void setDefaultMod(uint32_t mod);

	IntMod<0> inverse() const;
	IntMod<0> power(uint64_t exp) const;
	uint32_t get() const;

	static IntMod<0> fact(uint32_t num);
	static IntMod<0> invFact(uint32_t num);
	static IntMod<0> combi(uint32_t n, uint32_t r);
	static IntMod<0> homcombi(uint32_t n, uint32_t r);
	static IntMod<0> perm(uint32_t n, uint32_t r);
	static std::vector<IntMod<0>> calcInverseList(int max);
};

// グローバルの演算子はテンプレート版があれば動作するので特殊化をしない

using RuntimeMod = IntMod<0>;

#endif

#ifndef __INTERPOLATION_H__
#define __INTERPOLATION_H__

#include <vector>

template <class T>
std::vector<T> interpolatingPolynomial(const std::vector<T>& points, const std::vector<T>& values);
template <class T>
T interpolate(const std::vector<T>& values, T target);


/* 体 T 上の多項式補間 O(n^2) */
/* データ点 (points[i], values[i]) を全て通る多項式を返す (vector の i 番目の要素が x^i の係数) */
template <class T>
std::vector<T> interpolatingPolynomial(const std::vector<T>& points, const std::vector<T>& values) {
	assert(points.size() == values.size());
	int size = points.size();

	/* 多項式 p と多項式 (x + a) の積を返す */
	auto multiply = [](const std::vector<T>& p, T a) -> std::vector<T> {
		int length = p.size();
		std::vector<T> ret(length + 1);
		for (int i = 0; i < length; ++i) {
			ret[i] += a * p[i];
			ret[i + 1] += p[i];
		}
		return ret;
	};

	/* 多項式 p を多項式 (x + a) で割った商を返す */
	auto divide = [](const std::vector<T>& p, T a) -> std::vector<T> {
		int length = p.size();
		assert(length > 0);

		std::vector<T> ret(length - 1);
		T tmp = p[length - 1];
		for (int i = length - 1; i > 0; --i) {
			ret[i - 1] = tmp;
			tmp = p[i - 1] - tmp * a;
		}
		return ret;
	};
	
	/* product [0 <= i < n] (x - points[i]) を求める */
	std::vector<T> product(1, T(1));
	for (int i = 0; i < size; ++i) {
		product = multiply(product, -points[i]);
	}

	/* 答えは (values[i] / denom) * res の和 */
	std::vector<T> sum(size);
	for (int i = 0; i < size; ++i) {
		/* denom = product [k != i] (points[i] - points[k]) */
		T denom = 1;
		for (int k = 0; k < size; ++k) {
			if (k == i) continue;
			denom *= points[i] - points[k];
		}

		T coeff = values[i] / denom;

		/* res = product [k != i] (x - points[k]) */
		std::vector<T> res = divide(product, -points[i]);
		for (int j = 0; j < size; ++j) {
			sum[j] += coeff * res[j];
		}
	}
	return sum;
}


/* 体 T 上の多項式補間 O(n) or O(n log mod) */
/* x = 0, 1, ..., n - 1 に対する値 p(x) が既知のときに p(target) を求める */
template <class T>
T interpolate(const std::vector<T>& values, T target) {
	const int size = values.size();
	std::vector<T> left(size);		// left[i]  = (target - k) の [0, i) に関する積
	std::vector<T> right(size);		// right[i] = (target - k) の (i, n - 1] に関する積

	left[0] = 1;
	right[size - 1] = 1;
	for (int i = 0; i < size - 1; ++i) {
		left[i + 1] = left[i] * (target - i);
	}

	for (int i = size - 1; i > 0; --i) {
		right[i - 1] = right[i] * (target - i);
	}

	std::vector<T> fact(size);		// fact[i] = i の階乗 (添字に注意)
	fact[0] = 1;
	for (int i = 1; i < size; ++i) {
		fact[i] = fact[i - 1] * i;
	}

	// 答えは sum[0 <= i < n] { values[i] * product[k != i] { (target - k) / (i - k) } }
	T sum = T();
	for (int i = 0; i < size; ++i) {
		T sign = (size - 1 - i) % 2 == 0 ? 1 : -1;
		T numer = left[i] * right[i];
		T denom = fact[i] * fact[size - 1 - i];
		sum += sign * values[i] * numer / denom;
	}
	return sum;
}


#endif

#define _CRT_SECURE_NO_WARNINGS
#define _SCL_SECURE_NO_WARNINGS
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <cassert>
#include <iostream>
#include <string>
#include <vector>
#include <list>
#include <utility>
#include <algorithm>
#include <functional>
#include <cmath>
#include <stack>
#include <queue>
#include <set>
#include <map>
#include <iomanip>
#include <sstream>
#include <bitset>
#include <limits>
#include <numeric>
#include <valarray>
#include <fstream>
#include <array>
#include <random>
#include <unordered_set>
#include <unordered_map>

using namespace std;
using uint = uint32_t;
using LL = int64_t;
using ULL = uint64_t;
using PP = pair<LL, LL>;
template <typename T> using PriorityQ = priority_queue<T, vector<T>, greater<T> >;
#define REP(i, a, n) for(LL i = (a), i##_max_ = (n); i < i##_max_; ++i)
#define REM(i, a, n) for(LL i = (LL)(n) - 1, i##_min_ = (a); i >= i##_min_; --i)
#define FLOAT fixed << setprecision(16)
#define SPEEDUP { cin.tie(NULL); ios::sync_with_stdio(false); }
const int INF = 0x3FFFFFFF;
const LL INFLL = 0x3FFFFFFF3FFFFFFF;
const double INFD = 1.0e+308;
const double EPS = 1.0e-9;

void YesNo(bool b) { cout << (b ? "Yes" : "No") << endl; }
void YESNO(bool b) { cout << (b ? "YES" : "NO") << endl; }
template <class T, class U> istream& operator>>(istream& ist, pair<T, U>& right) { return ist >> right.first >> right.second; }
template <class T, class U> ostream& operator<<(ostream& ost, const pair<T, U>& right) { return ost << right.first << ' ' << right.second; }
template <class T, class TCompatible, size_t N> void Fill(T(&dest)[N], const TCompatible& val) { fill(dest, dest + N, val); }
template <class T, class TCompatible, size_t M, size_t N> void Fill(T(&dest)[M][N], const TCompatible& val) { for (int i = 0; i < M; ++i) Fill(dest[i], val); }
template <class T> T Next() { T buf; cin >> buf; return buf; }
istream& Ignore(istream& ist) { string s; ist >> s; return ist; }
bool Inside(int i, int j, int h, int w) { return i >= 0 && i < h && j >= 0 && j < w; }

#ifdef ONLY_MY_ENVIR
#include "Accumulator.h"
#include "Algebraic.h"
#include "BinaryMatrix.h"
#include "BinaryTree.h"
#include "Bipartite.h"
#include "BIT.h"
#include "Compressor.h"
#include "Decompositions.h"
#include "DiscreteLog.h"
#include "DynamicMod.h"
#include "Exponential.h"
#include "Factorization.h"
#include "FFT.h"
#include "FlowSolver.h"
#include "Graph.h"
#include "GraphUtil.h"
#include "Interpolation.h"
#include "IntMod.h"
#include "LazySegmentTree.h"
#include "LIS.h"
#include "Math.h"
#include "MathUtil.h"
#include "Matrix.h"
#include "MinCostFlowSolver.h"
#include "MinMax.h"
#include "Numbers.h"
#include "Optimize.h"
#include "Permutation.h"
#include "Polynomial.h"
#include "Position.h"
#include "Random.h"
#include "Range.h"
#include "Rational.h"
#include "RuntimeMod.h"
#include "SegmentTree.h"
#include "SegmentTree2D.h"
#include "Sets.h"
#include "Shortest.h"
#include "SlidingWindow.h"
#include "SpanningTree.h"
#include "StringSearching.h"
#include "SuffixArray.h"
#include "SwitchList.h"
#include "Timer.h"
#include "Tree.h"
#include "TreeUtil.h"
#include "UnionFind.h"
#include "Util.h"
#include "VectorUtil.h"
#endif

#ifdef __GNUC__
typedef __int128 LLL;
istream& operator>>(istream& ist, __int128& val) { LL tmp; ist >> tmp; val = tmp; return ist; }
ostream& operator<<(ostream& ost, __int128 val) { LL tmp = val; ost << tmp; return ost; }
#endif

int M;
vector<RuntimeMod> a;
vector<RuntimeMod> b;
int main() {
	cin >> M;
	RuntimeMod::setDefaultMod(M);
	REP(i, 0, M) a.push_back(i);
	REP(i, 0, M) b.push_back(Next<int>());
	auto res = interpolatingPolynomial(a, b);
	REP(i, 0, M) cout << res[i] << ' ';
	cout << endl;
	return 0;
}



Compiler output:
/usr/bin/ld: /tmp/ccVO45M2.o: in function `main':
a.cc:(.text+0x1da): undefined reference to `IntMod<0u>::setDefaultMod(unsigned int)'
/usr/bin/ld: a.cc:(.text+0x203): undefined reference to `IntMod<0u>::IntMod(long)'
/usr/bin/ld: a.cc:(.text+0x253): undefined reference to `IntMod<0u>::IntMod(int)'
/usr/bin/ld: /tmp/ccVO45M2.o: in function `interpolatingPolynomial<IntMod<0u> >(std::vector<IntMod<0u>, std::allocator<IntMod<0u> > > const&, std::vector<IntMod<0u>, std::allocator<IntMod<0u> > > const&)::{lambda(std::vector<IntMod<0u>, std::allocator<IntMod<0u> > > const&, IntMod<0u>)#1}::operator()(std::vector<IntMod<0u>, std::allocator<IntMod<0u> > > const&, IntMod<0u>) const':
a.cc:(.text._ZZ23interpolatingPolynomialI6IntModILj0EEESt6vectorIT_SaIS3_EERKS5_S7_ENKUlRKS2_IS1_SaIS1_EES1_E_clESB_S1_[_ZZ23interpolatingPolynomialI6IntModILj0EEESt6vectorIT_SaIS3_EERKS5_S7_ENKUlRKS2_IS1_SaIS1_EES1_E_clESB_S1_]+0xb3): undefined reference to `IntMod<0u>::operator+=(IntMod<0u> const&)'
/usr/bin/ld: a.cc:(.text._ZZ23interpolatingPolynomialI6IntModILj0EEESt6vectorIT_SaIS3_EERKS5_S7_ENKUlRKS2_IS1_SaIS1_EES1_E_clESB_S1_[_ZZ23interpolatingPolynomialI6IntModILj0EEESt6vectorIT_SaIS3_EERKS5_S7_ENKUlRKS2_IS1_SaIS1_EES1_E_clESB_S1_]+0xee): undefined reference to `IntMod<0u>::operator+=(IntMod<0u> const&)'
/usr/bin/ld: /tmp/ccVO45M2.o: in function `std::vector<IntMod<0u>, std::allocator<IntMod<0u> > > interpolatingPolynomial<IntMod<0u> >(std::vector<IntMod<0u>, std::allocator<IntMod<0u> > > const&, std::vector<IntMod<0u>, std::allocator<IntMod<0u> > > const&)':
a.cc:(.text._Z23interpolatingPolynomialI6IntModILj0EEESt6vectorIT_SaIS3_EERKS5_S7_[_Z23interpolatingPolynomialI6IntModILj0EEESt6vectorIT_SaIS3_EERKS5_S7_]+0x98): undefined reference to `IntMod<0u>::IntMod(int)'
/usr/bin/ld: a.cc:(.text._Z23interpolatingPolynomialI6IntModILj0EEESt6vectorIT_SaIS3_EERKS5_S7_[_Z23interpolatingPolynomialI6IntModILj0EEESt6vectorIT_SaIS3_EERKS5_S7_]+0xe7): undefined reference to `IntMod<0u>::operator-() const'
/usr/bin/ld: a.cc:(.text._Z23interpolatingPolynomialI6IntModILj0EEESt6vectorIT_SaIS3_EERKS5_S7_[_Z23interpolatingPolynomialI6IntModILj0EEESt6vectorIT_SaIS3_EERKS5_S7_]+0x17b): undefined reference to `IntMod<0u>::IntMod(int)'
/usr/bin/ld: a.cc:(.text._Z23interpolatingPolynomialI6IntModILj0EEESt6vectorIT_SaIS3_EERKS5_S7_[_Z23interpolatingPolynomialI6IntModILj0EEESt6vectorIT_SaIS3_EERKS5_S7_]+0x1e3): undefined reference to `IntMod<0u>::operator*=(IntMod<0u> const&)'
/usr/bin/ld: a.cc:(.text._Z23interpolatingPolynomialI6IntModILj0EEESt6vectorIT_SaIS3_EERKS5_S7_[_Z23interpolatingPolynomialI6IntModILj0EEESt6vectorIT_SaIS3_EERKS5_S7_]+0x245): undefined reference to `IntMod<0u>::operator-() const'
/usr/bin/ld: a.cc:(.text._Z23interpolatingPolynomialI6IntModILj0EEESt6vectorIT_SaIS3_EERKS5_S7_[_Z23interpolatingPolynomialI6IntModILj0EEESt6vectorIT_SaIS3_EERKS5_S7_]+0x2c1): undefined reference to `IntMod<0u>::operator+=(IntMod<0u> const&)'
/usr/bin/ld: /tmp/ccVO45M2.o: in function `std::ostream& operator<< <0u>(std::ostream&, IntMod<0u> const&)':
a.cc:(.text._ZlsILj0EERSoS0_RK6IntModIXT_EE[_ZlsILj0EERSoS0_RK6IntModIXT_EE]+0x18): undefined reference to `IntMod<0u>::get() const'
/usr/bin/ld: /tmp/ccVO45M2.o: in function `IntMod<0u> operator*<0u>(IntMod<0u> const&, IntMod<0u> const&)':
a.cc:(.text._ZmlILj0EE6IntModIXT_EERKS1_S3_[_ZmlILj0EE6IntModIXT_EERKS1_S3_]+0x1f): undefined reference to `IntMod<0u>::multiply(IntMod<0u> const&, IntMod<0u> const&)'
/usr/bin/ld: /tmp/ccVO45M2.o: in function `IntMod<0u> operator-<0u>(IntMod<0u> const&, IntMod<0u> const&)':
a.cc:(.text._ZmiILj0EE6IntModIXT_EERKS1_S3_[_ZmiILj0EE6IntModIXT_EERKS1_S3_]+0x1f): undefined reference to `IntMod<0u>::subtract(IntMod<0u> const&, IntMod<0u> const&)'
/usr/bin/ld: /tmp/ccVO45M2.o: in function `IntMod<0u> operator/<0u>(IntMod<0u> const&, IntMod<0u> const&)':
a.cc:(.text._ZdvILj0EE6IntModIXT_EERKS1_S3_[_ZdvILj0EE6IntModIXT_EERKS1_S3_]+0x1f): undefined reference to `IntMod<0u>::divide(IntMod<0u> const&, IntMod<0u> const&)'
/usr/bin/ld: /tmp/ccVO45M2.o: in function `void std::_Construct<IntMod<0u>>(IntMod<0u>*)':
a.cc:(.text._ZSt10_ConstructI6IntModILj0EEJEEvPT_DpOT0_[_ZSt10_ConstructI6IntModILj0EEJEEvPT_DpOT0_]+0x31): undefined reference to `IntMod<0u>::IntMod()'
collect2: error: ld returned 1 exit status