Submission #00314
ソースコード
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 145 146 147 148 149 150 151 152 153 154 155 | #include <bits/stdc++.h> using namespace std; #define fs first #define sc second #define pb emplace_back #define mp make_pair #define all(v) v.begin(), v.end() #define rall(v) v.rbegin(), v.rend() #define rep(i, n) for (int i = 0; i < (int)(n); ++i) using pii = pair< int , int >; using vi = vector< int >; using lint = long long ; const int inf = 1001001001; const lint linf = 1001001001001001001ll; const int mod = 1e9 + 7; const int dx[]{0, 1, 0, -1, -1, -1, 1, 1}, dy[]{1, 0, -1, 0, -1, 1, -1, 1}; template < typename T> inline bool chmin(T &a, T b) { if (a > b) { a = b; } return a > b; } template < typename T> inline bool chmax(T &a, T b) { if (a < b) { a = b; } return a < b; } template < typename T> inline void print( const T &x, string s = "\n" ) { cout << x << s; } template < typename T> inline void print( const vector<T> &v, string s = " " ) { if (!v.size()) puts ( "" ); rep(i, v.size()) cout << v[i] << (i + 1 == v.size() ? "\n" : s); } inline bool inside( int y, int x, int H, int W) { return 0 <= y && y < H && 0 <= x && x < W; } inline lint in() { lint x; std::cin>>x; return x; } typedef vector<vector<lint> > matrix; matrix e, A; int n; vi a, b; lint mod_pow(lint a, lint x, int p = 1e9 + 7) { lint ret = 1; for (lint i = 1ll << 60; i; i >>= 1) { ret *= ret; ret %= p; if (x & i) { ret *= a; ret %= p; } } return ret; } vector<lint> fact(1, 1); lint Fact(lint n, int p = 1e9 + 7) { if (n == 0) return 1; if (fact.size() <= n) { fact.push_back(Fact(n - 1) * n % p); } return fact[n]; } void makee( int n) { e.clear(); e.resize(n); rep(i, n) rep(j, n) { lint E = 1; lint Z = 0; e[i].push_back(i == j ? E : Z); } } void makeA( int n) { A.clear(); A.resize(n); rep(i, n) A[i].resize(n); A[0][0] = 1; A[0][1] = 0; A[1][0] = 1; A[1][1] = 10; } matrix product(matrix a, matrix b) { matrix ret; ret.resize(a.size()); for ( int i = 0; i < a.size(); i++) ret[i].resize(b[0].size()); for ( int i = 0; i < a.size(); i++) { for ( int j = 0; j < b[0].size(); j++) { for ( int k = 0; k < a[0].size(); k++) { ret[i][j] += a[i][k] * b[k][j]; if (k == ((k >> 3) << 3)) ret[i][j] %= mod; } ret[i][j] %= mod; } } return ret; } matrix po(matrix a, lint b) { if (b == 0) return e; matrix z = po(a, b / 2); z = product(z, z); if (b & 1) z = product(z, a); return z; } lint number( int a, lint b) { matrix x(2); x[0].pb(a); x[1].pb(0); matrix ans = product(po(A, b), x); return ans[1][0]; } signed main() { cin >> n; rep(i, n) { a.pb(in()); b.pb(in()); } vector< int > num[10]; rep(i, n) { num[a[i]].pb(b[i]); } rep(i, 10) sort(all(num[i])); makee(2); makeA(2); lint res1 = 0; lint res2 = 1; // 使われていない最小のインデックス int used[10] = {0}; // 1以上の最小の数字を前に for ( int i = 1; i < 10; ++i) { if (num[i].size()) { res1 = number(i, num[i][0]); used[i] = 1; break ; } } if (res1 == 0) { res2 = 1; } else { rep(i, 10) { lint b = 0; for ( int j = used[i]; j < num[i].size(); ++j) { b += num[i][j]; } res1 = (res1 * mod_pow(10, b)) % mod; res1 += number(i, b); res1 %= mod; res2 *= Fact(num[i].size() - used[i]); res2 %= mod; } } print(res1); print(res2); } |
ステータス
項目 | データ |
---|---|
問題 | 0003 - repdigit |
ユーザー名 | Luzzle0119 |
投稿日時 | 2017-07-07 22:41:50 |
言語 | C++11 |
状態 | Wrong Answer |
得点 | 0 |
ソースコード長 | 3471 Byte |
最大実行時間 | 63 ms |
最大メモリ使用量 | 7660 KB |
セット
セット | 得点 | Cases | |
---|---|---|---|
1 | ALL | 0 / 100 | * |
テストケース
ファイル名 | 状態 | 実行時間 | メモリ使用量 | # |
---|---|---|---|---|
01_sample1.in | AC | 13 ms | 480 KB |
1
|
01_sample2.in | AC | 12 ms | 448 KB |
1
|
01_sample3.in | AC | 13 ms | 544 KB |
1
|
02_handmake1.in | AC | 16 ms | 520 KB |
1
|
02_handmake2.in | WA | 14 ms | 492 KB |
1
|
02_handmake3.in | AC | 11 ms | 464 KB |
1
|
02_handmake4.in | AC | 13 ms | 436 KB |
1
|
02_handmake5.in | AC | 16 ms | 532 KB |
1
|
02_handmake6.in | AC | 13 ms | 504 KB |
1
|
02_handmake7.in | AC | 15 ms | 472 KB |
1
|
02_handmake8.in | WA | 13 ms | 568 KB |
1
|
02_handmake9.in | AC | 13 ms | 668 KB |
1
|
03_random1.in | WA | 11 ms | 636 KB |
1
|
03_random2.in | WA | 13 ms | 604 KB |
1
|
03_random3.in | WA | 10 ms | 452 KB |
1
|
03_random4.in | WA | 13 ms | 548 KB |
1
|
03_random5.in | WA | 13 ms | 648 KB |
1
|
03_random6.in | WA | 14 ms | 612 KB |
1
|
04_random1.in | AC | 13 ms | 584 KB |
1
|
04_random2.in | AC | 12 ms | 680 KB |
1
|
04_random3.in | AC | 14 ms | 652 KB |
1
|
04_random4.in | WA | 13 ms | 620 KB |
1
|
04_random5.in | WA | 13 ms | 716 KB |
1
|
04_random6.in | WA | 12 ms | 564 KB |
1
|
05_random1.in | AC | 52 ms | 2564 KB |
1
|
05_random2.in | AC | 58 ms | 2456 KB |
1
|
05_random3.in | AC | 55 ms | 2724 KB |
1
|
05_random4.in | AC | 63 ms | 2584 KB |
1
|
06_random1.in | AC | 60 ms | 7660 KB |
1
|
06_random3.in | WA | 60 ms | 2004 KB |
1
|