oj.uz

Submission #741181

# Submission time Handle Problem Language Result Execution time Memory
741181 2023-05-13T19:06:27 Z Filya Magneti (COCI21_magneti) C++14
110 / 110
 75 ms 27264 KB 
/////////////////////include/////////////////////
//#include <bits/stdc++.h>
#define _CRT_SECURE_NO_WARNINGS
#include <iostream>
#include <string>
#include <algorithm>
#include <cmath>
#include <cstring>
#include <cassert>
#include <set>
#include <map>
#include <unordered_map>
#include <vector>
#include <stack>
#include <queue>
#include <iomanip>
#include <bitset>
#include <stdio.h>
#include <climits>
#include <numeric>
using namespace std;
//#include <ext/pb_ds/assoc_container.hpp>
//#include <ext/pb_ds/tree_policy.hpp>
//using namespace __gnu_pbds;
//template <typename T>
//using ordered_set = tree <T, null_type, less<T>, rb_tree_tag, tree_order_statistics_node_update>;
/////////////////////define/////////////////////
#define ci(x) if(x) cout << "YES" << '\n'; else cout << "NO" << '\n';
#define cii(x) if(check(x))
#define MOD 1000000007
#define MOD2 998244353
#define oo 1e9
#define ool 1e18L
#define pii pair<int, int>
#define pll pair<long long, long long>
#define mii map<int, int>
#define vi vector<int>
#define vpi vector<pair<int, int>>
#define vll vector <ll>
#define ff first
#define ss second
#define mp make_pair
#define ll long long
#define ld long double
#define pb push_back
#define eb emplace_back
#define pob pop_back
#define lb lower_bound
#define ub upper_bound
#define bs binary_search
#define sz(x) (int((x).size()))
#define all(x) (x).begin(), (x).end()
#define alll(x) (x), (x) + n
#define clr(x) (x).clear();
#define fri(x) for(int i = 0; i < x; ++i)
#define frj(x) for(int j = 0; j < x; ++j)
#define frp(x) for(int p = 0; p < x; ++p)
#define frr(a, b) for(int i = a; i < b; ++i)
#define frrj(a, b) for(int j = a; j < b; ++j)
#define fra(x) for(int i = 0; i < x; ++i) cin >> a[i];
#define frb(x) for(int i = 0; i < x; ++i) cin >> b[i];
#define frs(x) for(auto it = x.begin(); it != x.end(); ++it)
#define fr(x) for(auto it : x) //el
#define fastio ios_base::sync_with_stdio(false); cin.tie(0);
#define dbg(x) cerr << #x << ": " << x << endl;
#define ce(x) cout << x << endl;
#define uniq(x) x.resize(unique(all(x)) - x.begin()); //make all one after sorting
#define blt __builtin_popcount
/////////////////////print array, vector, deque, set, multiset, pair, map /////////////////////
void print(long long t) {cerr << t;}
void print(int t) {cerr << t;}
void print(string t) {cerr << t;}
void print(char t) {cerr << t;}
void print(double t) {cerr << t;}
void print(long double t) {cerr << t;}
void print(unsigned long long t) {cerr << t;}
template <class T, class V> void print(pair <T, V> p) {cerr << "{"; print(p.first); cerr << ","; print(p.second); cerr << "}";}
template <class T, class V> void print(T v[],V n) {cerr << "["; for(int i = 0; i < n; i++) {print(v[i]); cerr << " "; } cerr << "]"; cout << endl;}
template <class T> void print(vector <T> v) {cerr << "[ "; for (T i : v) {print(i); cerr << " ";} cerr << "]"; cout << endl;}
template <class T> void print(set <T> v) {cerr << "[ "; for (T i : v) {print(i); cerr << " ";} cerr << "]"; cout << endl;}
template <class T> void print(multiset <T> v) {cerr << "[ "; for (T i : v) {print(i); cerr << " ";} cerr << "]"; cout << endl;}
template <class T> void print(stack <T> v) {cerr << "[ "; stack<T> s = v; while(s.size()) {T i = s.top(); print(i); s.pop(); cerr << " ";} cerr << "]"; cout << endl;}
template <class T> void print(queue <T> v) {cerr << "[ "; queue<T> s = v; while(s.size()) {T i = s.front(); print(i); s.pop(); cerr << " ";} cerr << "]"; cout << endl;}
template <class T> void print(deque <T> v) {cerr << "[ "; deque<T> s = v; while(s.size()) {T i = s.front(); print(i); s.pop_front(); cerr << " ";} cerr << "]"; cout << endl;}
template <class T, class V> void print(map <T, V> v) {cerr << "[ "; for (auto i : v) {print(i); cerr << " ";} cerr << "]"; cout << endl;}
template <class T, class V> void print(unordered_map <T, V> v) {cerr << "[ "; for (auto i : v) {print(i); cerr << " ";} cerr << "]"; cout << endl;}
/////////////////////code/////////////////////
const int N=1e5;
ll fact[N+2], invFact[N+2];
long long power(long long a, int n) {
	if (n == 0) return 1;
	if (n % 2 == 0) {
		long long r = power(a, n / 2);
		return (r * r) % MOD;
	}
	return a * power(a, n - 1) % MOD;
}
ll cnk(int n, int k) {
    return (((fact[n] * invFact[k]) % MOD) * invFact[n - k]) % MOD;
}
void precalc() {
    fact[0] = 1;
    for(int i = 1; i <= N; i++)
        fact[i] = (fact[i - 1] * i) % MOD;
    invFact[N] = power(fact[N], MOD - 2); //(a / b) % p = a * (b ^ (p - 2) % p), 1 / f % p  = f ^ (p - 2) % p
    for(int i = N; i >= 1; i--) //(invF(i-1)^(p-2))%p=((invF(i-1)*i)^(p-2)%p*i)%p=((invF(i)^(p-2)%p)*(i^p-2%p))*i%p=((invF(i)^p-2%p)*(1/i%p)*i)%p=(invF(i)^p-2%p)%p
        invFact[i - 1] = (invFact[i] * i) % MOD;
}

