結果
問題 | No.1139 Slime Race |
ユーザー |
![]() |
提出日時 | 2020-07-31 21:57:57 |
言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 79 ms / 2,000 ms |
コード長 | 3,512 bytes |
コンパイル時間 | 5,069 ms |
コンパイル使用メモリ | 345,888 KB |
実行使用メモリ | 6,944 KB |
最終ジャッジ日時 | 2024-07-02 08:21:45 |
合計ジャッジ時間 | 7,099 ms |
ジャッジサーバーID (参考情報) |
judge5 / judge4 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 22 |
ソースコード
#include <iostream>#include <random>#include <bits/stdc++.h>#include <boost/multiprecision/cpp_int.hpp>using namespace std;namespace mp = boost::multiprecision;using namespace mp;using ull = __int128;using ll = long long;using cll = cpp_int;#define rep(i,a) for(int i=0;i<a;i++)const int mod = 1000000007;const int MOD = 1000000007;const int MAX = 50000;long long fac[MAX],finv[MAX],inv[MAX];// テーブルを作る前処理void COMinit() {fac[0] = fac[1] = 1;finv[0] = finv[1] = 1;inv[1] = 1;for (int i = 2; i < MAX; i++){fac[i] = fac[i - 1] * i % MOD;inv[i] = MOD - inv[MOD%i] * (MOD / i) % MOD;finv[i] = finv[i - 1] * inv[i] % MOD;}}// 二項係数計算long long COM(int n, int k){if (n < k) return 0;if (n < 0 || k < 0) return 0;return fac[n] * (finv[k] * finv[n - k] % MOD) % MOD;}class mint {long long x;public:mint(long long x=0) : x((x%mod+mod)%mod) {}mint operator-() const {return mint(-x);}mint& operator+=(const mint& a) {if ((x += a.x) >= mod) x -= mod;return *this;}mint& operator-=(const mint& a) {if ((x += mod-a.x) >= mod) x -= mod;return *this;}mint& operator*=(const mint& a) {(x *= a.x) %= mod;return *this;}mint operator+(const mint& a) const {mint res(*this);return res+=a;}mint operator-(const mint& a) const {mint res(*this);return res-=a;}mint operator*(const mint& a) const {mint res(*this);return res*=a;}mint pow(ll t) const {if (!t) return 1;mint a = pow(t>>1);a *= a;if (t&1) a *= *this;return a;}// for prime modmint inv() const {return pow(mod-2);}mint& operator/=(const mint& a) {return (*this) *= a.inv();}mint operator/(const mint& a) const {mint res(*this);return res/=a;}friend ostream& operator<<(ostream& os, const mint& m){os << m.x;return os;}friend istream& operator>>(istream& lhs,mint& rhs) noexcept {lhs >> rhs.x;return lhs;}};ll modmul(ll x, ll y, ll mod){ll ret = 0;while(y){if(y & 1){ret = (ret + x) % mod;}x = x * 2 % mod;y >>= 1;}return ret;}ll modpow(ll x, ll n, ll mod){ll ret = 1;while(n){if(n & 1){ret = modmul(ret, x, mod);}x = modmul(x, x, mod);n >>= 1;}return ret;}int kkk(vector<int> n){ll b=n.size();if(b==1){return n.at(0);}else{ll ans=0;vector<vector<int>> data(b-1, vector<int>(b-1));for(int i=0;i<b-1;i++){for(int j=0;j<b-1;j++){if(i==j){data.at(i).at(j)=n.at(j)+n.at(j+1);}else if(i>j){data.at(i).at(j)=n.at(j);}else{data.at(i).at(j)=n.at(j+1);}}}vector<int> k(b-1);for(int i=0;i<b-1;i++){k.at(i)=n.at(i);}ans=(ans+(n.at(b-1)%MOD)*(kkk(k)%MOD)%MOD);for(int i=0;i<b-1;i++){ans=(ans+((n.at(i)%MOD)*(kkk(data.at(i))%MOD))%MOD)%MOD;}return ans;}}int main() {// 前処理COMinit();ll a,b;cin>>a>>b;vector<ll> n(a);vector<ll> m(a);for(int i=0;i<a;i++){cin>>n.at(i);}ll ans=0;for(int i=0;i<a;i++){cin>>m.at(i);ans=ans+m.at(i);}if(b%ans==0){cout<<b/ans<<endl;}else{cout<<b/ans+1<<endl;}}