QOJ.ac

QOJ

IDProblemSubmitterResultTimeMemoryLanguageFile sizeSubmit timeJudge time
#586091#9380. MatchUESTC_Snow_HalationAC ✓303ms6672kbC++143.0kb2024-09-24 02:10:502024-09-24 02:10:51

Judging History

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

  • [2024-09-24 02:10:51]
  • 评测
  • 测评结果:AC
  • 用时:303ms
  • 内存:6672kb
  • [2024-09-24 02:10:50]
  • 提交

answer

#include <algorithm>
#include <iostream>
#include <cstring>
#include <cstdio>
#include <cmath>
#define QWQ cout<<"QwQ"<<endl;
#define ll long long
#define poly vector<ll>
#include <vector>
#include <queue>
#include <stack>
#include <map>
using namespace std;
const ll N=101010;
const ll qwq=303030;
const ll inf=0x3f3f3f3f3f3f3f3f;
const ll p=998244353;

inline ll read() {
	ll sum = 0, f = 1; char c = getchar();
	while(c<'0' || c>'9') { if(c=='-') f = -1; c = getchar(); }
	while(c>='0'&&c<='9') { sum = sum * 10 + c - '0'; c = getchar(); }
	return sum * f;
}

ll n,K;
ll a[N],b[N];
poly ans;

ll f[N],ni[N];

inline ll ksm(ll aa,ll bb) {
    ll sum = 1;
    while(bb) {
        if(bb&1) sum = sum * aa %p;
        bb >>= 1; aa = aa * aa %p;
    }
    return sum;
}

inline ll CC(ll aa,ll bb) { return f[aa] * ni[bb] %p * ni[aa-bb] %p; }

void qiu() {
    f[0] = ni[0] = 1;
    for(ll i=1;i<=N-10;i++) f[i] = f[i-1] * i %p;
    ni[N-10] = ksm(f[N-10],p-2);
    for(ll i=N-11;i>=1;i--) ni[i] = ni[i+1] * (i+1) %p;
}

inline ll pipei(ll aa,ll bb,ll cc) {
	return CC(aa,cc) * CC(bb,cc) %p * f[cc] %p;
}

void outing(poly A) {
	for(ll v : A) cout<<v<<" ";
	cout<<endl;
}

poly DFS(ll la,ll ra,ll lb,ll rb,ll k) {
	// cout<<la<<","<<ra<<" "<<lb<<","<<rb<<" "<<k<<endl;
	// system("pause");
	poly res;
	if(la>ra || lb>rb) {
		res.resize(1);
		res[0] = 1;
		// outing(res);
		return res;
	}

	if(k==-1) {
		res.resize( min(ra-la+1, rb-lb+1)+1 );
		for(ll i=0;i<res.size();i++) res[i] = pipei(ra-la+1,rb-lb+1,i);
			// outing(res);
		return res;
	}

	ll l = la, r = ra, mid, mida = la-1, midb = lb-1;
	while(l<=r) {
		mid = l+r >> 1;
		if((a[mid]>>k)&1) r = mid-1;
		else mida = mid, l = mid+1;
	}
	l = lb, r = rb;
	while(l<=r) {
		mid = l+r >> 1;
		if((b[mid]>>k)&1) r = mid-1;
		else midb = mid, l = mid+1;
	}

	if((K>>k)&1) {
		poly A = DFS(la, mida, midb+1, rb, k-1);
		poly B = DFS(mida+1, ra, lb, midb, k-1);
		res.resize( A.size()+B.size()-1 );
		for(ll i=0;i<A.size();i++)
			for(ll j=0;j<B.size();j++)
				(res[i+j] += A[i] * B[j] %p) %= p;
			// outing(res);
		return res;
	}
	else {
		poly A = DFS(la, mida, lb, midb, k-1);
		poly B = DFS(mida+1, ra, midb+1, rb, k-1);
		res.resize( min(ra-la+1, rb-lb+1) + 1 );
		for(ll i=0;i<A.size();i++) {
			for(ll j=0;j<B.size();j++) {
				ll A0 = mida-la+1-i;
				ll B0 = midb-lb+1-i;
				ll A1 = ra-mida-j;
				ll B1 = rb-midb-j;
				for(ll x=0;x<=min(A0,B1);x++)
					for(ll y=0;y<=min(A1,B0);y++)
						(res[i+j+x+y] += A[i] * B[j] %p * pipei(A0,B1,x) %p * pipei(A1,B0,y) %p) %= p;
			}
		}
		// outing(res);
		return res;
	}
}

