QOJ.ac
QOJ
ID | 题目 | 提交者 | 结果 | 用时 | 内存 | 语言 | 文件大小 | 提交时间 | 测评时间 |
---|---|---|---|---|---|---|---|---|---|
#406265 | #7858. Basic Equation Solving | dnialh | AC ✓ | 319ms | 3940kb | C++20 | 12.9kb | 2024-05-07 03:47:41 | 2024-05-07 03:47:43 |
Judging History
answer
#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
typedef long double ld;
typedef pair<int, int> pi;
typedef pair<ll,ll> pl;
typedef vector<int> vi;
typedef vector<ll> vl;
typedef vector<pi> vpi;
typedef vector<pl> vpl;
#define F0R(i,n) for (int i = 0; i < n; i++)
#define FOR(i,a,b) for (int i = a; i <= b; i++)
#define F0Rd(i,a) for (int i = (a)-1; i >= 0; i--)
#define FORd(i,a,b) for (int i = (b); i >= (a); i--)
#define trav(a, x) for (auto& a : x)
#define rep(i, a, b) for(int i = a; i < (b); ++i)
#define f first
#define s second
#define mp make_pair
#define pb push_back
#define ins insert
#define lb lower_bound
#define ub upper_bound
#define sz(x) (int)x.size()
#define all(x) x.begin(), x.end()
const char nl = '\n';
const int MAX_N = 100011;
const ll INF = (1<<29) + 123;
const ll MOD = 998244353;
const ld PI = 4*atan((ld)1);
template <typename T> bool ckmin(T& a, const T& b) { return a > b ? a=b, 1 : 0; }
template <typename T> bool ckmax(T& a, const T& b) { return b > a ? a=b, 1 : 0; }
/**** Credit to chatgpt 4.0 ****/
// Stream operator for std::pair
template<typename T1, typename T2>
ostream& operator<<(ostream &out, const pair<T1, T2> &v) {
out << "(" << v.first << ", " << v.second << ")";
return out;
}
// Trait to check if a type is iterable
template<typename T, typename = void>
struct is_iterable : false_type {};
template<typename T>
struct is_iterable<T, void_t<decltype(begin(declval<T>())), decltype(end(declval<T>()))>> : true_type {};
// Stream operator for iterable types excluding std::string
template<typename TT>
typename enable_if<is_iterable<TT>::value && !is_same<TT, string>::value, ostream&>::type
operator<<(ostream& out, const TT& c) {
out << "{ ";
for (const auto& x : c) out << x << " ";
out << "}";
return out;
}
template<typename T>
ostream& operator<<(ostream& out, std::stack<T> container) {
std::vector<T> elements;
while (!container.empty()) {
elements.push_back(container.top());
container.pop();
}
std::reverse(elements.begin(), elements.end()); // Reverse to maintain order
return out << elements;
}
template<typename T>
ostream& operator<<(ostream& out, std::queue<T> container) {
std::vector<T> elements;
while (!container.empty()) {
elements.push_back(container.front());
container.pop();
}
return out << elements;
}
// Helper function to print std::priority_queue
template<typename T, typename Container, typename Compare>
ostream& operator<<(ostream& out, std::priority_queue<T, Container, Compare> pq) {
out << "{";
while (!pq.empty()) {
out << " " << pq.top();
pq.pop();
}
out << " }";
return out;
}
#ifdef DBG
void dbg_out() { cerr << endl; }
template<typename Head, typename... Tail>
void dbg_out(Head H, Tail... T) {
cerr << ' ' << H;
dbg_out(T...);
}
#define dbg(...) cerr << #__VA_ARGS__ << ":", dbg_out(__VA_ARGS__);
#define dbg_array(a, n) cerr << #a << ": { "; for(int i = 0; i < n; i++) cerr << a[i] << " "; cerr << "}\n";
#else
#define dbg(...)
#define dbg_array(a, n)
#endif
mt19937_64 rng(chrono::steady_clock::now().time_since_epoch().count());
const int MX = 3e5+5;
pair<pair<string, string>, char> p[11];
int n;
struct State {
int gr[26][26]; // 0 is not set, 1 is equal, 2 is i > j, 3 is i < j
int more[26];
int less[26];
State() {
F0R(i, 26) F0R(j, 26) gr[i][j] = 0;
F0R(i, 26) more[i] = -1;
F0R(i, 26) less[i] = 10;
}
};
bool isDigit(char c) {
return c >= '0' && c <= '9';
}
ll ans = 0;
vi comp;
int vis[26];
void dfs(State& s, int v) {
if (vis[v]) return;
comp.pb(v);
vis[v] = 1;
F0R(i, 26) {
if (s.gr[v][i] != 0) dfs(s, i);
}
}
int par[26], SZ[26];
int find(int v) {
return (v==par[v]?v:par[v] = find(par[v]));
}
void merge(int a, int b) {
a = find(a), b = find(b);
if (a == b) return;
if (SZ[a] < SZ[b]) swap(a, b);
par[b] = a;
SZ[a] += SZ[b];
}
int vis2[26];
int dfs2(vector<vi>& gr, int v, int mark) {
if (vis2[v] == mark) return 0;
vis2[v] = mark;
int ret = (1<<v);
trav(u, gr[v]) {
ret |= dfs2(gr, u, mark);
}
return ret;
}
ll calc(State& s) {
dbg("calcing ans for state");
// F0R(i, 26) {
// F0R(j, 26) cerr << s.gr[i][j] << " ";
// cerr << nl;
// }
// F0R(i, 26) cerr << s.more[i] << " "; cerr << nl;
// F0R(i, 26) cerr << s.less[i] << " "; cerr << nl;
F0R(i, 26) vis[i] = 0;
// now for the graph of 26
// find each component
ll ret = 1;
F0R(i, 26) {
if (vis[i] == 1) continue;
F0R(i, 26) {
par[i] = i;
SZ[i] = 1;
}
comp.clear();
dfs(s, i);
// now for this component, merge things
trav(u, comp) {
trav(v, comp) {
if (u == v) continue;
if (s.gr[u][v] == 1) {
merge(u, v);
}
}
}
// now relabel it to be 0, 1, 2...
vi labels(26, -1);
int cnt = 0;
trav(u, comp) if (find(u) == u) {
dbg((char)(u+'A'), cnt);
labels[u] = cnt++;
}
vector<vi> gr(cnt);
trav(u, comp) {
trav(v, comp) {
if (find(u) == find(v) && s.gr[u][v] >= 2) return 0;
int a = labels[find(u)], b = labels[find(v)];
if (s.gr[u][v] == 2) {
gr[a].pb(b);
dbg("edge", a, b);
} else if (s.gr[u][v] == 3) {
gr[b].pb(a);
dbg("edge", b, a);
}
}
}
dbg(i, cnt);
// find all subtree masks for each number
int subtreeMask[cnt];
F0R(i, 26) vis2[i] = -1;
F0R(i, cnt) {
subtreeMask[i] = dfs2(gr, i, i);
dbg(i, subtreeMask[i]);
}
bool downwardClosed[1<<cnt]; F0R(i, (1<<cnt)) downwardClosed[i] = 0;
F0R(i, 1<<cnt) {
int mask = 0;
F0R(j, cnt) if (i&(1<<j)) mask |= subtreeMask[j];
downwardClosed[mask] = 1;
dbg(mask);
}
// create validity mask for each level
int mergedMore[cnt]; F0R(i, cnt) mergedMore[i] = -1;
int mergedLess[cnt]; F0R(i, cnt) mergedLess[i] = 10;
trav(v, comp) {
ckmax(mergedMore[labels[find(v)]], s.more[v]);
ckmin(mergedLess[labels[find(v)]], s.less[v]);
}
int invalidMask[10] = {};
F0R(v, 10) {
F0R(i, cnt) {
if (v <= mergedMore[i] || v >= mergedLess[i]) invalidMask[v] |= (1<<i);
}
dbg(v, invalidMask[v]);
}
// find masks that don't have one be a parent of another
int noParent[1<<cnt];
F0R(i, 1<<cnt) {
bool ok = 1;
F0R(j, cnt) {
if (i&(1<<j))
trav(k, gr[j]) if (i&(1<<k)) ok = 0;
}
noParent[i] = ok;
dbg(i, noParent[i]);
}
int dp[10][1<<cnt];
F0R(v, 10) {
F0R(i, 1<<cnt) {
dp[v][i] = 0;
if (!downwardClosed[i]) continue;
// dbg(v, i, noParent[i], (i&invalidMask[v]));
// loop over submasks of i on previous level
if (v > 0) {
for (int mask = i; mask; mask = (mask-1)&i) {
if (dp[v-1][mask] == 0) continue;
int newMask = (i^mask);
if (newMask&invalidMask[v]) continue;
if (!noParent[newMask]) continue;
if (!downwardClosed[mask]) continue;
dbg(v, i, mask, newMask, dp[v-1][mask]);
dp[v][i] = (dp[v][i]+dp[v-1][mask])%MOD;
}
}
// check 0 separately
if (noParent[i] && (i&invalidMask[v]) == 0) {
dbg(v, i, 0, i, 1);
dp[v][i] = (dp[v][i]+1)%MOD;
}
dbg(v, i, dp[v][i]);
}
}
dbg(i, dp[9][(1<<cnt)-1]);
ret = (1LL*ret*dp[9][(1<<cnt)-1])%MOD;
}
dbg(ret);
return ret;
}
void solve(int k, State s) {
if (k == n) {
// solve this graph
ans = (ans+calc(s))%MOD;
return;
}
string a = p[k].f.f;
string b = p[k].f.s;
int m = sz(a);
dbg("solve", k, a, b);
char op = p[k].s;
if (op == '=') {
State nextState = s;
// then everything is equal
bool ok = 1;
F0R(i, m) {
char x = a[i];
char y = b[i];
if (isDigit(x) && isDigit(y)) {
if (x != y) {
ok = 0;
break;
}
} else if (isDigit(x) && !isDigit(y)) {
// then y-'A' is equal to x
// less than x+1, more than x-1
ckmax(nextState.more[y-'A'], (x-'0'-1));
ckmin(nextState.less[y-'A'], (x-'0'+1));
} else if (!isDigit(x) && isDigit(y)) {
// then x-'A' is equal to y
// less than y+1, more than y-1
ckmax(nextState.more[x-'A'], (y-'0'-1));
ckmin(nextState.less[x-'A'], (y-'0'+1));
} else {
if (nextState.gr[x-'A'][y-'A'] > 1) {
ok = 0;
break;
}
nextState.gr[x-'A'][y-'A'] = 1;
nextState.gr[y-'A'][x-'A'] = 1;
}
}
if (!ok) return;
solve(k+1, nextState);
return;
}
if (op == '<') {
swap(a, b);
}
// otherwise, we need to loop through to find the not equal stuff
F0R(j, m) {
// j is the one where the operator applies
State nextState = s;
bool ok = 1;
F0R(i, j) {
char x = a[i];
char y = b[i];
if (isDigit(x) && isDigit(y)) {
if (x != y) {
ok = 0;
break;
}
} else if (isDigit(x) && !isDigit(y)) {
// then y-'A' is equal to x
// less than x+1, more than x-1
ckmax(nextState.more[y-'A'], (x-'0'-1));
ckmin(nextState.less[y-'A'], (x-'0'+1));
} else if (!isDigit(x) && isDigit(y)) {
// then x-'A' is equal to y
// less than y+1, more than y-1
ckmax(nextState.more[x-'A'], (y-'0'-1));
ckmin(nextState.less[x-'A'], (y-'0'+1));
} else {
if (nextState.gr[x-'A'][y-'A'] > 1) {
ok = 0;
break;
}
nextState.gr[x-'A'][y-'A'] = 1;
nextState.gr[y-'A'][x-'A'] = 1;
}
}
if (!ok) continue;
char x = a[j];
char y = b[j];
dbg(j, x, y);
if (isDigit(x) && isDigit(y)) {
if (!(x > y)) {
ok = 0;
continue;
}
} else if (isDigit(x) && !isDigit(y)) {
ckmin(nextState.less[y-'A'], (x-'0'));
} else if (!isDigit(x) && isDigit(y)) {
ckmax(nextState.more[x-'A'], (y-'0'));
} else {
if (nextState.gr[x-'A'][y-'A'] == 1 || nextState.gr[x-'A'][y-'A'] == 3) {
continue;
}
dbg("setting directed edge", x, y);
nextState.gr[x-'A'][y-'A'] = 2; // >
nextState.gr[y-'A'][x-'A'] = 3; // <
}
solve(k+1, nextState);
}
}
int main() {
ios_base::sync_with_stdio(false);
cin.tie(0); cout.tie(0);
cin >> n;
string constraints[n];
F0R(i, n) cin >> constraints[i];
// split by the operator
string opset = "><=";
F0R(i, n) {
string left = "", right = "";
char op;
int j = 0;
while (opset.find(constraints[i][j]) == -1) {
left += constraints[i][j];
j++;
}
op = constraints[i][j];
j++;
while (j < sz(constraints[i])) {
right += constraints[i][j];
j++;
}
// pad the smaller side with leading zeros
if (sz(left) < sz(right)) {
left = string(sz(right) - sz(left), '0') + left;
} else if (sz(left) > sz(right)) {
right = string(sz(left) - sz(right), '0') + right;
}
dbg(left, op, right);
p[i] = {{left, right}, op};
}
// now solve all cases
State initial;
solve(0, initial);
cout << ans << nl;
return 0;
}
这程序好像有点Bug,我给组数据试试?
详细
Test #1:
score: 100
Accepted
time: 1ms
memory: 3580kb
input:
1 P=NP
output:
766136394
result:
ok single line: '766136394'
Test #2:
score: 0
Accepted
time: 1ms
memory: 3580kb
input:
1 2000CNY>3000USD
output:
0
result:
ok single line: '0'
Test #3:
score: 0
Accepted
time: 1ms
memory: 3580kb
input:
4 AB>CD E<A BC>FF EF>F1
output:
23645065
result:
ok single line: '23645065'
Test #4:
score: 0
Accepted
time: 0ms
memory: 3592kb
input:
2 BC>DD BD<EA
output:
27271695
result:
ok single line: '27271695'
Test #5:
score: 0
Accepted
time: 1ms
memory: 3592kb
input:
3 CE>ED CC>BA BB<AC
output:
426829091
result:
ok single line: '426829091'
Test #6:
score: 0
Accepted
time: 7ms
memory: 3828kb
input:
10 KG<EI EJ>DA EB<IH EB>JG KF<CF JC>FC IC<BJ FI>HH KD>AH AE>GJ
output:
87744507
result:
ok single line: '87744507'
Test #7:
score: 0
Accepted
time: 10ms
memory: 3820kb
input:
10 EK<GM EL<DC DH>IH EF>BL IM<LL EH<JA DJ<AL GL>MB DB>FM AI<HA
output:
665533468
result:
ok single line: '665533468'
Test #8:
score: 0
Accepted
time: 28ms
memory: 3940kb
input:
10 OD<FK FJ>NL NH>KB KM>CA CI>JH CI<AH CE>GI CO<EG FA>HA FA<IJ
output:
878923575
result:
ok single line: '878923575'
Test #9:
score: 0
Accepted
time: 319ms
memory: 3696kb
input:
10 YH>UQ UQ>FD YZ>MK FY<GO YV<QW UV<VJ UZ>EB EQ>LX VP>ZF LZ>TS
output:
867624189
result:
ok single line: '867624189'
Test #10:
score: 0
Accepted
time: 279ms
memory: 3684kb
input:
10 YH<UL UD<FY FK<MU MM<GO GG<QW QJ<VQ VZ<EB EG<LX LZ<ZP ZV<TS
output:
57935948
result:
ok single line: '57935948'
Test #11:
score: 0
Accepted
time: 1ms
memory: 3752kb
input:
6 EDDC>AB5A B<C E9A9B>CACAA DE2>A0D DBCDAC>AED3D5 AAA>BB5
output:
169889581
result:
ok single line: '169889581'
Test #12:
score: 0
Accepted
time: 1ms
memory: 3776kb
input:
9 C<B A>B FDF2<FBDB DB>B4 CF>DA EF4<D1A B8<A5 B3>BF FFA<D5B
output:
0
result:
ok single line: '0'
Test #13:
score: 0
Accepted
time: 6ms
memory: 3596kb
input:
5 SP6<GCT J0RFZ<ZZLUX UDY7<UEVX C1CQ>FXTG SOCT07<MEABU8
output:
603602671
result:
ok single line: '603602671'
Test #14:
score: 0
Accepted
time: 3ms
memory: 3712kb
input:
7 F>M G8F<KC5 F06<E8G H5J<BJE M8CDE<DIGMC AE08>EFI7 DM>CI
output:
821791712
result:
ok single line: '821791712'
Test #15:
score: 0
Accepted
time: 1ms
memory: 3780kb
input:
10 PS1>O9O G76>F8S J<S SB>Y4 WS<VM E<N ZR<CV G8T>XPJ J<A KT<LS
output:
97272892
result:
ok single line: '97272892'
Test #16:
score: 0
Accepted
time: 0ms
memory: 3680kb
input:
4 K1TVV0>TOB4QTH E5U5C9>QGDEGU Q9LW3SK>LWFRP DXUQM=V4N4
output:
467745652
result:
ok single line: '467745652'
Test #17:
score: 0
Accepted
time: 1ms
memory: 3704kb
input:
5 BC5F<AC3F FA4<D48306EDD EFDD>FDABE CF5C<AFDDB FAF<C387
output:
808992671
result:
ok single line: '808992671'
Test #18:
score: 0
Accepted
time: 1ms
memory: 3676kb
input:
1 ABCDEFGHIJKLMNOPQRSTUVWX>BCDEFGHIJKLMNOPQRSTUVWXY
output:
835948861
result:
ok single line: '835948861'
Test #19:
score: 0
Accepted
time: 1ms
memory: 3820kb
input:
3 A=A 00109=109 XX=Z
output:
276262510
result:
ok single line: '276262510'
Test #20:
score: 0
Accepted
time: 0ms
memory: 3588kb
input:
2 ABCDEFGHIJKL=CDEFGHIJKLMN OPQRSTUVWXYZ=RSTUVWXYZOPQ
output:
100000
result:
ok single line: '100000'
Test #21:
score: 0
Accepted
time: 12ms
memory: 3832kb
input:
9 N=A8 TT<QO3G LS>JV TSG>U5F D<A934 FK<HKG O>S1 GT<BBCX SG>S
output:
929594610
result:
ok single line: '929594610'
Test #22:
score: 0
Accepted
time: 0ms
memory: 3800kb
input:
0
output:
673653469
result:
ok single line: '673653469'
Test #23:
score: 0
Accepted
time: 0ms
memory: 3792kb
input:
3 AB<CD AC<BD AD<BC
output:
219041723
result:
ok single line: '219041723'
Extra Test:
score: 0
Extra Test Passed