QOJ.ac
QOJ
ID | Problem | Submitter | Result | Time | Memory | Language | File size | Submit time | Judge time |
---|---|---|---|---|---|---|---|---|---|
#251868 | #7775. 【模板】矩阵快速幂 | hos_lyric# | 10 | 186ms | 4124kb | C++14 | 6.1kb | 2023-11-15 11:15:25 | 2024-07-04 02:25:22 |
Judging History
answer
#include <cassert>
#include <cmath>
#include <cstdint>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <algorithm>
#include <bitset>
#include <complex>
#include <deque>
#include <functional>
#include <iostream>
#include <limits>
#include <map>
#include <numeric>
#include <queue>
#include <random>
#include <set>
#include <sstream>
#include <string>
#include <unordered_map>
#include <unordered_set>
#include <utility>
#include <vector>
using namespace std;
// using Int = long long;
template <class T1, class T2> ostream &operator<<(ostream &os, const pair<T1, T2> &a) { return os << "(" << a.first << ", " << a.second << ")"; };
template <class T> ostream &operator<<(ostream &os, const vector<T> &as) { const int sz = as.size(); os << "["; for (int i = 0; i < sz; ++i) { if (i >= 256) { os << ", ..."; break; } if (i > 0) { os << ", "; } os << as[i]; } return os << "]"; }
template <class T> void pv(T a, T b) { for (T i = a; i != b; ++i) cerr << *i << " "; cerr << endl; }
template <class T> bool chmin(T &t, const T &f) { if (t > f) { t = f; return true; } return false; }
template <class T> bool chmax(T &t, const T &f) { if (t < f) { t = f; return true; } return false; }
#define COLOR(s) ("\x1b[" s "m")
#ifndef LIBRA_OTHER_INT128_H_
#define LIBRA_OTHER_INT128_H_
#include <stdio.h>
#include <iostream>
constexpr unsigned __int128 toUInt128(const char *s) {
unsigned __int128 x = 0;
for (; *s; ++s) x = x * 10 + (*s - '0');
return x;
}
constexpr __int128 toInt128(const char *s) {
if (*s == '-') return -toInt128(s + 1);
__int128 x = 0;
for (; *s; ++s) x = x * 10 + (*s - '0');
return x;
}
unsigned __int128 inUInt128() {
static char buf[41];
scanf("%s", buf);
return toUInt128(buf);
}
__int128 inInt128() {
static char buf[41];
scanf("%s", buf);
return toInt128(buf);
}
void out(unsigned __int128 x) {
static char buf[41];
int len = 0;
do { buf[len++] = '0' + static_cast<int>(x % 10); } while (x /= 10);
for (int i = len; --i >= 0; ) putchar(buf[i]);
}
void out(__int128 x) {
if (x < 0) {
putchar('-');
out(-static_cast<unsigned __int128>(x));
} else {
out(static_cast<unsigned __int128>(x));
}
}
std::ostream &operator<<(std::ostream &os, unsigned __int128 x) {
static char buf[41];
int len = 0;
do { buf[len++] = '0' + static_cast<int>(x % 10); } while (x /= 10);
for (int i = len; --i >= 0; ) os << buf[i];
return os;
}
std::ostream &operator<<(std::ostream &os, __int128 x) {
if (x < 0) {
os << '-' << -static_cast<unsigned __int128>(x);
} else {
os << static_cast<unsigned __int128>(x);
}
return os;
}
#endif // LIBRA_OTHER_INT128_H_
using Int = __int128;
constexpr Int INF = toInt128("1002003004005006007006005004003002001");
constexpr unsigned MO = 998244353;
char buf[110];
int S;
int N, M;
string K;
vector<int> A, B;
vector<Int> C;
namespace small {
vector<int> run() {
cerr<<"[small::run]"<<endl;
const int K0 = stoi(K);
vector<Int> crt(N, INF);
crt[0] = 0;
for (int k = 0; k < K0; ++k) {
vector<Int> nxt(N, INF);
for (int i = 0; i < M; ++i) {
chmin(nxt[B[i]], crt[A[i]] + C[i]);
}
crt.swap(nxt);
}
vector<int> ans(N, -1);
for (int u = 0; u < N; ++u) if (crt[u] < INF) {
ans[u] = crt[u] % MO;
}
return ans;
}
} // small
namespace large {
vector<int> run() {
cerr<<"[large::run]"<<endl;
vector<vector<int>> G(N);
for (int i = 0; i < M; ++i) {
G[A[i]].push_back(i);
}
vector<pair<pair<Int, int>, int>> cycs;
for (int s = 0; s < N; ++s) {
Int pm = 1;
int qm = 0;
vector<Int> crt(N, INF);
crt[s] = 0;
for (int k = 1; k <= N; ++k) {
vector<Int> nxt(N, INF);
for (int i = 0; i < M; ++i) {
chmin(nxt[B[i]], crt[A[i]] + C[i]);
}
crt.swap(nxt);
Int p = crt[s];
int q = k;
if (p < INF) {
if (pm * q > p * qm) {
pm = p;
qm = q;
}
}
}
if (qm) {
cycs.emplace_back(make_pair(pm, qm), s);
}
}
sort(cycs.begin(), cycs.end());
const int cycsLen = cycs.size();
// cerr<<"cycs = "<<cycs<<endl;
// TODO
const Int K0 = toInt128(K.c_str());
vector<Int> mns(N, INF);
for (int j0 = 0, j1 = 0; j0 < cycsLen; j0 = j1) {
for (; j1 < cycsLen && cycs[j0].first == cycs[j1].first; ++j1) {}
const Int p = cycs[j0].first.first;
const int q = cycs[j0].first.second;
vector<int> on(N, 0);
for (int j = j0; j < j1; ++j) {
on[cycs[j].second] = 1;
}
const Int c0 = (K0 / q) * p;
const int r0 = K0 % q;
vector<Int> here(N, INF);
vector<Int> crt(N << 1, INF);
crt[0 << 1 | on[0]] = 0;
for (int k = 0; k < N + q * N; ++k) {
if (k % q == r0) {
for (int u = 0; u < N; ++u) if (crt[u << 1 | 1] < INF) {
chmin(here[u], crt[u << 1 | 1] - (k / q) * p);
}
}
vector<Int> nxt(N << 1, INF);
for (int i = 0; i < M; ++i) {
chmin(nxt[B[i] << 1 | on[B[i]]], crt[A[i] << 1 | 0] + C[i]);
chmin(nxt[B[i] << 1 | 1], crt[A[i] << 1 | 1] + C[i]);
}
crt.swap(nxt);
}
for (int u = 0; u < N; ++u) if (here[u] < INF) {
chmin(mns[u], c0 + here[u]);
}
}
vector<int> ans(N, -1);
for (int u = 0; u < N; ++u) if (mns[u] < INF) {
ans[u] = mns[u] % MO;
}
return ans;
}
} // large
int main() {
for (int numCases; ~scanf("%d%d", &S, &numCases); ) { for (int caseId = 1; caseId <= numCases; ++caseId) {
scanf("%d%d%s", &N, &M, buf);
K = buf;
A.resize(M);
B.resize(M);
C.resize(M);
for (int i = 0; i < M; ++i) {
scanf("%d%d", &A[i], &B[i]);
--A[i];
--B[i];
C[i] = inInt128();
}
vector<int> ans;
if (K.size() <= 9 && stoi(K) <= N * N) {
ans = small::run();
} else {
ans = large::run();
}
for (int u = 0; u < N; ++u) {
if (u) printf(" ");
printf("%d", ans[u]);
}
#ifdef LOCAL
printf(" ");
#endif
puts("");
}
#ifndef LOCAL
break;
#endif
}
return 0;
}
Details
Tip: Click on the bar to expand more detailed information
Subtask #1:
score: 10
Accepted
Test #1:
score: 10
Accepted
time: 3ms
memory: 4124kb
input:
1 1 100 101 899539897889989959 74 35 910832669819965536 35 85 910832669819965536 85 88 910832669819965536 88 30 910832669819965536 30 58 910832669819965536 58 60 910832669819965536 60 34 910832669819965536 34 8 910832669819965536 8 67 910832669819965536 67 89 910832669819965536 89 32 910832669819965...
output:
395495792 395495781 395495783 395495839 395495793 395495789 395495754 395495832 395495845 395495755 395495823 395495773 395495753 395495800 395495782 395495763 395495847 395495761 395495794 395495791 395495786 395495821 395495798 395495765 395495772 395495770 395495757 395495819 395495843 395495828 ...
result:
ok 100 numbers
Test #2:
score: 0
Accepted
time: 157ms
memory: 3840kb
input:
1 1 100 200 998858598565699977 89 61 596014036562538539 89 84 921297646113897322 61 84 946923234442637386 61 35 641628261157284465 84 35 979893473772327497 84 78 700172488379560611 35 78 963617193748189613 35 54 951598888254521423 78 54 680825215292116806 78 21 737055858973038555 54 21 7491794406112...
output:
590375247 265938345 203065828 597548045 369717762 226160283 377877020 360218254 956162456 408060901 387231165 759578975 67601808 790211315 608425007 343195480 177353482 436533546 717630459 417099733 542227025 861764246 913806375 587268602 989846681 435016550 66609901 817090566 256847656 844441854 94...
result:
ok 100 numbers
Test #3:
score: 0
Accepted
time: 179ms
memory: 3784kb
input:
1 1 100 181 348568663892999968 25 19 990622898175774733 19 94 871060999389241529 94 24 969317630558501400 24 43 908457844888427461 43 52 816088481082287266 52 62 978618931332609685 62 99 761714433396732044 99 85 741344935503895668 85 64 964684335126604843 64 69 988098065125373655 69 31 7506975506815...
output:
916998469 916998469 916998469 76035207 62461893 916998469 389136594 916998469 916998469 173423529 164423356 822964468 626456020 916998469 744111524 916998469 398953850 916998469 342238577 916998469 255074799 784015663 916998469 740933556 587088671 811719512 916998469 916998469 916998469 916998469 14...
result:
ok 100 numbers
Test #4:
score: 0
Accepted
time: 186ms
memory: 3840kb
input:
1 1 100 189 295064635124730243 18 50 754672892083203214 50 88 962632394366987404 88 15 906700334097319336 15 26 967741400981618572 26 91 996214498763867892 91 35 882157548994344280 35 68 983621159612138407 68 51 563935036482744182 51 75 991205513962219551 75 72 974025375183814852 72 11 7979447663592...
output:
663199381 739882534 663199381 28600701 663199381 944601671 836329160 894091561 629507606 663199381 246830507 663199381 491987421 663199381 802123884 663199381 663199381 663199381 414785533 989396289 663199381 663199381 663199381 663199381 663199381 663199381 663199381 663199381 663199381 663199381 4...
result:
ok 100 numbers
Test #5:
score: 0
Accepted
time: 24ms
memory: 3936kb
input:
1 254 40 74 997173688939799978 38 6 890721839505665075 6 10 992308491267087930 10 29 960202932780090595 29 20 952827125924298715 20 34 868314670055961466 34 31 756448635709788087 31 14 857625921909632516 14 18 917667459973696862 18 21 985939328882662624 21 1 734882468602343649 1 11 66102593854575036...
output:
177014577 177014577 177014577 885341552 472856470 177014577 363547548 177014577 499847464 653076748 177014577 177014577 177014577 177014577 487939796 177014577 213466543 586729345 244952763 177014577 177014577 177014577 177014577 890105934 177014577 177014577 890105934 177014577 177014577 798890006 ...
result:
ok 3575 numbers
Test #6:
score: 0
Accepted
time: 24ms
memory: 3848kb
input:
1 356 32 47 967844399484634837 4 30 776954643355911997 30 20 811634053140142741 20 22 747630229183579429 22 2 806282875388761050 2 26 719793351534499411 26 17 797537828929335673 17 24 890423236992687627 24 21 970792227007588899 21 8 850078803097295262 8 15 958474507028658347 15 1 972636122087215360 ...
output:
890097469 525779071 636798453 776362497 776362497 687961593 158033324 776362497 345910504 380622623 239804834 440670451 137231885 985041116 222869127 137231885 705696901 637534644 347889826 696528073 291555427 146553026 776362497 624486185 137231885 642408114 520519927 137231885 438373632 263924254 ...
result:
ok 4784 numbers
Subtask #2:
score: 0
Time Limit Exceeded
Test #7:
score: 0
Time Limit Exceeded
input:
2 1 300 598 8179377797889487867988994778539839593376697796496698959964978969 1 2 977880533270721156 2 1 977880533270721156 2 3 977880533270721156 3 2 977880533270721156 3 4 977880533270721156 4 3 977880533270721156 4 5 977880533270721156 5 4 977880533270721156 5 6 977880533270721156 6 5 977880533270...
output:
result:
Subtask #3:
score: 0
Skipped
Dependency #2:
0%
Subtask #4:
score: 0
Wrong Answer
Dependency #1:
100%
Accepted
Test #19:
score: 0
Wrong Answer
time: 6ms
memory: 3916kb
input:
4 1 100 101 6888995999928874698772868926699656683388498575797893294688976887 25 90 495511874996847106 90 84 495511874996847106 84 82 495511874996847106 82 40 495511874996847106 40 97 495511874996847106 97 5 495511874996847106 5 24 495511874996847106 24 16 495511874996847106 16 19 495511874996847106 ...
output:
-242181964 -242181971 -242181914 -242181955 -242181926 -242181981 -242181950 -242181900 -242181972 -242181962 -242181933 -242181903 -242181974 -242181957 -242181910 -242181924 -242181939 -242181952 -242181923 -242181951 -242181963 -242181922 -242181986 -242181925 -242181932 -242181997 -242181992 -24...
result:
wrong answer 1st numbers differ - expected: '662900138', found: '-242181964'
Subtask #5:
score: 0
Time Limit Exceeded
Dependency #1:
100%
Accepted
Test #25:
score: 15
Accepted
time: 143ms
memory: 3888kb
input:
5 1 300 301 969767789936486493 164 284 964646444984408140 284 241 964646444984408140 241 281 964646444984408140 281 138 964646444984408140 138 242 964646444984408140 242 112 964646444984408140 112 217 964646444984408140 217 170 964646444984408140 170 31 964646444984408140 31 300 964646444984408140 3...
output:
562333388 562333371 562333450 562333457 562333181 562333366 562333433 562333276 562333204 562333354 562333361 562333374 562333436 562333405 562333369 562333286 562333360 562333318 562333396 562333251 562333480 562333220 562333333 562333460 562333359 562333295 562333293 562333335 562333402 562333226 ...
result:
ok 300 numbers
Test #26:
score: -15
Time Limit Exceeded
input:
5 1 300 600 798876399989994933 7 196 978372754397099680 7 150 850366341978113658 196 150 741178931696536015 196 241 918555502737513857 150 241 755464499814711391 150 249 715712249601810459 241 249 834572033520725671 241 172 840925258261612828 249 172 765221764158211117 249 92 987381804975984305 172 ...
output:
result:
Subtask #6:
score: 0
Skipped
Dependency #3:
0%