2024百度之星初赛第一场 跑步

跑步的终止是当所以有都同时到达终点，不难想到，这时一号跑的圈数是1~n的最小公倍数，设这个数为p，所以第i个人跑的圈数则是p/i。第i个人与第j个人碰面的次数为p/i-p/j。

这个题的难点有两个，第一个是在O（n）的时间求出1~n的最小公倍数p。

我们可以用质数筛，筛出1~n的所以质数，然后在求每一个质数的最大次数，这个时间复杂的接近O（n）。

第二个是求出每一个人跑的圈数。

这个我们需要用的乘法逆元的，一开始我直接用最小公倍数p乘i的逆元，每次都是用的这个公式：

// 求解 (a / b) % M  
ll ans = a * ksm(b, M - 2) % M;

结果超时了QWQ

后来有个大佬给我说可以O（n）求出1~n的所有数的逆元，代码如下：

void getly(int n)
{
	ly[0]=ly[1]=1;
	for(int i=2;i<=n;i++)
		ly[i]=1ll*(mod-mod/i)*ly[mod%i]%mod;
}

最后计算总的碰面次数就可以得到答案了

#include <bits/stdc++.h>
#define endl '\n'
#define int long long

using namespace std;

const int mod=998244353,N=1e7+5;

bool isprime[N];
int cnt,primes[N];
void getprimes(int n)//欧拉筛求质数
{
	for(int i=2;i<=n;i++)
	{
		if(!isprime[i]) primes[cnt++]=i;
		for(int j=0;primes[j]*i<=n;j++)
		{
			isprime[primes[j]*i]=1;
			if(i%primes[j]==0) break;
		}
	}
}

int n,a[N],s[N],ly[N];

void getly()//求1~n的逆元
{
	ly[0]=ly[1]=1;
	for(int i=2;i<=n;i++)
		ly[i]=1ll*(mod-mod/i)*ly[mod%i]%mod;
}

signed main()
{
	ios::sync_with_stdio(0),cin.tie(0),cout.tie(0);
	cin>>n;
	getprimes(n);
	int p=1;
	for(int i=0;i<cnt;i++)
	{
		int tmp=n;
		while(tmp/primes[i]) tmp/=primes[i],p=1ll*p*primes[i]%mod;
	}

	getly();
	
	for(int i=1;i<=n;i++)
		a[i]=1ll*p*ly[i]%mod;
	for(int i=n;i>=1;i--)
		s[i]=(a[i]+s[i+1])%mod;

	int ans=0;
	for(int i=1;i<=n;i++)
		(ans+=1ll*(n-i)*a[i]+mod-s[i+1])%=mod;//这里得加一个mod，因为减了后可能为负
	
	cout<<ans;
	
	return 0;
}