codeforces 894C Marco and GCD Sequence

本文介绍了一种解决特定数学问题的方法,即如何构造一个数列使得该数列的所有子区间的最大公约数(GCD)恰好构成给定集合。文章详细阐述了通过找到集合内元素的最大公约数,并检查其是否存在于原始集合中的方式来解决此问题。

http://www.elijahqi.win/archives/1629

In a dream Marco met an elderly man with a pair of black glasses. The man told him the key to immortality and then disappeared with the wind of time.

When he woke up, he only remembered that the key was a sequence of positive integers of some length n, but forgot the exact sequence. Let the elements of the sequence be a1, a2, …, an. He remembered that he calculated gcd(ai, ai + 1, …, aj) for every 1 ≤ i ≤ j ≤ n and put it into a set S. gcd here means the greatest common divisor.

Note that even if a number is put into the set S twice or more, it only appears once in the set.

Now Marco gives you the set S and asks you to help him figure out the initial sequence. If there are many solutions, print any of them. It is also possible that there are no sequences that produce the set S, in this case print -1.

Input
The first line contains a single integer m (1 ≤ m ≤ 1000) — the size of the set S.

The second line contains m integers s1, s2, …, sm (1 ≤ si ≤ 106) — the elements of the set S. It’s guaranteed that the elements of the set are given in strictly increasing order, that means s1 < s2 < … < sm.

Output
If there is no solution, print a single line containing -1.

Otherwise, in the first line print a single integer n denoting the length of the sequence, n should not exceed 4000.

In the second line print n integers a1, a2, …, an (1 ≤ ai ≤ 106) — the sequence.

We can show that if a solution exists, then there is a solution with n not exceeding 4000 and ai not exceeding 106.

If there are multiple solutions, print any of them.

Examples
Input
4
2 4 6 12
Output
3
4 6 12
Input
2
2 3
Output
-1
Note
In the first example 2 = gcd(4, 6), the other elements from the set appear in the sequence, and we can show that there are no values different from 2, 4, 6 and 12 among gcd(ai, ai + 1, …, aj) for every 1 ≤ i ≤ j ≤ n.

题意:要求我们构造一个数列 满足这里面所有区间的最小公因数都在这里出现过 并且不会有没有出现过的 首先我把所有gcd做一个gcd求出他们的最大公因数 然后 我判断这个最大公因数 是否在给定的数列中出现 如果没有 则一定是输出-1 否则则我在原数列中没两个中间加一个他们的最大公因数即可

#include<cstdio>
#define N 1100
inline char gc(){
    static char now[1<<16],*S,*T;
    if (T==S){T=(S=now)+fread(now,1,1<<16,stdin);if (T==S) return EOF;}
    return *S++;
}
inline int read(){
    int x=0;char ch=gc();
    while (ch<'0'||ch>'9') ch=gc();
    while (ch<='9'&&ch>='0'){x=x*10+ch-'0';ch=gc();}
    return x;
}
inline int gcd(int x,int y){
    if (y==0) return x;return gcd(y,x%y);
}
int a[N],n;
int main(){
    freopen("cf.in","r",stdin);
    n=read();int gg=0;bool flag=0;
    for (int i=1;i<=n;++i) a[i]=read(),gg=gcd(gg,a[i]);
    for (int i=1;i<=n;++i) if (a[i]==gg) {flag=1;break;}
    if (!flag) {printf("-1");return 0;}printf("%d\n",(n<<1)-1);
    for (int i=1;i<n;++i) printf("%d %d ",a[i],gg);printf("%d",a[n]);
    return 0;
}
### 关于Codeforces中的GCD问题 在Codeforces平台上存在多个涉及最大公约数(GCD)概念的问题。其中一道具有代表性的题目是编号为1025B的“Weakened Common Divisor”,该题由著名数学家Ildar引入了一个新的概念——弱化公因数(WCD),即对于一系列整数对列表而言的一种特殊性质[^2]。 具体到这道题目的描述如下:给出一个长度为\(n\)的数组\(a\),目标是在所有元素上加上同一个常量\(d\)之后能够找到至少两个不同的位置其值的最大公约数大于等于2,并且要使这个加上的常量尽可能小。此题的关键在于通过计算相邻两数之差来间接获取可能存在的公共因子,进而利用这些信息推导出满足条件所需的最小增量\[d\][^4]。 为了高效解决这类基于GCD的问题,在算法设计方面通常会采用一些特定技巧: - **差分遍历**:通过对原始序列做适当变换简化问题结构; - **快速求解GCD**:借助欧几里得算法迅速定位潜在候选者; - **优化查找过程**:针对所得结果进一步筛选最优方案; 下面是一个Python版本的解决方案片段用于演示如何处理上述提到的任务逻辑: ```python from math import gcd from itertools import pairwise def min_operations_to_weak_gcd(nums): diff_gcd = 0 for prev, curr in pairwise(nums): diff_gcd = gcd(diff_gcd, abs(curr - prev)) if diff_gcd == 1: return -1 factors = get_factors(diff_gcd) result = float('inf') target_modulo = nums[0] % diff_gcd for factor in factors: candidate = ((target_modulo + diff_gcd - (nums[0] % factor)) % factor) result = min(result, candidate) return int(result) def get_factors(n): """Helper function to generate all divisors.""" res = [] i = 1 while i*i <= n: if n % i == 0: res.append(i) if i != n // i: res.append(n//i) i += 1 return sorted(res)[::-1] ```
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值