2024百度之星 跑步

原题链接:码题集OJ-跑步

题目大意:一个n个人在绕圈跑,第i个人跑一圈的时间是i分钟,每二个人位置相同就会打一次招呼,如果同时来到终点,他们就会停下来,请问会打多少次招呼?

思路:首先可以想到这n个人会跑他们的最小公倍数的圈数之后停下来。最小公倍数用ores代替,如何求最小公倍数呢?一个数肯定是由一堆质数相乘得到的,所以只要求出1-n中每个质数的最高次幂就可以了,例如说要求1 2 3 4 5 6 7 8 9 10的最小公倍数,那么其实就是求1 1 1 1 5 1 7 8 9 1的最小公倍数,因为8=2*2*2,那么2的这个质数本身就不重要了。p字母为质数,那么这个质数的最高次幂就是\log n /\log p

因为跑的快的不会被跑的慢的人追到,那么可以想到一个必定超时的方案,那就是用二层for来枚举。对于第i个人,他前面的所有人都会被他追到,所有第i个人的贡献就是ores*n/i-ores*n/(i+1)+ores*n/i-ores*n/(i+2)......双重循环枚举就可以了,但是数据范围明显会超时,可以想到在双重枚举的过程中肯定会有很多不必要的计算,一个人可以追前面的人,也可以被后面的人追上,如果是追前面的人,那么n/i是减数,一共有(n-i)个人可以被追上,如果是被追上那么n/i就是被减数,一共有(i-1)个人,那么减数减去被减数的数量乘上当前数跑的圈数就是打招呼的数量也就是(n-2*i+1)*(n/i)。例如说1 2 3,如果双重循环计算,第一个人的贡献是:6-1/2*6+6-1/3*6,第二个人的贡献是:1/2*6-1/3*6,如果单独计算,那么第一个的贡献就是:6*(3-2*1+1),第二个人的贡献就是:6/2*(3-2*2+1)

那这样题目就很明显了,但是因为数据会很大要取模,所以需要算出从1-n的所有数的逆元。对于1-n的逆元可以线性的求出。

inv数组表示逆元

p=q*i+r 

二边同时mod p

0=q*i+r

二边同时乘上i的逆元和r的逆元

0=q*r^{-1}+i^{-1}

移项变形

i^{-1}=-q*r^{-1}=-p/i*pmodi^{-1} 

q=p/i,r^-1=(p%i)^-1,因为是mod意义下的计算,所以右边可以加上p*(p%i)^-1.

最终就是inv[i]=(mod-mod/i)*inv[mod%i]%mod.

//冷静,冷静,冷静
//调不出来就重构
#pragma GCC optimize(2)
#pragma GCC optimize("O3")
#include<bits/stdc++.h>
#define endl '\n'
using namespace std;
typedef long long ll;
typedef long double ld;
typedef pair<ll,ll> pii;
const int N=1e7+10,mod=998244353;
ll inv[N],prime[N],n;
bool vis[N]; 
ll ksm(ll a,ll b)
{
	ll ans=1;
	do
	{
		if(b&1)ans*=a;
		a*=a;b>>=1;
		a%=mod;ans%=mod;
	}while(b); 
	return ans;
}
int main()
{
	ios::sync_with_stdio(0),cin.tie(0),cout.tie(0);
	inv[1]=1;
	cin>>n;
	ll cnt=0,ores=1;
	for(int i=2;i<=n;i++)
	{
		inv[i]=(mod-mod/i)*inv[mod%i]%mod;//求每个数的逆元 
		if(!vis[i])prime[cnt++]=i,ores=ores*ksm(i,log(n)/log(i))%mod;//求n范围内的质数的最高次幂的乘积 
		for(int j=0;j<cnt&&i*prime[j]<=n;j++)
		{
			vis[i*prime[j]]=1;
			if(i%prime[j]==0)break;
		}
	}
	ll ans=0;
	for(int i=1;i<=n;i++)
	{
		ll op=ores*inv[i]%mod;//这个人跑的圈数 
		ans=(ans+op*(n-2*i+1)%mod+mod)%mod; 
	}
	cout<<ans;
    return 0;
}

相关推荐

  1. 2024题目记录

    2024-06-07 01:24:04       85 阅读
  2. 2022题目记录

    2024-06-07 01:24:04       35 阅读
  3. 2024第三场第一题 数星星

    2024-06-07 01:24:04       27 阅读
  4. Create2024AI开发者大会记录

    2024-06-07 01:24:04       34 阅读

最近更新

  1. docker php8.1+nginx base 镜像 dockerfile 配置

    2024-06-07 01:24:04       94 阅读
  2. Could not load dynamic library ‘cudart64_100.dll‘

    2024-06-07 01:24:04       100 阅读
  3. 在Django里面运行非项目文件

    2024-06-07 01:24:04       82 阅读
  4. Python语言-面向对象

    2024-06-07 01:24:04       91 阅读

热门阅读

  1. 黑马es数据同步mq解决方案

    2024-06-07 01:24:04       20 阅读
  2. 如何快速入门使用Vue.js

    2024-06-07 01:24:04       32 阅读
  3. 在Linux/Ubuntu/Debian系统中使用 `tar` 压缩文件

    2024-06-07 01:24:04       26 阅读
  4. Attention as an RNN

    2024-06-07 01:24:04       28 阅读
  5. 【C++】C++程序的四个区和智能指针的实现

    2024-06-07 01:24:04       33 阅读
  6. AWS需要实名吗?

    2024-06-07 01:24:04       28 阅读
  7. Android基础-Log管理

    2024-06-07 01:24:04       26 阅读
  8. Git命令行

    2024-06-07 01:24:04       26 阅读
  9. linux:scp用法详解

    2024-06-07 01:24:04       26 阅读
  10. 从零学算法30

    2024-06-07 01:24:04       30 阅读
  11. 区块链学习记录01

    2024-06-07 01:24:04       31 阅读