結果

問題 No.603 hel__world (2)
ユーザー Sebastian KingSebastian King
提出日時 2017-12-03 14:00:15
言語 C++11
(gcc 11.4.0)
結果
TLE  
実行時間 -
コード長 2,551 bytes
コンパイル時間 1,605 ms
コンパイル使用メモリ 174,740 KB
実行使用メモリ 30,276 KB
最終ジャッジ日時 2024-05-06 05:24:10
合計ジャッジ時間 7,034 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge4
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 12 ms
25,760 KB
testcase_01 AC 11 ms
18,944 KB
testcase_02 AC 12 ms
18,944 KB
testcase_03 AC 13 ms
18,944 KB
testcase_04 AC 11 ms
18,944 KB
testcase_05 AC 12 ms
18,816 KB
testcase_06 AC 12 ms
19,024 KB
testcase_07 AC 12 ms
18,944 KB
testcase_08 AC 12 ms
18,816 KB
testcase_09 AC 13 ms
18,980 KB
testcase_10 AC 13 ms
18,944 KB
testcase_11 AC 203 ms
30,196 KB
testcase_12 AC 93 ms
29,192 KB
testcase_13 AC 204 ms
30,276 KB
testcase_14 AC 13 ms
19,072 KB
testcase_15 TLE -
testcase_16 -- -
testcase_17 -- -
testcase_18 -- -
testcase_19 -- -
testcase_20 -- -
testcase_21 -- -
testcase_22 -- -
testcase_23 -- -
testcase_24 -- -
testcase_25 -- -
testcase_26 -- -
testcase_27 -- -
testcase_28 -- -
testcase_29 -- -
testcase_30 -- -
testcase_31 -- -
testcase_32 -- -
権限があれば一括ダウンロードができます
コンパイルメッセージ
main.cpp: In function ‘ll gao(ll, std::vector<int>)’:
main.cpp:81:47: warning: narrowing conversion of ‘b.std::vector<long long int>::size()’ from ‘std::vector<long long int>::size_type’ {aka ‘long unsigned int’} to ‘int’ [-Wnarrowing]
   81 |                 Q.push((node){y + 1, x, b.size()});
      |                                         ~~~~~~^~

ソースコード

diff # 

#include<bits/stdc++.h>

using namespace std;

typedef long long ll;
typedef pair<int, int> PII;

const int MM = 1e6 + 3;
const double eps = 1e-8;
const int MAXN = 2e6 + 10;
using int64 = long long;
using int128 = __int128;
const int N = 1e6 + 10;
const int mod = 1e6 + 3;
const int128 denom = 2e20;

int64 fac[mod], inv[mod];

int64 binom(int64 n, int64 m) {
  if (n < m) return 0;
  if (n == m || m == 0) return 1;
  int64 a = n % mod, b = m % mod;
  if (a < b) return 0;
  return fac[a] * inv[fac[b]] * inv[fac[a - b]] % mod * binom(n / mod, m / mod) % mod;
}

ll a[100];
vector<int> s[100];
string str;

int n, m;

void prework(){
  fac[0] = fac[1] = inv[1] = 1;
  for (int i = 2; i < mod; ++i) {
    fac[i] = i * fac[i - 1] % mod;
    inv[i] = (mod - mod / i) * inv[mod % i] % mod;
  }
}

void read(){}

struct node{
	__int128 y, x;
	int num;
};

bool operator < (const node & a, const node & b){
	return a.y * (b.y - b.x) < b.y * (a.y - a.x);
}

ll gao(ll S, vector<int> a){
	//cout << S << ' ' << a.size() << endl;
	ll __sm = 0;
	for (auto &x : a) __sm += x;
	sort(a.begin(), a.end());
	if (__sm > S) return 0ll;
	if (__sm == S) return 1ll;
	if (a.size() == 0) return 1ll;
	if (a.size() == 1)
		return binom(S, a[0]);
	__int128 limit = 1e20;
	__int128 l = limit + 1, r = limit * (a.back() + 2) + 1;
	while(l + 1 < r){
		__int128 m = (r + l) / 2;
		__int128 tot = 0;
		for (auto &x : a){
			tot += x * m / (m - limit);
			if (tot > S) break;
		}
		if (tot > S)
			l = m;
		else
			r = m;
	//	cout << (ll)tot << ' ' << S << endl;
	}
	vector<ll> b;
	priority_queue<node> Q;
	for (auto &x : a){
		ll y = (ll)(x * r / (r - limit));
		Q.push((node){y + 1, x, b.size()});
		b.push_back(y);
		S -= y;
	}
	assert(S >= 0);
	//cout << "!!!" << S << endl;
	while(S){
		S--;
		node p = Q.top();
		Q.pop();
		b[p.num]++;
		p.y++;
		Q.push(p);
	}
	ll ret = 1;
	for (int i = 0; i < a.size(); i++)
		ret = ret * binom(b[i], a[i]) % mod;
	//cout << "?? " << S << ' ' << a.size() << ' ' << ret << endl;
	return ret;
}

void solve(int casi){
//	cout << "Case #" << casi << ": ";
	for (int i = 0; i < 26; i++)
		cin >> a[i];
	cin >> str;
	str += '+';
	for (int i = 1, j = 0; i < str.length(); i++){
		if (str[i] != str[j]){
			s[str[j]-'a'].push_back(i - j);
			j=i;
		}
	}
	ll ans = 1;
	for (int i = 0; i < 26; i++)
		ans = ans * gao(a[i], s[i]) % MM;
	cout << ans << endl;
}

void printans(){

}


int main(){
	std::ios::sync_with_stdio(false);
	prework();
	int T = 1;
//	cin>>T;
	for(int i = 1; i <= T; i++){
		read();
		solve(i);
		printans();
	}
	return 0;
}
0