Submission #3777000
Source Code Expand
#define _CRT_SECURE_NO_WARNINGS
#include "bits/stdc++.h"
//#include "libs.h"
#include <random>
#include <unordered_map>
#include <unordered_set>
//#include <opencv2/core.hpp>
//#include <opencv2/highgui.hpp>
//#include <opencv2/imgproc.hpp>
using namespace std;
//呪文
#define DUMPOUT cerr
#define dump(...) DUMPOUT<<" ";DUMPOUT<<#__VA_ARGS__<<" :["<<__LINE__<<":"<<__FUNCTION__<<"]"<<endl;DUMPOUT<<" ";dump_func(__VA_ARGS__)
typedef unsigned uint; typedef long long ll; typedef unsigned long long ull; typedef pair<int, int> pii; typedef pair<ll, ll> pll; typedef pair<double, double> pdd; typedef pair<string, string> pss;
template <typename _KTy, typename _Ty> ostream& operator << (ostream& o, const pair<_KTy, _Ty>& m) { o << "{" << m.first << ", " << m.second << "}"; return o; }
template <typename _KTy, typename _Ty> ostream& operator << (ostream& o, const map<_KTy, _Ty>& m) { if (m.empty()) { o << "{ }"; return o; } o << "{" << *m.begin(); for (auto itr = ++m.begin(); itr != m.end(); itr++) { o << ", " << *itr; } o << "}"; return o; }
template <typename _KTy, typename _Ty> ostream& operator << (ostream& o, const unordered_map<_KTy, _Ty>& m) { if (m.empty()) { o << "{ }"; return o; } o << "{" << *m.begin(); for (auto itr = ++m.begin(); itr != m.end(); itr++) { o << ", " << *itr; } o << "}"; return o; }
template <typename _Ty> ostream& operator << (ostream& o, const vector<_Ty>& v) { if (v.empty()) { o << "{ }"; return o; } o << "{" << v.front(); for (auto itr = ++v.begin(); itr != v.end(); itr++) { o << ", " << *itr; } o << "}"; return o; }
template <typename _Ty> ostream& operator << (ostream& o, const set<_Ty>& s) { if (s.empty()) { o << "{ }"; return o; } o << "{" << *(s.begin()); for (auto itr = ++s.begin(); itr != s.end(); itr++) { o << ", " << *itr; } o << "}"; return o; }
template <typename _Ty> ostream& operator << (ostream& o, const unordered_set<_Ty>& s) { if (s.empty()) { o << "{ }"; return o; } o << "{" << *(s.begin()); for (auto itr = ++s.begin(); itr != s.end(); itr++) { o << ", " << *itr; } o << "}"; return o; }
template <typename _Ty> ostream& operator << (ostream& o, const stack<_Ty>& s) { if (s.empty()) { o << "{ }"; return o; } stack<_Ty> t(s); o << "{" << t.top(); t.pop(); while (!t.empty()) { o << ", " << t.top(); t.pop(); } o << "}"; return o; }
template <typename _Ty> ostream& operator << (ostream& o, const list<_Ty>& l) { if (l.empty()) { o << "{ }"; return o; } o << "{" << l.front(); for (auto itr = ++l.begin(); itr != l.end(); ++itr) { o << ", " << *itr; } o << "}"; return o; }
template <typename _KTy, typename _Ty> istream& operator >> (istream& is, pair<_KTy, _Ty>& m) { is >> m.first >> m.second; return is; }
template <typename _Ty> istream& operator >> (istream& is, vector<_Ty>& v) { for (size_t i = 0; i < v.size(); i++) is >> v[i]; return is; }
namespace aux { // print tuple
template<typename Ty, unsigned N, unsigned L> struct tp { static void print(ostream& os, const Ty& v) { os << get<N>(v) << ", "; tp<Ty, N + 1, L>::print(os, v); } };
template<typename Ty, unsigned N> struct tp<Ty, N, N> { static void print(ostream& os, const Ty& v) { os << get<N>(v); } };
}
template<typename... Tys> ostream& operator<<(ostream& os, const tuple<Tys...>& t) { os << "{"; aux::tp<tuple<Tys...>, 0, sizeof...(Tys) - 1>::print(os, t); os << "}"; return os; }
template<typename A, size_t N, typename T> inline void Fill(A(&array)[N], const T &val) { std::fill((T*)array, (T*)(array + N), val); }
void dump_func() { DUMPOUT << endl; }
template <class Head, class... Tail> void dump_func(Head&& head, Tail&&... tail) { DUMPOUT << head; if (sizeof...(Tail) == 0) { DUMPOUT << " "; } else { DUMPOUT << ", "; } dump_func(std::move(tail)...); }
#define PI 3.14159265358979323846
#define EPS 1e-10
#define FOR(i,a,n) for(int i=(a);i<(n);++i)
#define REP(i,n) FOR(i,0,n)
#define all(j) (j).begin(), (j).end()
#define SZ(j) ((int)(j).size())
#define fake false
constexpr ll MOD = 1000000007;
ll modfac[1010];
ll modfacinv[1010];
// 拡張ユークリッド互除法 (入力 : 整数 a, n, 出力 : GCD(a, n) と xa + yb = GCD(a, n) を満たす整数 x, y)
template<typename Ty> void ExEuclid(Ty a, Ty b, Ty& gcd, Ty& x, Ty& y) {
Ty r1, r2, r3, x1, x2, x3, y1, y2, y3, q;
r1 = a; x1 = 1; y1 = 0;
r2 = b; x2 = 0; y2 = 1;
while (r2 != 0) {
r3 = r1 % r2;
q = r1 / r2;
x3 = x1 - q * x2;
y3 = y1 - q * y2;
r1 = r2; r2 = r3;
x1 = x2; x2 = x3;
y1 = y2; y2 = y3;
}
gcd = r1; x = x1; y = y1;
if (gcd < 0) { gcd *= -1; x *= -1; y *= -1; };
}
// a と m は互いに素が前提
template<typename Ty> Ty invmod(Ty a, Ty m) {
Ty gcd, x, q;
ExEuclid(a, -m, gcd, x, q);
if (gcd != 1) return -1;
return x;
}
// n^e mod m
template<typename Ty> Ty modpow(Ty n, Ty e, Ty m) {
Ty result = 1;
while (e > 0) {
if (e & 1) result = (result * n) % m;
e >>= 1;
n = (n * n) % m;
}
return result;
}
int main() {
cin.tie(0);
ios::sync_with_stdio(false);
modfac[0] = 1;
for (ll i = 1; i <= 1000; i++)
modfac[i] = modfac[i - 1] * i % MOD;
for (ll i = 0; i <= 1000; i++)
modfacinv[i] = invmod(modfac[i], MOD);
ll N, A, B, C, D;
cin >> N >> A >> B >> C >> D;
vector<vector<ll>> dp(B + 1, vector<ll>(N + 1, 0));
REP(i, A) dp[i][0] = 1;
for (ll i = A; i <= B; i++) {
for (ll j = 0; j <= N; j++) {
dp[i][j] += dp[i - 1][j];
for (ll k = C; k <= D; k++) {
ll n = j - k * i;
if (n < 0) break;
ll c = dp[i - 1][n] * modfac[N - n] % MOD;
c = c * /*invmod(modfac[N - j], MOD)*/ modfacinv[N - j] % MOD;
c = c * /*invmod(modfac[k], MOD)*/ modfacinv[k] % MOD;
c = c * invmod(modpow(modfac[i], k, MOD), MOD) % MOD;
dp[i][j] = (dp[i][j] + c) % MOD;
}
}
}
cout << (dp[B][N] + MOD) % MOD << endl;
return 0;
}
Submission Info
Submission Time |
|
Task |
E - Grouping |
User |
komori3 |
Language |
C++14 (GCC 5.4.1) |
Score |
600 |
Code Size |
5866 Byte |
Status |
AC |
Exec Time |
835 ms |
Memory |
8192 KB |
Judge Result
Set Name |
Sample |
All |
Score / Max Score |
0 / 0 |
600 / 600 |
Status |
|
|
Set Name |
Test Cases |
Sample |
sample_01.txt, sample_02.txt, sample_03.txt, sample_04.txt |
All |
sample_01.txt, sample_02.txt, sample_03.txt, sample_04.txt, subtask_1_many_01.txt, subtask_1_many_02.txt, subtask_1_many_03.txt, subtask_1_many_04.txt, subtask_1_max_01.txt, subtask_1_max_02.txt, subtask_1_min_01.txt, subtask_1_randa_01.txt, subtask_1_randa_02.txt, subtask_1_randb_01.txt, subtask_1_randb_02.txt |
Case Name |
Status |
Exec Time |
Memory |
sample_01.txt |
AC |
2 ms |
256 KB |
sample_02.txt |
AC |
2 ms |
256 KB |
sample_03.txt |
AC |
835 ms |
8192 KB |
sample_04.txt |
AC |
2 ms |
256 KB |
subtask_1_many_01.txt |
AC |
580 ms |
896 KB |
subtask_1_many_02.txt |
AC |
485 ms |
7680 KB |
subtask_1_many_03.txt |
AC |
582 ms |
896 KB |
subtask_1_many_04.txt |
AC |
586 ms |
3968 KB |
subtask_1_max_01.txt |
AC |
14 ms |
7552 KB |
subtask_1_max_02.txt |
AC |
9 ms |
3584 KB |
subtask_1_min_01.txt |
AC |
2 ms |
256 KB |
subtask_1_randa_01.txt |
AC |
3 ms |
1920 KB |
subtask_1_randa_02.txt |
AC |
2 ms |
256 KB |
subtask_1_randb_01.txt |
AC |
5 ms |
3328 KB |
subtask_1_randb_02.txt |
AC |
4 ms |
512 KB |