QOJ.ac
QOJ
ID | 题目 | 提交者 | 结果 | 用时 | 内存 | 语言 | 文件大小 | 提交时间 | 测评时间 |
---|---|---|---|---|---|---|---|---|---|
#88615 | #2574. Fancy Arrays | s8194272 | WA | 135ms | 25136kb | C++14 | 3.3kb | 2023-03-16 19:14:22 | 2023-03-16 19:14:23 |
Judging History
answer
#include<iostream>
#include<algorithm>
#include<cstring>
#include<cstdio>
#include<bitset>
#include<random>
#include<cmath>
#include<ctime>
#include<queue>
#include<map>
#include<set>
#include<assert.h>
#define int long long
#define fi first
#define se second
#define max Max
#define min Min
#define abs Abs
#define lc (x<<1)
#define rc (x<<1|1)
#define mid ((l+r)>>1)
#define pb(x) push_back(x)
#define lowbit(x) ((x)&(-(x)))
#define fan(x) ((((x)-1)^1)+1)
#define mp(x,y) make_pair(x,y)
#define clr(f,n) memset(f,0,sizeof(int)*(n))
#define cpy(f,g,n) memcpy(f,g,sizeof(int)*(n))
#define SZ(x) ((int)(x.size()))
#define INF 0x3f3f3f3f
using namespace std;
inline int read()
{
int ans=0,f=1;
char c=getchar();
while(c>'9'||c<'0'){if(c=='-')f=-1;c=getchar();}
while(c>='0'&&c<='9'){ans=(ans<<1)+(ans<<3)+c-'0';c=getchar();}
return ans*f;
}
inline void write(int x)
{
if(x<0) putchar('-'),x=-x;
if(x/10) write(x/10);
putchar((char)(x%10)+'0');
}
template<typename T>inline T Abs(T a){return a>0?a:-a;};
template<typename T,typename TT>inline T Min(T a,TT b){return a<b?a:b;}
template<typename T,typename TT> inline T Max(T a,TT b){return a<b?b:a;}
const int N=205,mod=998244353;
int n,m,q,k,a[N],b[N],c[N],C[N][N],cnt;
struct Matrix
{
int n,m,c[N][N];
Matrix()
{
memset(c,0,sizeof(c));
}
Matrix(int a,int b)
{
n=a;m=b;
memset(c,0,sizeof(c));
}
Matrix operator * (const Matrix p)
{
Matrix res(n,p.n);
assert(m==p.n);
for(int i=1;i<=n;++i)
for(int j=1;j<=p.m;++j)
for(int k=1;k<=p.n;++k)
res.c[i][j]=(res.c[i][j]+c[i][k]*p.c[k][j]%mod)%mod;
return res;
}
}pw[60],ans;
vector<int> now;
vector<vector<int> > state;
void dfs(int u)
{
if(u==cnt+1)
{
state.push_back(now);++k;
return;
}
for(int i=0;i<=b[u];++i)
{
now.push_back(i);
dfs(u+1);
now.pop_back();
}
}
inline int q_pow(int a,int b)
{
int c=1;
while(b)
{
if(b&1) c=a*c%mod;
a=a*a%mod;b>>=1;
}
return c;
}
signed main()
{
m=read();q=read();
for(int i=0;i<=200;++i)
for(int j=0;j<=i;++j)
if(j==0) C[i][j]=1;
else C[i][j]=(C[i-1][j-1]+C[i-1][j])%mod;
for(int i=2;i*i<=m;++i)
if(m%i==0)
{
++k;
while(m%i==0)
m/=i,c[k]++;
}
if(m!=1) c[++k]=1;
sort(c+1,c+1+k);
for(int i=1;i<=k;++i)
{
if(c[i]!=c[i-1])
a[++cnt]=c[i],b[cnt]=1;
else b[cnt]++;
}
k=0;dfs(1);
pw[0].n=pw[0].m=k;
for(int i=1;i<=k;++i)
for(int j=1;j<=k;++j)
{
int g[2]={1,0};
for(int p=1;p<=cnt;++p)
{
int all=C[b[p]][state[j-1][p-1]]*q_pow(a[p],state[j-1][p-1])%mod;
int w0=(b[p]-state[i-1][p-1]<state[j-1][p-1]?0:C[b[p]-state[i-1][p-1]][state[j-1][p-1]]*q_pow(a[p],state[j-1][p-1])%mod),w1=(all-w0+mod)%mod;
g[1]=(g[1]*all+g[0]*w1)%mod;
g[0]=g[0]*w0%mod;
}
pw[0].c[j][i]=g[1];
}
for(int i=1;i<=59;++i)
pw[i]=pw[i-1]*pw[i-1];
while(q--)
{
n=read()-1;
ans.n=k;ans.m=1;
for(int i=1;i<=k;++i)
{
int tmp=1;
for(int j=1;j<=cnt;++j)
tmp=tmp*C[b[j]][state[i-1][j-1]]%mod*q_pow(a[j],state[i-1][j-1])%mod;
ans.c[i][1]=tmp;
}
for(int i=0;i<=59;++i)
if((n>>i)&1)
ans=pw[i]*ans;
int Ans=0;
for(int i=1;i<=k;++i)
Ans=(Ans+ans.c[i][1])%mod;
write(Ans);puts("");
}
return 0;
}
詳細信息
Test #1:
score: 100
Accepted
time: 17ms
memory: 24992kb
input:
12 3 1 2 3
output:
6 21 91
result:
ok 3 number(s): "6 21 91"
Test #2:
score: 0
Accepted
time: 135ms
memory: 24976kb
input:
1 150 471816347971198367 934144370769132530 85747619384378846 928941512582005725 154937870030720168 947932149793416512 27783441557851811 522085897018258944 254251197759739965 280173028039582607 323577718378116194 390211126917894813 349211961997885462 482844442408522462 582732208453073301 94800734555...
output:
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
result:
ok 150 numbers
Test #3:
score: 0
Accepted
time: 133ms
memory: 24984kb
input:
2 150 879653409269605014 957081824205994700 92943925332284309 70508831927780168 72367833784810922 57052500883916366 260855517197770739 493364569696106472 261906268272035425 712282959082227662 35005533487670014 740269757357303611 472541044721679500 231355986572948422 563516773952248704 38987628675191...
output:
1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
result:
ok 150 numbers
Test #4:
score: -100
Wrong Answer
time: 122ms
memory: 25136kb
input:
4 150 833174642454220423 721913650877167279 111257970647375842 922819627392160450 408011919008881312 938552585499192014 401181394137854787 154596954164557809 43303362814617574 450360165684736834 713407776281798115 265067947883317301 820681723927726574 17493726254665319 431343457571478167 51814600647...
output:
951575782 619990425 288768710 441870637 893677154 985627908 157875830 986732592 174416053 114900859 254948884 159412823 561233616 53041748 706692641 47502474 727130388 310080803 419736807 113411171 446032540 824422898 963360110 255866093 260373474 918260634 418338816 773394191 124312932 279569235 66...
result:
wrong answer 1st numbers differ - expected: '468840309', found: '951575782'