信友队算法随笔——数学1

最新推荐文章于 2024-07-01 09:54:05 发布

TangyixiaoQAQ

最新推荐文章于 2024-07-01 09:54:05 发布

阅读量982

点赞数 17

文章标签：算法 c++

本文链接：https://blog.csdn.net/DCMyyds/article/details/135004668

版权

A. 优美子数列2

题目描述：

数学家小 Q 得到了一个长度为 $n $的数列 ${a_n}$ 。

小 Q 的幸运数字是 $k$ ，所以他认为，若一个子数列 $a_l, a_{l+1}, …, a_r $的和为 $k $的倍数，则该子数列是优美子数列。

小 Q 现在想考考你， ${a_n}$ 里有多少个优美子数列呢？

输入格式：
第一行输入两个正整数 $n, k$ 。

接下来一行输入 $n $个数，代表 $a_1 $到 $a_n$ 。

输出格式：

输出一行一个正整数，表示优美子数列的个数。

样例输入：

5 3

1 2 3 4 5

样例输出：

数据规模：

对于 $ 60% $ 的数据， $n \leq 1000$ ；

对于 $100$ 的数据， $1≤n≤2×10^5$ ， $1≤k≤10^7$ ， $10^9≤a_i≤10^9$ 。

WA60：

#include<bits/stdc++.h>
using namespace std;
const int N=2e5+10;
int n,k,ar[N],br[100*N],ans;
map<int,int> mp;
int main(){
	cin>>n>>k;
	for(int i=1;i<=n;i++) {
		cin>>ar[i];
		ar[i]+=ar[i-1];
		int sum=(ar[i]%k+k)%k;
		br[sum]++;
		if(sum==0)ans+=br[sum];
		else ans+=(br[sum]-1);
	}
	cout<<ans;
	return 0;
}

AC Code：

#include<bits/stdc++.h>
using namespace std;
const int N=2e5+10;
long long n,k,ar[N],br[100*N],ans;
map<int,int> mp;
int main(){
	cin>>n>>k;
	for(int i=1;i<=n;i++) {
		cin>>ar[i];
		ar[i]+=ar[i-1];
		int sum=(ar[i]%k+k)%k;
		br[sum]++;
		if(sum==0)ans+=br[sum];
		else ans+=(br[sum]-1);
	}
	cout<<ans;
	return 0;
}

B.蒜头君的数轴

题目描述：

今天蒜头君拿到了一个数轴，上边有 n 个点，但是蒜头君嫌这根数轴不够优美，想要通过加一些点让它变优美，所谓优美是指考虑相邻两个点的距离，最多只有一对点的距离与其它的不同。

蒜头君想知道，他最少需要加多少个点使这个数轴变优美。

输入格式：

输入第一行为一个整数 $n (1≤n≤10^5)$ ，表示数轴上的点数。

第二行为 n 个不重复的整数 $x_1, x_2, ..., x_n (−10^9≤x_i≤10^9)$ ，表示这些点的坐标，点坐标乱序排列。

输出格式：

输出一行，为一个整数，表示蒜头君最少需要加多少个点使这个数轴变优美。

样例输入：

1 3 7 15

样例输出：

数据规模：

对于 100% 的数据， $1≤n≤10^5$ ， $10^9≤x_i≤10^9$ 。

AC Code:

#include <bits/stdc++.h>
#define int long long
using namespace std;
const int N= 1e5+10;
int n,ar[N],br[N],l[N],r[N],res ,ans;
inline void solve() {
	cin>>n;
	for(int i=1; i<=n; i++) cin>>ar[i];
	sort(ar+1,ar+1+n);
	if(n<=2) {
		cout<<"0\n";
		return;
	}
	for(int i=2; i<=n; i++) {
		br[i-1] = ar[i] - ar[i-1];
		ans += br[i-1];
	}
	l[1] = br[1];
	r[n-1] =br[n-1];
	for(int i=2; i<=n-1; i++) l[i] = __gcd(l[i-1],br[i]);
	for(int i=n-2; i>=1; i--) r[i] =__gcd(r[i+1],br[i]);
	res = min((ans-br[1])/r[2],(ans-br[n-1])/l[n-2]);
	for(int i=2; i<=n-2; i++) {
		res= min(res,(ans-br[i])/__gcd(l[i-1],r[i+1]));
	}
	cout<<res-(n-2)<<"\n";
}
signed main() {
	ios::sync_with_stdio(false);
	cin.tie(0);
	solve();
	return 0;
}

思路：

先看，“点坐标乱序排列“，因为是数轴，所以就要保证从左到右有序。
也就是是用sort从小到大排序

不难看出，其实当 $n \leq 2$ 时，我们并不需要添加顶点，所以直接输出0即可

相邻的点的距离一定是原来点距离的 $gcd() $ 中的最大值，可以预处理出 $gcd() $的前缀和后缀，直接用前缀 $gcd() $后缀找到最大值

C. 最长整除序列

题目描述

给一个正整数 $x$ ,让你构造出一个长度为 $m + 1$ 的序列： $a_0,a_1,a_2,...a_m$
其中 $a_0=1,a_m=x,a_i<a_{i+1}$ ，并且 $a_i|a_{i+1}$

