QOJ.ac
QOJ
ID | Problem | Submitter | Result | Time | Memory | Language | File size | Submit time | Judge time |
---|---|---|---|---|---|---|---|---|---|
#232178 | #7613. Inverse Problem | CharlieVinnie | AC ✓ | 22737ms | 227056kb | C++17 | 3.4kb | 2023-10-29 22:58:48 | 2023-10-29 22:58:50 |
Judging History
answer
#include "bits/stdc++.h"
#ifdef DEBUG
#include "PrettyDebug.hpp"
#else
#define debug(...) [](auto...){}(__VA_ARGS__)
#define debuga(...) [](auto...){}(__VA_ARGS__)
#endif
#define For(i,a,b) for(int i=a;i<=b;i++)
#define Rev(i,a,b) for(int i=a;i>=b;i--)
#define Fin(file) freopen(file,"r",stdin)
#define Fout(file) freopen(file,"w",stdout)
#define assume(expr) ((!!(expr))||(exit((fprintf(stderr,"Assumption Failed: %s on Line %d\n",#expr,__LINE__),-1)),false))
#define range basic_string_view
using namespace std; typedef long long ll;
constexpr int mod=1e9+7,B=5;
inline int qpow(int x,int y){int res=1;while(y){if(y&1)res=1ll*res*x%mod;x=1ll*x*x%mod;y>>=1;}return res;}
bitset<mod> vis; int T,Q[11],val[205],ival[205]; vector<int> lis[205],lis2[205]; vector<pair<int,int>> Ans[11];
void solve(int n){
debug(n);
for(int i=0,t=1;i<=n-2;i++) val[i]=t,t=1ll*t*(n-2-i)%mod;
For(i,0,n-2) ival[i]=qpow(val[i],mod-2);
For(i,0,n-2) lis[i].clear(),lis2[i].clear();
function<void(int,int,int)> dfs1=[&](int d,int tot,int x){
lis[tot].push_back(x);
For(i,d,min(B,n-2-tot)) dfs1(i,tot+i,1ll*x*ival[i]%mod);
};
function<void(int,int,int)> dfs2=[&](int d,int tot,int x){
lis2[tot].push_back(x);
For(i,d,n-2-tot) dfs2(i,tot+i,1ll*x*val[i]%mod);
};
dfs1(1,0,1); dfs2(B+1,0,1);
For(id,1,T) if(Ans[id].empty()){
int o=1ll*qpow(n,mod-2)%mod*qpow(n-1,mod-2)%mod*Q[id]%mod;
For(s,0,n-2){
for(int x:lis[s]) vis.set(1ll*o*x%mod);
int X=0;
for(int x:lis2[n-2-s]) if(vis.test(x)) { X=x; break; }
for(int x:lis[s]) vis.reset(1ll*o*x%mod);
if(!X) continue;
vector<int> ans;
function<bool(int,int,int)> dfs3=[&](int d,int tot,int x){
if(tot==s&&x==X) return true;
For(i,d,min(B,s-tot)){
ans.push_back(i); if(dfs3(i,tot+i,1ll*x*ival[i]%mod)) return true; else ans.pop_back();
}
return false;
};
function<bool(int,int,int)> dfs4=[&](int d,int tot,int x){
if(tot==n-2-s&&x==X) return true;
For(i,d,n-2-s-tot){
ans.push_back(i); if(dfs4(i,tot+i,1ll*x*val[i]%mod)) return true; else ans.pop_back();
}
return false;
};
assert(dfs3(1,0,o)); assert(dfs4(B+1,0,1));
ans.resize(n); debug(ans);
for(int i=1,p=2;i<=n;i++){
For(_,1,ans[i-1]+(i==1)) Ans[id].emplace_back(i,p++);
}
break;
}
}
}
signed main(){
atexit([](){cerr<<"Time = "<<clock()<<" ms"<<endl;});
// int T; cin>>T; while(T--) { int o; cin>>o; if(o==1) { cout<<1<<'\n'; continue; }; int n=2; while(!solve(o,n)) n++; }
cin>>T; For(i,1,T) { cin>>Q[i]; if(Q[i]==1) Ans[i].emplace_back(0,0); }
for(int n=2;;n++){
solve(n); int ok=1; For(i,1,T) if(Ans[i].empty()) { ok=0; break; } if(ok) break;
}
For(i,1,T){
if(Ans[i][0].first==0) cout<<1<<'\n';
else{
cout<<Ans[i].size()+1<<'\n'; for(auto [x,y]:Ans[i]) cout<<x<<' '<<y<<'\n';
}
}
return 0;
}
// START TYPING IF YOU DON'T KNOW WHAT TO DO
// STOP TYPING IF YOU DON'T KNOW WHAT YOU'RE DOING
// CONTINUE, NON-STOPPING, FOR CHARLIEVINNIE
// Started Coding On: October 28 Sat, 07 : 43 : 27
Details
Tip: Click on the bar to expand more detailed information
Test #1:
score: 100
Accepted
time: 4ms
memory: 116212kb
input:
4 2 360 1 509949433
output:
2 1 2 5 1 2 1 3 2 4 2 5 1 10 1 2 1 3 2 4 3 5 4 6 5 7 6 8 7 9 8 10
result:
ok OK (4 test cases)
Test #2:
score: 0
Accepted
time: 10199ms
memory: 227056kb
input:
9 185396120 468170792 837583517 696626231 338497514 762842660 800028852 928391161 733524004
output:
14 1 2 1 3 2 4 3 5 4 6 5 7 6 8 7 9 8 10 9 11 9 12 10 13 10 14 122 1 2 1 3 2 4 3 5 4 6 5 7 6 8 6 9 7 10 7 11 8 12 8 13 9 14 9 15 9 16 10 17 10 18 10 19 11 20 11 21 11 22 12 23 12 24 12 25 12 26 13 27 13 28 13 29 13 30 14 31 14 32 14 33 14 34 14 35 15 36 15 37 15 38 15 39 15 40 15 41 15 42 16 43 16 44...
result:
ok OK (9 test cases)
Test #3:
score: 0
Accepted
time: 22737ms
memory: 221968kb
input:
10 338497514 733524004 447182954 415993605 453460670 50499055 648088551 986982752 907925397 315315230
output:
124 1 2 1 3 2 4 3 5 4 6 5 7 6 8 7 9 8 10 9 11 10 12 11 13 11 14 12 15 12 16 12 17 13 18 13 19 13 20 13 21 13 22 14 23 14 24 14 25 14 26 14 27 14 28 15 29 15 30 15 31 15 32 15 33 15 34 15 35 16 36 16 37 16 38 16 39 16 40 16 41 16 42 17 43 17 44 17 45 17 46 17 47 17 48 17 49 18 50 18 51 18 52 18 53 18...
result:
ok OK (10 test cases)
Test #4:
score: 0
Accepted
time: 2337ms
memory: 143000kb
input:
10 1 2 3 4 5 6 7 8 9 10
output:
1 2 1 2 102 1 2 1 3 2 4 3 5 4 6 4 7 5 8 5 9 6 10 6 11 7 12 7 13 8 14 8 15 9 16 9 17 9 18 9 19 9 20 9 21 10 22 10 23 10 24 10 25 10 26 10 27 10 28 10 29 11 30 11 31 11 32 11 33 11 34 11 35 11 36 11 37 11 38 12 39 12 40 12 41 12 42 12 43 12 44 12 45 12 46 12 47 12 48 13 49 13 50 13 51 13 52 13 53 13 5...
result:
ok OK (10 test cases)
Test #5:
score: 0
Accepted
time: 611ms
memory: 131504kb
input:
10 269199917 392009324 753889928 751355133 472639410 132096559 331333826 40414701 72847302 475706026
output:
55 1 2 1 3 1 4 1 5 1 6 1 7 1 8 1 9 1 10 1 11 1 12 1 13 1 14 1 15 1 16 1 17 1 18 1 19 1 20 1 21 1 22 1 23 1 24 1 25 1 26 1 27 1 28 1 29 1 30 1 31 1 32 1 33 1 34 1 35 1 36 1 37 1 38 1 39 1 40 1 41 1 42 1 43 1 44 1 45 1 46 1 47 1 48 1 49 1 50 1 51 1 52 1 53 1 54 1 55 84 1 2 1 3 2 4 3 5 4 6 5 7 6 8 7 9 ...
result:
ok OK (10 test cases)
Extra Test:
score: 0
Extra Test Passed