ll a[55], dp[51][51][60002]; //i, j hat xumb ka, p hat tex enq zbaxacrel(datark texery radiusneri hety)
int main() {
    fastio;
    precalc();
    int n, l; 
    cin >> n >> l;
    fra(n) 
    sort(alll(a));
    dp[0][1][1] = 1; //1 voch te 3 vorovhetev heto vor urish ban avelacnenq inqy ed taracqy arden kzbaxacni
    fri(n-1) {
        for(ll j = n; j >= 1; j--) {
            for(int p = 1; p < l; p++) {
                if(!dp[i][j][p]) continue;
                dp[i + 1][j + 1][p + 1] += dp[i][j][p];
                dp[i + 1][j + 1][p + 1] %= MOD;
                dp[i + 1][j][p + a[i+1]] += dp[i][j][p] * 2 * j;
                dp[i + 1][j][p + a[i+1]] %= MOD;
                dp[i + 1][j-1][p + 2 * a[i+1] - 1] += dp[i][j][p] * (((j * (j - 1)) % MOD)); //blockeric 2y vercnenq u kpcnenq irar es tvov(hertakanutyuny kap uni)
                dp[i + 1][j-1][p + 2 * a[i+1] - 1] %= MOD;
            }
        }
    }
    ll ans = 0;
    for(int d = 1; d <= l; d++) {
        // dbg(d) dbg(dp[n-1][1][d]) dbg(n+l-d) dbg(l-d)
        ans += dp[n-1][1][d] * cnk(n+l-d, l-d); //l-d hat tex a mnacel
        ans %= MOD;
    }
    cout << ans;
    return 0;
}

//	           ♥ ♥ ♥  ♥  ♥    ♥   ♥    ♥
//	           ♥      ♥  ♥     ♥ ♥    ♥ ♥
//	           ♥ ♥ ♥  ♥  ♥      ♥    ♥   ♥
//	           ♥      ♥  ♥      ♥   ♥ ♥ ♥ ♥
//	           ♥      ♥  ♥ ♥ ♥  ♥  ♥       ♥
//
//        God loves Fil, Fil accepts God's will

Subtask #1
10.0 / 10.0

# Verdict Execution time Memory Grader output
1 Correct 49 ms 9484 KB Output is correct
2 Correct 5 ms 6996 KB Output is correct
3 Correct 4 ms 2260 KB Output is correct
4 Correct 2 ms 1876 KB Output is correct
5 Correct 5 ms 3284 KB Output is correct
6 Correct 14 ms 4408 KB Output is correct
7 Correct 23 ms 4816 KB Output is correct
8 Correct 2 ms 1876 KB Output is correct
9 Correct 12 ms 3668 KB Output is correct
10 Correct 2 ms 1876 KB Output is correct