$a ∣ b$ 的意思是 $b$ 能被 $a$ 整出，比如 $3∣6$

求能构造出的序列中 $m$ 的最大值,以及 $m$ 等于最大值时的方案数。

输入格式

第一行包含一个正整数 $T(1<T<10^5)$ 。

接下来 $T$ 行，每行包括一个正整数 $x(2<x<10^6)$

输出格式

对于每组数据,输出一行,包含2个整数: $m$ 的最大值,以及 $m$ 等于最大值时的方案数。两个整数按空格隔开。

样例

Input 1

100

Output 1

1 1

2 1

2 2

4 6

思路

求余后组合

AC Code

#include <bits/stdc++.h>
using namespace std;
int f[25][25];
inline void init() {
	f[0][0]=1;
	for(int i=1; i<24; i++) {
		f[i][0] = 1;
		for(int j=1; j<=i; j++) {
			f[i][j] = f[i-1][j] + f[i-1][j-1];
		}
	}
}
inline void solve() {
	int x;
	cin>>x;
	int res = 0,ans = 1;
	for(int i=2; i*i<=x; i++) {
		if(x%i==0) {
			int tmp=0;
			while(x%i==0) {
				x/=i;
				tmp++;
			}
			res += tmp;
			ans *= f[res][tmp];
		}
	}
	if(x>1) {
		res +=1;
		ans *= f[res][1];
	}
	cout<<res<<" "<<ans<<"\n";
}
signed main() {
	ios::sync_with_stdio(false);
	cin.tie(0);
	init();
	int _;
	cin>>_;
	while(_--) {
		solve();
	}
	return 0;
}

D. 世纪之花

题目描述：

“丛林变得焦躁不安…”

世纪之花触手正在向丛林的各个角落蔓延。

现在丛林中已经有 $x$ 个触手，每过一分钟世纪之花会长出一些新触手，新触手的数量等于当前触手数的最小质因子。

勇者准备出发，去击败世纪之花。勇者想要知道，再经过几分钟，世纪之花的触手数量就不小于 $y$ 了呢？

输入格式：

第一行包含一个整数 $T$ ，表示数据组数。

接下来 $T$ 行，每行两个正整数 $x$ 和 $y$ 。

输出格式：

输出 $T$ 行，每行一个整数表示答案。

样例输入：

2 23

9 20

5 100

6 89

样例输出：

数据规模：

$1 \leq T \leq 1000$ ， $2 \leq X \leq 10$ ， $20≤Y≤10^9$

2023.05.31 UPD：数据已加强

思路：

虽然有的人说过，特判能过。

但我还是倾向于「正义」的做法

其实很简单（我喜欢简洁！

因为 $x \leq 10$ ，
所以直接判断2，3，5，7即可

AC Code

#include <bits/stdc++.h>
#define int long long
using namespace std;
int ar[4] = {2,3,5,7};
inline void solve() {
	int x,y;
	cin>>x>>y;
	if(x>=y) {
		cout<<"0\n";
		return;
	}
	int ans =1;
	for(int i=0; i<4; i++) {
		if(x%ar[i]==0) {
			x+=ar[i];
			break;
		}
	}
	if(y>x) {
		ans+=(y-x+1)/2;
	}
	cout<<ans<<"\n";
}
signed main() {
	ios::sync_with_stdio(false);
	cin.tie(0);
	int _;
	cin>>_;
	while(_--) {
		solve();
	}
	return 0;
}

E. 参天大树

题目描述：

丛林中矗立着一棵参天大树，高度为 $y$ 。

大树 $2$ ~ $p$ 高度处，每个位置都有一只蚱蜢。

一只蚱蜢如果在 x 高度处，那么它可以跳到 $2 x$ 、 $3 x$ 、 $4 x$ 、 $......$ 等任意一个 $x$ 的倍数处。

你想在 $2$ ~ $y$ 高度范围内，找到一个尽可能高且不会有任何蚂蚱能跳到的位置。

注意你和蚂蚱的位置都必须是整数，如果你找不到任何一个合法的位置的话，就输出 $- 1$ 。

输入格式：

第一行包含两个正整数 $p$ 和 $y$ 。

输出格式：

输出一行一个整数表示答案。

样例输入：

3 6

样例输出：

样例输入：

3 4

样例输出：

-1

数据规模：

$2$ $\leq$ $p$ $\leq$ $y$ $\leq$ $10^9$

思路

直接筛即可

AC Code

#include <bits/stdc++.h>
#define int long long
using namespace std;
int p,y;
inline bool check(int x) {
	for(int i=2; i<=sqrt(x); i++) {
		if(x%i==0&&i>=2&&i<=p) return false;
	}
	return true;
}
inline void solve() {
	cin>>p>>y;
	int z = sqrt(y);
	for(int i=y; i>=p+1; i--) {
		if(check(i)) {
			cout<<i<<"\n";
			return;
		}
	}
	cout<<"-1\n";
}
signed main() {
	ios::sync_with_stdio(false);
	cin.tie(0);
	solve();
	return 0;
}

