QOJ.ac
QOJ
ID | Problem | Submitter | Result | Time | Memory | Language | File size | Submit time | Judge time |
---|---|---|---|---|---|---|---|---|---|
#480720 | #7613. Inverse Problem | ucup-team052 | TL | 0ms | 0kb | C++14 | 2.7kb | 2024-07-16 17:47:15 | 2024-07-16 17:47:16 |
answer
#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
const int md = 1e9 + 7;
const int mod = 1e9 + 7;
inline int add(int x, int y) {
if (x + y >= md) return x + y - md;
return x + y;
}
inline int sub(int x, int y) {
if (x < y) return x - y + md;
return x - y;
}
inline int mul(int x, int y) {
return 1ull * x * y % md;
}
inline int mul(int x,int y,int z) {return mul(mul(x,y),z);}
inline int qpow(int x, int y) {
int ans = 1;
while (y) {
if (y & 1) ans = mul(ans, x);
y >>= 1; x = mul(x, x);
}
return ans;
}
int Inv(int x) {return qpow(x,mod-2);}
struct Node
{
int siz;
vector<pair<int,int>> eg;
Node(int s=0)
{
assert(s<=1);
siz=s;
}
};
const int L=200;
int n,dv[1005][1005];
map<int,vector<int>> f[L+1];
void init()
{
for(int i=0;i<=n;i++) f[i].clear();
f[0][1]=vector<int>(0);
for(int i=1;i<=2*n/3+2;i++)
{
for(int j=1;j<=i;j++)
{
int adw=dv[n-2][j];
for(auto [w,v]:f[i-j])
{
int gw=mul(adw,w);
if(f[i].find(gw)==f[i].end())
{
f[i][gw]=f[i-j][w];
f[i][gw].push_back(j);
}
}
}
// printf("%d %d\n",i,(int)f[i].size());
}
}
int Q,q[15],ok[15],sok;
vector<int> ans[15];
int ansn[15];
int main() {
#ifdef xay5421
freopen("a.in","r",stdin);
#endif
for(int i=1;i<=1000;i++)
{
dv[i][0]=1;
for(int j=1;j<=1000;j++) dv[i][j]=mul(dv[i][j-1],i-j+1);
}
cin>>Q;
for(int i=1;i<=Q;i++)
{
cin>>q[i];
if(q[i]==1) ok[i]=1,sok++;
if(q[i]==2)
{
ok[i]=1,sok++;
}
}
for(n=3;n<=L&&sok<Q;n++)
{
fprintf(stderr,"n = %d\n",n);
init();
int adw=mul(n,n-1);
for(int qid=1;qid<=Q;qid++)
{
// q[qid] sum=n-2
for(int big=1;big<=n-2&&!ok[qid];big++)
{
int cw=dv[n-2][big];
for(int L=0;L<=big&&!ok[qid];L++)
{
int R=n-2-big-L;
if(!(R>=0)) continue;
for(auto [w,v]:f[L])
{
int way=mul(adw,w,cw);
int ned=mul(Inv(way),q[qid]);
if(f[R].find(ned)!=f[R].end())
{
ans[qid]=f[R][ned];
for(int i:v) ans[qid].push_back(i);
ans[qid].push_back(big);
ok[qid]=1;
ansn[qid]=n;
break;
}
}
}
}
}
}
for(int i=1;i<=Q;i++)
{
if(q[i]==1) printf("1\n");
else if(q[i]==2) printf("2\n1 2\n");
else
{
printf("%d\n",ansn[i]);
while(ans[i].size()<ansn[i]) ans[i].push_back(0);
// for(int j=0;j<(int)ans[i].size();j++) ans[i][j]++;
for(int j=0;j<ansn[i]-1;j++)
{
for(int c=0;c<(int)ans[i].size();c++) for(int d=0;d<(int)ans[i].size();d++)
{
if(ans[i][c]==0&&ans[i][d]>0)
{
ans[i][c]--,ans[i][d]--;
printf("%d %d\n",c+1,d+1);
}
}
}
}
}
return 0;
}
Details
Tip: Click on the bar to expand more detailed information
Test #1:
score: 0
Time Limit Exceeded
input:
4 2 360 1 509949433