Subtask #2
20.0 / 20.0

# Verdict Execution time Memory Grader output
1 Correct 4 ms 2644 KB Output is correct
2 Correct 3 ms 1876 KB Output is correct
3 Correct 4 ms 2772 KB Output is correct
4 Correct 2 ms 2004 KB Output is correct
5 Correct 4 ms 2772 KB Output is correct
6 Correct 2 ms 2260 KB Output is correct
7 Correct 5 ms 3540 KB Output is correct
8 Correct 3 ms 2388 KB Output is correct
9 Correct 2 ms 2084 KB Output is correct
10 Correct 2 ms 1876 KB Output is correct

Subtask #3
30.0 / 30.0

# Verdict Execution time Memory Grader output
1 Correct 5 ms 4888 KB Output is correct
2 Correct 4 ms 3284 KB Output is correct
3 Correct 5 ms 4820 KB Output is correct
4 Correct 3 ms 3412 KB Output is correct
5 Correct 4 ms 4820 KB Output is correct
6 Correct 2 ms 3028 KB Output is correct
7 Correct 4 ms 4948 KB Output is correct
8 Correct 3 ms 2324 KB Output is correct
9 Correct 3 ms 3028 KB Output is correct
10 Correct 2 ms 2260 KB Output is correct

Subtask #4
50.0 / 50.0

# Verdict Execution time Memory Grader output
1 Correct 49 ms 9484 KB Output is correct
2 Correct 5 ms 6996 KB Output is correct
3 Correct 4 ms 2260 KB Output is correct
4 Correct 2 ms 1876 KB Output is correct
5 Correct 5 ms 3284 KB Output is correct
6 Correct 14 ms 4408 KB Output is correct
7 Correct 23 ms 4816 KB Output is correct
8 Correct 2 ms 1876 KB Output is correct
9 Correct 12 ms 3668 KB Output is correct
10 Correct 2 ms 1876 KB Output is correct
11 Correct 4 ms 2644 KB Output is correct
12 Correct 3 ms 1876 KB Output is correct
13 Correct 4 ms 2772 KB Output is correct
14 Correct 2 ms 2004 KB Output is correct
15 Correct 4 ms 2772 KB Output is correct
16 Correct 2 ms 2260 KB Output is correct
17 Correct 5 ms 3540 KB Output is correct
18 Correct 3 ms 2388 KB Output is correct
19 Correct 2 ms 2084 KB Output is correct
20 Correct 2 ms 1876 KB Output is correct
21 Correct 5 ms 4888 KB Output is correct
22 Correct 4 ms 3284 KB Output is correct
23 Correct 5 ms 4820 KB Output is correct
24 Correct 3 ms 3412 KB Output is correct
25 Correct 4 ms 4820 KB Output is correct
26 Correct 2 ms 3028 KB Output is correct
27 Correct 4 ms 4948 KB Output is correct
28 Correct 3 ms 2324 KB Output is correct
29 Correct 3 ms 3028 KB Output is correct
30 Correct 2 ms 2260 KB Output is correct
31 Correct 61 ms 20880 KB Output is correct
32 Correct 38 ms 16612 KB Output is correct
33 Correct 58 ms 21120 KB Output is correct
34 Correct 13 ms 6300 KB Output is correct
35 Correct 60 ms 20424 KB Output is correct
36 Correct 5 ms 3028 KB Output is correct
37 Correct 75 ms 21708 KB Output is correct
38 Correct 14 ms 10580 KB Output is correct
39 Correct 65 ms 27264 KB Output is correct
40 Correct 16 ms 7252 KB Output is correct
41 Correct 59 ms 24144 KB Output is correct
42 Correct 3 ms 2260 KB Output is correct
43 Correct 50 ms 14428 KB Output is correct
44 Correct 8 ms 4820 KB Output is correct
45 Correct 49 ms 15948 KB Output is correct
46 Correct 2 ms 1876 KB Output is correct
47 Correct 18 ms 11236 KB Output is correct
48 Correct 14 ms 9460 KB Output is correct
49 Correct 3 ms 2388 KB Output is correct
50 Correct 3 ms 2260 KB Output is correct