F. 小夏的大师旅程

题目描述

小夏已经踏上了成为最伟大的口袋妖怪大师的旅程。为了得到他的第一个口袋妖怪，他去了Zulu教授的实验室。由于小夏是Zulu教授最喜欢的学生，Zulu允许他从实验室里取出任意数量的口袋妖怪。

但是Zulu警告他，每个小精灵都有一个力量值，例如k(k>1)个小精灵在一起，它们的力量值为 ${s1, s2, s3,…, sk}$ ，如果 $g c d (s 1, s 2, s 3, \dots, s k) = 1$ ，它们之间就会互相打架。

小夏作为一个聪明的人，不希望他的口袋妖怪互相斗争。然而，他也想最大化他从实验室里带走的神奇宝贝的数量。你能帮小夏找出他能带走的最大数量的口袋妖怪吗？

注意：口袋妖怪不能与自己战斗。

输入格式

输入包含两行。第一行一个整数 $n(1<=n<=10^5)$ ，表示实验室中的小精灵总数。第二行n个用空格隔开的整数，第i个整数代表第i个小精灵的力量值 $s_i(1<=s_i<=10^5)$ 。

输出格式

一行包含一个整数，表示能拿走的小精灵数量最大值。

样例

Input 1

2 3 4

Output 1

Input 2

2 3 4 6 7

Output 2

数据范围

输入的小精灵总数n的范围是 $1$ 到 $10^5$ ，小精灵的力量值 $s_i$ 的范围是 $1$ 到 $10^5$ 。

样例解释

对于样例1，小夏可以选择第2号和第4号小精灵，它们的力量值为 ${2, 4}$ ， $g c d (2, 4) = 2$ ，它们不会互相打架。所以最多可以拿走2个小精灵。

对于样例2，小夏可以选择第1号、第3号和第4号小精灵，它们的力量值为 ${2, 4, 6}$ ， $g c d (2, 4, 6) = 2$ ，它们不会互相打架。所以最多可以拿走3个小精灵。

思路

rt即可

AC Code

#include <bits/stdc++.h>
using namespace std;
const int N = 1e5+10;
int ar[N],n;
inline void check(int x) {
	for(int i=2; i*i<=x; i++) {
		if(x%i==0) {
			ar[i]++;
			while(x%i==0) x/=i;
		}
	}
	if(x>1) ar[x]++;
}
inline void solve() {
	cin>>n;
	for(int i=1,a; i<=n; i++) {
		cin>>a;
		check(a);
	}
	int res =1;
	for(int i=2; i<=N; i++) {
		res = max(res,ar[i]);
	}
	cout<<res<<"\n";
}
signed main() {
	ios::sync_with_stdio(false);
	cin.tie(0);
	solve();
	return 0;
}

G. 连续互质

题目描述：

小信有两个正整数 $(a, b)$ ，其中 $a < b$ 。

现在小信想找一个最大的正整数 $x$ , $x > 0$ ，使得对于所有 $k = 0$ ， $x$ $-$ $1$ 满足 $gcd(a+k,b+k) $都为 $1$ 。

如果没有满足条件的 $x$ ，输出 $0$ ；如果 $x$ 无限大，输出 $-$ $1$ 。

输入格式：

第一行包含一个整数
$T$ ，表示测试数据组数。

对于每组测试数据，包含两个整数 $
a,
b$。

输出格式：

对于每组测试数据，输出一个整数表示答案。

样例1输入：

3 9

1 3

8 9

10009 20000

样例1输出：

-1

约定与提示：

对于100%的数据，
$1≤T≤10^6$ ， $1≤a<b≤10^7$ 。

对于样例1的第一组测试数据：

$g c d (3, 9) = 3$ ，已经不满足条件，所以输出
0
。

对于样例1的第二组测试数据：

$g c d (1, 3) = 1$ ， $g c d (2, 4) = 2$ ，所以输出1。

思路：

先筛后捡

AC Code

#include <bits/stdc++.h>
//#define int long long
using namespace std;
const int N=1e7+10,mod = 1e9+7;
bool st[N];
int a,b,idx;
int ar[N];
vector<int> prime;
inline void init() {
	for(int i=2; i<N; i++) {
		if(!st[i]) prime.push_back(i),ar[i] = i,idx++;
		for(int j=0; j<idx&&prime[j]*i<N; j++) {
			st[prime[j]*i]= true;
			ar [prime[j]*i] = prime[j];
			if(i%prime[j]==0)break;
		}
	}
}
inline void solve() {
	cin>>a>>b;
	if(__gcd(a,b)>1) cout<<"0\n";
	else if(b-a==1) cout<<"-1\n";
	else {
		int num =b-a;
		int res= 1e9;
		while(num>1) {
			int k=(ar[num] - a%ar[num]);
			res = min(res,k);
			num/=ar[num];
		}
		cout<<res<<"\n";
	}
}
signed main() {
	ios::sync_with_stdio(false);
	cin.tie(0);
	init();
	int _;
	cin>>_;
	while(_--) {
		solve();
	}
	return 0;
}