QOJ.ac

QOJ

ID题目提交者结果用时内存语言文件大小提交时间测评时间
#105106#4793. QnpDeterminantWA 1071ms7140kbC++141.8kb2023-05-12 23:18:212023-05-12 23:18:22

Judging History

你现在查看的是最新测评结果

  • [2023-08-10 23:21:45]
  • System Update: QOJ starts to keep a history of the judgings of all the submissions.
  • [2023-05-12 23:18:22]
  • 评测
  • 测评结果:WA
  • 用时:1071ms
  • 内存:7140kb
  • [2023-05-12 23:18:21]
  • 提交

answer

#include<bits/stdc++.h>
using namespace std;typedef long long ll;
const int N=70007;const ll p=1445141919821ll,p2=1e9+7;
int a[11],st,minn,fl,cnt;ll pow0[N],pow1[N],ans,m,inv[N],f[N],f_[N];double s[N];
void sol(){
	ans=st=0;for(int i=0;i<10;++i)scanf("%d",&a[i]),st+=a[i];
	scanf("%lld",&m);double qu=0;ll wz=1;
	for(int i=0;i<10;++i)qu-=s[a[i]],wz=(__int128)wz*f_[a[i]]%p;
	for(int i=0;i<10;++i)if(a[i]){
		if(qu+s[st-a[i]]+s[a[i]]>12.1||m<=(__int128)wz*f[st-a[i]]*f[a[i]]%p){
			qu+=s[a[i]];ans=(ans*pow0[a[i]]+pow1[a[i]]*i)%p2;
			wz=(__int128)wz*f[a[i]]%p;st-=a[i];a[i]=0;
		}
		else{minn=i;break;}
	}
	while(st){
		for(int i=minn;i<10;++i)if(a[i]){
			if(qu+s[st-1]+s[a[i]]-s[a[i]-1]>12.1){fl=i;break;}
			ll q=(__int128)wz*f[st-1]*a[i]%p;if(m<=q){fl=i;break;}m-=q;
		}
		if(fl==minn&&a[minn]>8&&(qu+s[st-8]+s[a[minn]]-s[a[minn]-8]>12.1||
		m<=(__int128)wz*f[st-8]*f[a[minn]]%p*f_[a[minn]-8]%p)){
			int l=0,r=a[minn],mid;
			while(l<r){
				mid=(l+r+1)/2;if(qu+s[st-mid]+s[a[minn]]-s[a[minn]-mid]>12.1||
				m<=(__int128)wz*f[st-mid]*f[a[minn]]%p*f_[a[minn]-mid]%p)l=mid;
				else r=mid-1;
			}
			qu+=s[a[minn]]-s[a[minn]-l];wz=(__int128)wz*f[a[minn]]*f_[a[minn]-l]%p;
			ans=(ans*pow0[l]+pow1[l]*minn)%p2;st-=l;a[minn]-=l;
		}
		else{ans=(ans*10+fl)%p2;qu+=s[a[fl]]-s[a[fl]-1];wz=(__int128)wz*a[fl]%p;st--;a[fl]--;}
		if(!a[minn]){for(int i=0;i<10;++i)if(a[i]){minn=i;break;}}
	}
	printf("%lld\n",ans);return;
}
int main(){
	inv[1]=1;for(int i=2;i<N;++i)inv[i]=(__int128)(p-p/i)*inv[p%i]%p;
	f[0]=f_[0]=pow0[0]=1;for(int i=1;i<N;++i){
		f[i]=(__int128)f[i-1]*i%p,f_[i]=(__int128)f_[i-1]*inv[i]%p;
		s[i]=s[i-1]+log10(i);pow1[i]=(pow1[i-1]*10+1)%p2;
		pow0[i]=pow0[i-1]*10%p2;
	}
	int t;scanf("%d",&t);t=min(t,2500);while(t--)sol();return 0;
}

详细

Test #1:

score: 100
Accepted
time: 3ms
memory: 7112kb

input:

6
1 1 0 0 0 0 0 0 0 0 1
1 1 0 0 0 0 0 0 0 0 2
1 1 1 0 0 0 0 0 0 0 1
1 1 1 0 0 0 0 0 0 0 2
1 1 1 0 0 0 0 0 0 0 5
1 2 0 0 0 0 0 0 0 0 2

output:

1
10
12
21
201
101

result:

ok 6 numbers

Test #2:

score: -100
Wrong Answer
time: 1071ms
memory: 7140kb

input:

5000
7142 6835 7022 6981 6929 7008 7038 7015 6885 7145 659213485437
7015 7033 6963 7136 7053 7072 6847 6923 6953 7005 82053183749
7013 7003 6969 7000 7011 7137 7030 6823 7021 6993 817812793310
6893 7008 6963 7086 7012 6922 7128 7094 7028 6866 143084249211
7020 7021 6997 6961 7017 7032 7013 7028 6987...

output:

190295171
625739245
198000759
198944878
575159465
972833519
64095209
49623132
936823516
4168747
3116218
789303150
240806129
772063597
812079040
453627969
543807021
604471487
143706250
632200817
838349641
517706711
785181565
476902075
327278139
613377939
161377153
968347598
454781889
685567117
331121...

result:

wrong answer Answer contains longer sequence [length = 5000], but output contains 2500 elements