int main() {
	// freopen("data.in", "r", stdin);
	// freopen("K.in", "w", stdout);
	qiu();
	n = read(); K = read();
	for(ll i=1;i<=n;i++) a[i] = read();
	for(ll i=1;i<=n;i++) b[i] = read();
	sort(a+1,a+n+1);
	sort(b+1,b+n+1);
	ans = DFS(1,n,1,n,61);
	for(ll i=1;i<=n;i++) cout<<ans[i]<<"\n";
	return 0;
}

/*

9 5
1 2 2 4 5 7 8 9 10
1 2 3 4 5 8 8 8 9

*/

这程序好像有点Bug,我给组数据试试?

Details

Tip: Click on the bar to expand more detailed information

Test #1:

score: 100
Accepted
time: 0ms
memory: 6672kb

input:

9 5
1 7 2 8 4 9 2 5 10
1 3 2 4 5 8 8 8 9

output:

51
1034
10768
62195
200965
348924
294444
97344
7200

result:

ok 9 lines

Test #2:

score: 0
Accepted
time: 303ms
memory: 5156kb

input:

200 569102225177443347
1103663383682482176 1103381908705771520 1099441259031822336 1098878309078401024 1089871109823660032 1080863910568919040 1074108511127863296 1071856711314178048 1069041961547071488 1068479011593650176 1067353111686807552 1062849512059437056 1049338713177325568 10448351135499550...

output:

20506
208107723
47878304
53020813
972282728
933586852
658157196
670189811
957980024
366179738
217980591
967482558
833450149
987731802
260904367
5263881
600332344
906061351
658256294
93700706
421323952
178075016
219871690
986880524
848776106
191185484
641917326
576497440
908609746
349728876
606714342...

result:

ok 200 lines

Test #3:

score: 0
Accepted
time: 18ms
memory: 5860kb

input:

200 1098723432450995412
972777519512027136 936748722493063168 864691128455135232 860187528827764736 857935729014079488 856246879153815552 855683929200394240 851180329573023744 848928529759338496 847802629852495872 847310048643252224 847239679899074560 846676729945653248 828662331436171264 8280993814...

output:

8320
34035512
428014253
916072411
696504298
748377440
32424677
188402417
233587075
130639672
759476380
546285625
187068736
159002787
131866334
381530906
133126344
373727612
923311054
725293681
718162548
39511135
322320638
44709653
156884882
285926751
787179409
282967016
153092344
615721347
950855850...

result:

ok 200 lines

Test #4:

score: 0
Accepted
time: 0ms
memory: 6608kb

input:

200 1
1 1 1 1 1 1 1 1 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 ...

output:

3072
4412928
940025853
665370743
752672705
581549490
41990996
541401698
170451802
26584141
220983766
844620126
64506869
137621326
418866920
351049174
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 200 lines

Test #5:

score: 0
Accepted
time: 3ms
memory: 6528kb

input:

200 0
1 1 1 1 1 1 1 1 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 ...

output:

40000
792020000
319908096
286630939
515085548
95626140
472749362
800008042
594973486
139901984
847967250
125664877
254232250
419486325
204481923
979660340
41873095
540955890
585129047
205968776
196086667
674391130
737319172
669580034
973843015
207125612
331414909
119844754
24879477
990932807
8476927...

result:

ok 200 lines

Extra Test:

score: 0
Extra Test Passed