QOJ.ac

QOJ

IDProblemSubmitterResultTimeMemoryLanguageFile sizeSubmit timeJudge time
#92317#618. 多项式乘法cc000#100 ✓920ms82564kbC++141.6kb2023-03-30 15:48:452023-03-30 15:48:48

Judging History

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

  • [2023-08-10 23:21:45]
  • System Update: QOJ starts to keep a history of the judgings of all the submissions.
  • [2023-03-30 15:48:48]
  • 评测
  • 测评结果:100
  • 用时:920ms
  • 内存:82564kb
  • [2023-03-30 15:48:45]
  • 提交

answer

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
const int maxn=5000020;
namespace POLY
{
	const int mod=998244353,g=3,ginv=332748118;
	ll ksm(ll x,ll y)
	{
		ll ret=1;
		while(y)
		{
			if(y&1) ret=ret*x%mod;
			x=x*x%mod; y>>=1;
		}
		return ret;
	}
	int pos[maxn];
	int pre(int n,int m)
	{
		int lim=1;
		while((1<<lim)<=n+m) lim++;
		for(int i=0;i<(1<<lim);i++) pos[i]=(pos[i>>1]>>1)|((i&1)<<(lim-1));
		return lim;		
	}
	void NTT(ll *f,int n,int opt)
	{
		for(int i=0;i<n;i++) if(i<pos[i]) swap(f[i],f[pos[i]]);
		for(int i=1;i<n;i<<=1)
		{
			ll step=ksm(opt>0?g:ginv,(mod-1)/(2*i));
			for(int j=0;j<n;j+=i+i)
			{
				ll cur=1;
				for(int k=j;k<j+i;k++)
				{
					ll fk=f[k];ll hk=f[k+i]*cur%mod;
					f[k]=(fk+hk)%mod; f[k+i]=(fk-hk+mod)%mod;
					cur=cur*step%mod;
				}
			}
		}
		if(opt<0)
		{
			ll ninv=ksm(n,mod-2);
			for(int i=0;i<n;i++) f[i]=f[i]*ninv%mod;
		}
	}
	void MUL(ll *s,ll *F,ll *G,int n,int m)
	{
		static ll A[maxn],B[maxn];
		int lim=pre(n,m);
		for(int i=0;i<=n;i++) A[i]=F[i];
		for(int i=0;i<=m;i++) B[i]=G[i];
		for(int i=n+1;i<(1<<lim);i++) A[i]=0;
		for(int i=m+1;i<(1<<lim);i++) B[i]=0;
		NTT(A,1<<lim,1); NTT(B,1<<lim,1);
		for(int i=0;i<(1<<lim);i++) s[i]=A[i]*B[i]%mod;
		NTT(s,1<<lim,-1);
	}
}
ll f[maxn],g[maxn],h[maxn];
int n,m;
int main()
{
	//freopen("C.in","r",stdin);
	scanf("%d%d",&n,&m);
	for(int i=0;i<=n;i++) scanf("%lld",&f[i]);
	for(int i=0;i<=m;i++) scanf("%lld",&g[i]);
	POLY::MUL(h,f,g,n,m);
	for(int i=0;i<=n+m;i++)
		printf("%lld ",h[i]);
}

Details

Tip: Click on the bar to expand more detailed information

Test #1:

score: 20
Accepted
time: 6ms
memory: 13876kb

input:

96 96
600395131 184265451 942971382 534262851 830366882 542271170 294355449 501371170 797809599 964826049 276651245 375755165 662619442 941920605 328216963 507795473 460271147 874920847 818231910 156789488 590591583 732194508 793983630 93566697 836155866 305319153 432040686 621119061 835023373 57138...

output:

683858396 5532883 499734624 910262414 221004044 924081841 392466229 64190174 260661815 939986106 283456690 260629512 990528995 704246427 991946815 236857583 903415172 900324859 938555797 225258152 874945420 516870315 74759441 769850097 353889928 300397164 63689540 115003940 872945378 407694641 91843...

result:

ok 193 numbers

Test #2:

score: 20
Accepted
time: 6ms
memory: 13976kb

input:

4992 4994
471194390 313639917 705341086 119536186 430124603 244978845 185588341 13731324 707132801 88167972 927324568 846658454 523684029 5133605 767200778 502476309 539772401 778154025 266136872 183516351 260704325 49303370 475056182 928574546 740424153 277920667 708439854 746983628 839491869 53579...

output:

700935456 677302967 772159864 479386810 109686665 263919131 29567167 960045078 636326916 585682137 409426717 14510019 441964472 92801447 551536199 216995135 59736203 790078879 55883568 796138076 265361608 66124731 150347029 93682849 205256362 672081205 86396898 573029352 541084997 293480941 90518071...

result:

ok 9987 numbers

Test #3:

score: 20
Accepted
time: 40ms
memory: 16384kb

input:

29995 29992
417238081 53580806 733071257 224121793 786137422 127072245 129083351 988357079 246853229 150935424 596994106 975796660 838029970 619117898 328485797 948485083 574261409 79312345 596346086 489787404 929520168 515647000 211731260 50868568 811515357 428215135 498099163 644485329 802849075 3...

output:

115270920 49832720 758693293 745763567 322999821 510579248 697424729 850661201 678364508 817667211 668544763 136619207 562899653 692811546 351397117 768369036 573254435 891143982 717302438 707939747 41743610 540709722 240732780 931265491 38731999 642520590 630812534 632188732 342954490 225414102 836...

result:

ok 59988 numbers

Test #4:

score: 20
Accepted
time: 101ms
memory: 23372kb

input:

100000 99993
812398607 947396010 797321381 883949986 56052416 586258761 193247973 611124334 773505112 142179482 565466227 140875825 79890768 893500101 553768089 648879319 480419657 915530184 799329430 494818755 793895824 851865180 459534006 259473419 610037701 472768430 868914058 887444584 588850309...

output:

821875273 646409297 701893040 744951544 891720486 338002304 134405948 686576985 653633849 704180950 763960458 160339533 773107048 630019221 467173934 675237413 824356289 394352126 870024535 473719536 246319541 372709664 656104889 677100818 890131281 374587639 160832628 144239351 450760970 646488586 ...

result:

ok 199994 numbers

Test #5:

score: 20
Accepted
time: 920ms
memory: 82564kb

input:

999993 999994
388529697 811245378 165909114 295553883 667981275 78502012 400874009 139394758 249494489 4636487 997712665 259780805 431039016 716944209 709300152 356513646 823185021 699568300 650937921 859190797 899514799 785648601 933470757 627225124 349752104 471458923 456404256 48134357 315599086 ...

output:

199012842 735467570 660520906 870291510 102406003 509914017 591503608 692425397 149848591 232605296 411728228 285507919 90090498 682749099 507720817 425946949 937188332 619041823 738654334 153862895 272311969 793838225 260785140 350903642 151151058 631242104 304026658 123734332 23714740 438936743 77...

result:

ok 1999988 numbers