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

A. 优美子数列2

题目描述:

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

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

小 Q 现在想考考你, a n {a_n} an 里有多少个优美子数列呢?

输入格式:
第一行输入两个正整数 n , k n,k n,k

接下来一行输入 $n $个数,代表 $a_1 $到 a n a_n an

输出格式:

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

样例输入:

5 3

1 2 3 4 5

样例输出:

7

数据规模:

对于 $ 60% $ 的数据, n ≤ 1000 n≤1000 n1000

对于 100 100% 100 的数据, 1 ≤ n ≤ 2 × 1 0 5 1≤n≤2×10^5 1n2×105 1 ≤ k ≤ 1 0 7 1≤k≤10^7 1k107 − 1 0 9 ≤ a i ≤ 1 0 9 -10^9≤a_i≤10^9 109ai109

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 ≤ 1 0 5 ) n (1≤n≤10^5) n(1n105),表示数轴上的点数。

第二行为 n 个不重复的整数 x 1 , x 2 ​ , . . . , x n ​ ( − 1 0 9 ≤ x i ​ ≤ 1 0 9 ) x_1, x_2​, ..., x_n​ (−10^9≤x_i​≤10^9) x1,x2,...,xn(109xi109),表示这些点的坐标,点坐标乱序排列

输出格式:

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

样例输入:

4

1 3 7 15

样例输出:

1

数据规模:

对于 100% 的数据, 1 ≤ n ≤ 1 0 5 1≤n≤10^5 1n105 − 1 0 9 ≤ x i ≤ 1 0 9 -10^9≤x_i≤10^9 109xi109

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 ≤ 2 n≤2 n2时,我们并不需要添加顶点,所以直接输出0即可

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

C. 最长整除序列

题目描述

给一个正整数 x x x,让你构造出一个长度为 m + 1 m+1 m+1的序列: a 0 , a 1 , a 2 , . . . a m a_0,a_1,a_2,...a_m a0,a1,a2,...am
其中 a 0 = 1 , a m = x , a i < a i + 1 a_0=1,a_m=x,a_i<a_{i+1} a0=1,am=x,ai<ai+1,并且 a i ∣ a i + 1 a_i|a_{i+1} aiai+1

a ∣ b a|b ab的意思是 b b b能被 a a a整出,比如 3 ∣ 6 3|6 3∣6

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

输入格式

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

接下来 T T T行,每行包括一个正整数 x ( 2 < x < 1 0 6 ) x(2<x<10^6) x(2<x<106)

输出格式

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

样例

Input 1

5

2

3

4

10

100

Output 1

1 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 x x 个触手,每过一分钟世纪之花会长出一些新触手,新触手的数量等于当前触手数的最小质因子。

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

输入格式:

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

接下来 T T T 行,每行两个正整数 x x x y y y

输出格式:

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

样例输入:

4

2 23

9 20

5 100

6 89

样例输出:

11

5

46

42

数据规模:

1 ≤ T ≤ 1000 1≤T≤1000 1T1000 2 ≤ X ≤ 10 2≤X≤10 2X10 20 ≤ Y ≤ 1 0 9 20≤Y≤10^9 20Y109

2023.05.31 UPD:数据已加强

思路:

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

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

其实很简单(我喜欢简洁!

因为 x ≤ 10 x≤10 x10
所以直接判断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 y y

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

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

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

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

输入格式:

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

输出格式:

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

样例输入:

3 6

样例输出:

5

样例输入:

3 4

样例输出:

-1

数据规模:

2 2 2 ≤ ≤ p p p ≤ ≤ y y y ≤ ≤ 1 0 9 10^9 109

思路

直接筛即可

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)个小精灵在一起,它们的力量值为 s 1 , s 2 , s 3 , … , s k {s1, s2, s3,…, sk} s1,s2,s3,,sk,如果 g c d ( s 1 , s 2 , s 3 , … , s k ) = 1 gcd(s1, s2, s3,…, sk)=1 gcd(s1,s2,s3,,sk)=1,它们之间就会互相打架。

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

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

输入格式

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

输出格式

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

样例

Input 1

3

2 3 4

Output 1

2

Input 2

5

2 3 4 6 7

Output 2

3

数据范围

输入的小精灵总数n的范围是 1 1 1 1 0 5 10^5 105,小精灵的力量值 s i s_i si的范围是 1 1 1 1 0 5 10^5 105

样例解释

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

对于样例2,小夏可以选择第1号、第3号和第4号小精灵,它们的力量值为 2 , 4 , 6 {2, 4, 6} 2,4,6 g c d ( 2 , 4 , 6 ) = 2 gcd(2, 4, 6)=2 gcd(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) (a,b),其中 a < b a<b a<b

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

如果没有满足条件的 x x x,输出 0 0 0;如果 x x x 无限大,输出 − − 1 1 1

输入格式:

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

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

输出格式:

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

样例1输入:

4

3 9

1 3

8 9

10009 20000

样例1输出:

0

1

-1

79

约定与提示:

对于100%的数据,
1 ≤ T ≤ 1 0 6 1≤T≤10^6 1T106 1 ≤ a < b ≤ 1 0 7 1≤a<b≤10^7 1a<b107

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

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

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

g c d ( 1 , 3 ) = 1 gcd(1,3)=1 gcd(1,3)=1 g c d ( 2 , 4 ) = 2 gcd(2,4)=2 gcd(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;
}

总结

在数论中,开 l o n g long long l o n g long long十分重要。但也要认真审题,看清数据范围。

数论最主要的是数学的公式,所以要仔细分析,眼见未必为实

  • 17
    点赞
  • 18
    收藏
    觉得还不错? 一键收藏
  • 1
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值