POJ 2115 C Looooops(扩展欧几里德)

原创 2016年05月31日 01:22:42

题目链接:
POJ 2115 C Looooops
题意:
有一个循环语句,变量从A开始,每次都要+C,终止条件是变量==B,每次累加都要mod 2^k.问循环最多执行几次?或者永远循环,
分析:
令K = 2^k.设循环执行了x次则:(A + C * x) mod K = B,也就是A + C * x = B + K * y.
移项可得:C * x - K * y = B - A.那么就是解这个不定方程的最小正整数解x。扩展欧几里德呀~

#include <iostream>
#include <cstdio>
#include <cstring>
#include <string>
#include <algorithm>
#include <climits>
#include <cmath>
#include <ctime>
#include <cassert>
#define IOS ios_base::sync_with_stdio(0); cin.tie(0);
using namespace std;
typedef long long ll;

int k;
ll A, B, C, K;

ll ex_gcd(ll a, ll b, ll& x, ll& y)
{
    if(b == 0) {
        x = 1, y = 0;
        return a;
    }
    ll d = ex_gcd(b, a % b, y, x);
    y -= a / b * x;
    return d;
}

int main()
{
    while(~scanf("%lld%lld%lld%d", &A, &B, &C, &k) && (A || B || C || k)){
        if(C == 0){
            if(A == B) printf("0\n");
            else printf("FOREVER\n");
            continue;
        }
        K = (ll)1 << k;
        ll x, y, d;
        d = ex_gcd(C, K, x, y);
        if((B - A) % d != 0) printf("FOREVER\n");
        else {
            x = x * (B - A) / d;
            ll r = K / d;
            x = (x % r + r) % r;
            printf("%lld\n", x);
        }
    }
    return 0;
}
版权声明:缥缈玉京人,想语然、京兆眉妩。

相关文章推荐

POJ 2115- C Looooops(扩展欧几里德)

C Looooops Time Limit:1000MS     Memory Limit:65536KB     64bit IO Format:%I64d & %I64u Subm...

POJ 2115 C Looooops(扩展欧几里德 + 求解模线性方程)

C Looooops Description A Compiler Mystery: We are given a C-language style for loop of typ...

POJ 2115-C Looooops-扩展欧几里德算法

C Looooops Time Limit: 1000MS   Memory Limit: 65536K Total Submissions: 21526   Accepted: 5868...

C Looooops(欧几里德+poj2115)

C Looooops Description A Compiler Mystery: We are given a C-language style for loop of type for (va...

poj2115 C Looooops (欧几里德)

C Looooops Time Limit: 1000MS   Memory Limit: 65536K Total Submissions: 20515   Ac...
  • d_x_d
  • d_x_d
  • 2015-09-23 21:04
  • 937

欧几里德 poj2115 C Looooops

关于拓展欧几里德,别人已经写的很好的了,我也自己写一下,方便以后自己复习。 其最基础的思想就是 gcd(a,b)=gcd(b,a mod b),其中 gcd的意思是求最大公约数。。 拓展欧几里德是...

poj2115 Looooops(扩展欧几里德模板)

poj2115 Looooops(扩展欧几里德模板) A Compiler Mystery: We are given a C-language style for loop of type for...

POJ 2115 Looooops 扩展欧几里德 Java

利用扩展欧几里德算法求解模线性方程!!!【点击蓝色字体,查看算法详情】import java.io.BufferedReader; import java.io.BufferedWriter; imp...

poj2115 Looooops 扩展欧几里德

x=[ (B-A+2^k)%2^k ]/C C*X=(B-A+2^k)%2^k 一开始 方程不知道咋办 后来 看到别人统一一下变量 顿时 愚蠢的我感觉清晰了一点点  令a=C, b=B-A, MOD=...

PKU 2115 C Looooops(拓展欧几里德)

原题链接:Here! 题意: Description: 对于循环语句:for(variable=A;variable!=B;variable+=C)statement; 已知所有的数均要m...
内容举报
返回顶部
收藏助手
不良信息举报
您举报文章:深度学习:神经网络中的前向传播和反向传播算法推导
举报原因:
原因补充:

(最多只允许输入30个字)