UVA12169

思路:枚举a,根据扩展欧几里得算法算出b,之后检查此组a和b是否正确。

#include<iostream>
using namespace std;
typedef long long LL;

const int maxn = 100 + 10;
const int maxv = 10001;
int x[2*maxn], T;

void gcd(LL a, LL b, LL& d, LL& x, LL& y) {  //扩展欧几里得辗转相除法
    if(!b) { d = a; x = 1; y = 0; }
    else { gcd(b, a % b, d, y, x); y -= x*(a/b); }
}

int main() {
    cin >> T;
    for(int i = 1; i < 2*T; i += 2) cin >> x[i];
    for(LL a = 0; a < maxv; a++) {
        LL b, k, d, t = x[3] - a*a*x[1];  //决定a, t必须为LL
        gcd(a+1, maxv, d, b, k);
        if(t % d) continue;
        b = b*t/d;  //决定b必须为LL
        b = (b % maxv + maxv) % maxv;  //保证b为正且小于maxv的值
        //cout << b << endl;
        //cout << a << endl;
        bool ok = true;
        for(int i = 2; i <= 2*T; i++)
            if(i % 2) {
                if(x[i] != (a * x[i - 1] + b) % maxv) {  //Check
                    ok = false;
                    break;
                }
            } else x[i] = (a * x[i - 1] + b) % maxv;  //赋值
        if(!ok) continue;
        for(int i = 2; i <= 2*T; i += 2) cout << x[i] << endl;
        break;
    }
    return 0;
}
  • 1
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值