結果

問題 No.2206 Popcount Sum 2
コンテスト
ユーザー vjudge1
提出日時 2026-06-08 22:56:28
言語 C++17
(gcc 15.2.0 + boost 1.89.0)
コンパイル:
g++-15 -O2 -lm -std=c++17 -Wuninitialized -DONLINE_JUDGE -o a.out _filename_
実行:
./a.out
結果
AC  
実行時間 427 ms / 4,000 ms
コード長 1,329 bytes
記録
記録タグの例:
初AC ショートコード 純ショートコード 純主流ショートコード 最速実行時間
コンパイル時間 1,551 ms
コンパイル使用メモリ 216,240 KB
実行使用メモリ 14,592 KB
最終ジャッジ日時 2026-06-08 22:56:39
合計ジャッジ時間 10,228 ms
ジャッジサーバーID
（参考情報） 		judge3_0 / judge1_0
純コード判定待ち
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 1
other AC * 18
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

#include<bits/stdc++.h>
#define int long long
#define mod 998244353
using namespace std;
int t,fac[200005],inv[200005],wop[200005],inv2,ans[200005],blk;
struct node{
	int n,m,id;
}q[200005];
int qpow(int a,int b){
	int r=1;
	while(b){
		if(b&1) r=r*a%mod;
		a=a*a%mod,b>>=1;
	}
	return r;
}
int caxy(int n,int m){
	if(m<0||m>n) return 0;
	return fac[n]*inv[m]%mod*inv[n-m]%mod;
}
bool cmp(node a,node b){
	int ba=a.n/blk,bb=b.n/blk;
	if(ba!=bb) return ba<bb;
	return (ba&1)?a.m>b.m:a.m<b.m;
}
signed main(){
	ios::sync_with_stdio(0);
	cin.tie(0);cout.tie(0);
	fac[0]=1;
	for(int i=1;i<200005;i++) fac[i]=fac[i-1]*i%mod;
	inv[200004]=qpow(fac[200004],mod-2);
	for(int i=200003;i>=0;i--) inv[i]=inv[i+1]*(i+1)%mod;
	wop[0]=1;
	for(int i=1;i<200005;i++) wop[i]=(wop[i-1]<<1)%mod;
	inv2=(mod+1)>>1;cin>>t;
	for(int i=1;i<=t;i++){
		int n,m;cin>>n>>m;
		q[i]={n-1,m-1,i};
	}
	blk=max((int)1,(int)(200000/sqrt(t)));
	sort(q+1,q+t+1,cmp);
	int n=0,m=0,now=1;
	for(int i=1;i<=t;i++){
		int qn=q[i].n,qm=q[i].m;
		while(n<qn) now=((now<<1)-caxy(n,m)+mod)%mod,n++;
		while(n>qn){
			n--;
			now=(now+caxy(n,m))%mod;
			now=now*inv2%mod;
		}
		while(m<qm) m++,now=(now+caxy(n,m))%mod;
		while(m>qm) now=(now-caxy(n,m)+mod)%mod,m--;
		ans[q[i].id]=(wop[q[i].n+1]-1+mod)%mod*now%mod;
	}
	for(int i=1;i<=t;i++) cout<<ans[i]<<"\n";
	return 0;
}
0