【XJOI娱乐赛】D.抗疫斗争

题目

http://115.236.49.52:83/contest/1351

思路

首先考虑如何计算hm
如果m为奇数，那么显然hm=1（你一个我一个，先手胜）
如果m为偶数，那么我们可以选一个数k，使得(m/k)%2==1，显然k最小为lowbit(m)
所以hm=lowbit(m)

下面我们来考虑计算f(n)的前缀和，则

们考虑枚举hi的取值，并统一计算其贡献。方便起见，我们记


复杂度O(sqrt(n))

代码

#include<bits/stdc++.h>
#define lowbit(x) ((x)&(-(x)))
#define ll long long
using namespace std;
const int mod = 998244353,N=1e6+77;
void inc(int &x,int y) 
{
	x=x+y>=mod?x+y-mod:x+y;
}
ll n;
int yjy,val[N];
int work(ll x)
{
	ll yjy=0;
	int sq=sqrt(x);
	for(int i=1; i<=sq; i++) yjy+=(x/i);
	yjy=(2*yjy-1ll*sq*sq)%mod;
	return yjy;
}
int main()
{
	scanf("%lld",&n);
	for(ll a=1; a<=n; a<<=1)
	{
		inc(yjy,1ll*(a-(a>>1))%mod*work(n/a)%mod);
	}
	cout<<yjy;
}