IT数码 购物 网址 头条 软件 日历 阅读 图书馆
TxT小说阅读器
↓语音阅读,小说下载,古典文学↓
图片批量下载器
↓批量下载图片,美女图库↓
图片自动播放器
↓图片自动播放器↓
一键清除垃圾
↓轻轻一点,清除系统垃圾↓
开发: C++知识库 Java知识库 JavaScript Python PHP知识库 人工智能 区块链 大数据 移动开发 嵌入式 开发工具 数据结构与算法 开发测试 游戏开发 网络协议 系统运维
教程: HTML教程 CSS教程 JavaScript教程 Go语言教程 JQuery教程 VUE教程 VUE3教程 Bootstrap教程 SQL数据库教程 C语言教程 C++教程 Java教程 Python教程 Python3教程 C#教程
数码: 电脑 笔记本 显卡 显示器 固态硬盘 硬盘 耳机 手机 iphone vivo oppo 小米 华为 单反 装机 图拉丁
 
   -> C++知识库 -> Lost Arithmetic Progression【Codeforces Round #785 (Div. 2)】 -> 正文阅读

[C++知识库]Lost Arithmetic Progression【Codeforces Round #785 (Div. 2)】

作者:recommend-item-box type_blog clearfix

题意
给出两个增等差数列B、C,请问有多少个等差数列A,使得A和B的交集是C。特殊地,当A有无限多个,输出-1。

思路
以下st,d,num,ed分别表示首项,公差,项数,末项。其中,ed = st + (num - 1)*ed.
1、首先应当判断什么时候不存在这样的A.

当存在C中的元素不在B中时,A一定不存在。即我们需要判断:是否所有C中元素在B中都有
首先,应当满足C首项大于等于B首项,C末项小于等于B末项。
其次,stc应当为stb + k*db。
最后,dc应当为db的倍数。

以上任一条件不满足,则存在C中元素不在B中。输出0.

2、判断什么时候输出-1.
在满足1的所有条件下,如果C数列的stc往前一项小于stb,则A数列可以无限往负无穷方向延展;同样地,如果C数列的edc往后一项大于edb,则A数列可以无限往正无穷方向延展。

这两种情况满足任一,输出-1,即有无限多种可能的A。

3、其他情况(明确计数).
由于C是A、B的交集,显然需要满足dc为da和db的最小公倍数,否则,A、B交集可能包含C中没有的元素,或者不完全包含C中的元素。
由于dc是da的(最小公)倍数,所以用 d c \sqrt{d_c} dc? ?的时间复杂度枚举dc的所有因数i(即枚举da),判断枚举的da和db的最小公倍数是否为dc即可。对于这样的da,它提供的答案数为 ( d c d a ) 2 (\frac{d_c}{d_a})^2 (da?dc??)2.

代码

#include<stdio.h>
#include<iostream>
#include<string.h>
#include<stdlib.h>
#include<algorithm>
#include<vector>
#include<queue>
#include<map>
#include<math.h>
using namespace std;
typedef long long ll;
ll stb,db,numb;
ll stc,dc,numc;

const int mod = 1e9+7;

ll gcd(ll a,ll b){
	return b == 0?a:gcd(b,a%b);
}

void solve(){
	scanf("%lld%lld%lld",&stb,&db,&numb);
	scanf("%lld%lld%lld",&stc,&dc,&numc);
	ll edb = stb + (numb - 1)*db;
	ll edc = stc + (numc - 1)*dc;
	if((stc < stb || edc > edb) || dc%db != 0 || ((stc - stb)%db != 0)){
		printf("0\n");
		return ;
	}
	
	if(stc - dc < stb || edc + dc > edb){
		printf("-1\n");
		return ;
	}
	
	ll res = 0;
	for(ll i = 1;i*i <= dc;i++){
		//i is da
		if(dc%i != 0)
			continue;
			
		if(i*db == gcd(i,db)*dc){
			res = (res + dc/i*dc/i)%mod;
		}
		
		if(i*i != dc){
			ll tda = dc/i;
			if(tda*db == gcd(tda,db)*dc){
				res = (res + dc/tda*dc/tda)%mod;
			}
		}
	}
	
	printf("%lld\n",res);
}

int main(){
	int t;
	scanf("%d",&t);
	while(t--){
		solve();
	}
	return 0;
}

  C++知识库 最新文章
【C++】友元、嵌套类、异常、RTTI、类型转换
通讯录的思路与实现(C语言)
C++PrimerPlus 第七章 函数-C++的编程模块(
Problem C: 算法9-9~9-12:平衡二叉树的基本
MSVC C++ UTF-8编程
C++进阶 多态原理
简单string类c++实现
我的年度总结
【C语言】以深厚地基筑伟岸高楼-基础篇(六
c语言常见错误合集
上一篇文章      下一篇文章      查看所有文章
加:2022-05-07 11:00:54  更:2022-05-07 11:01:50 
 
开发: C++知识库 Java知识库 JavaScript Python PHP知识库 人工智能 区块链 大数据 移动开发 嵌入式 开发工具 数据结构与算法 开发测试 游戏开发 网络协议 系统运维
教程: HTML教程 CSS教程 JavaScript教程 Go语言教程 JQuery教程 VUE教程 VUE3教程 Bootstrap教程 SQL数据库教程 C语言教程 C++教程 Java教程 Python教程 Python3教程 C#教程
数码: 电脑 笔记本 显卡 显示器 固态硬盘 硬盘 耳机 手机 iphone vivo oppo 小米 华为 单反 装机 图拉丁

360图书馆 购物 三丰科技 阅读网 日历 万年历 2024年11日历 -2024/11/30 12:07:44-

图片自动播放器
↓图片自动播放器↓
TxT小说阅读器
↓语音阅读,小说下载,古典文学↓
一键清除垃圾
↓轻轻一点,清除系统垃圾↓
图片批量下载器
↓批量下载图片,美女图库↓
  网站联系: qq:121756557 email:121756557@qq.com  IT数码