【数论】密码——扩展欧几里得

题目描述

有一个密码箱,0到n-1中的某些整数是它的密码。且满足:如果a和b都是它的密码,那么(a+b)%n也是它的密码(a,b可以相等,%表示整除取余数),某人试了k次密码,前k-1次都失败了,最后一次成功了。
问:该密码箱最多有多少不同的密码。

输入

第一行两个整数分别表示n,k(1≤k≤250000,k≤n≤1014)。第二行为k个用空格隔开的非负整数,表示每次试的密码。数据保证存在合法解。

输出

输出一行一个数,表示结果。

样例输入 Copy

42 5
28 31 10 38 24

样例输出 Copy

14

【分析】:参考文章

由题意得:若x是密码,则x的倍数一定是密码,kx%n一定是密码,即kx-tn=?这个正数一定是密码。

拓展欧几里得指出:kx-tn=gcd(x,n)一定存在可行的k和t使得式子成立,故gcd(x,n)一定是密码。

推论:若x不是密码,则gcd(x,n)一定不是(因为gcd(x,n)是x的因子),所以gcd(x,n)的因子也一定不是密码。

由推论可知 对于1<=i<k,gcd(a[i],n)的因子一定不是密码。

由于gcd(a[k],n)是密码,所以其因子只要不是gcd(a[i],n) (1<=i<k)的因子,就可以作为密码。

找到这样的最小的那个因子,将其作为密码,故其倍数全部是密码。
 

加上去重,避免超时

unique参考文章

unique的作用是“去掉”容器中相邻元素的重复元素,这里去掉要加一个引号,为什么呢,是因为它实质上是一个伪去除,它会把重复的元素添加到容器末尾,而返回值是去重之后的尾地址(是地址!!)

//#pragma GCC optimize(1)
//#pragma GCC optimize(2)
//#pragma GCC optimize(3,"Ofast","inline")
#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
ll a[250000];
ll n,k;
由于gcd(a[k],n)是密码,所以其因子只要不是gcd(a[i],n) (1<=i<k)的因子,就可以作为密码。
bool check(ll x)
{
    for(int i=1; i<=k; i++)
    {
        if(a[i]%x==0) return 0;
    }
    return 1;
}
ll gcd(ll a,ll b){ b==0?a:gcd(b,a%b); }
int main()
{
    scanf("%lld%lld",&n,&k);
    for(int i=1; i<=k; i++)
    {
        scanf("%lld",a+i);
        a[i]=__gcd(a[i],n); ///这里用内置的gcd就行啦
        // printf("%lld ",a[i]);
    }
    ll ak=a[k];
    sort(a+1,a+k);
    k=unique(a+1,a+k)-a-1;
    //一般使用前需要对容器进行排序,这样才能实现对整个数组去重。
    //cout<<endl;
    ll ans=0;
    for(int i=1; i*i<=ak; i++)
    {
        if(ak%i==0)
        {
            if(check(i))
            {
                ans=n/i;
                break;
            }
            else if(check(ak/i))
            {
                ans=n/ak*i;
            }
        }
    }
    printf("%lld",ans);
    return 0;
}

  • 0
    点赞
  • 2
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值