zoukankan      html  css  js  c++  java
  • poj 2115 C Looooops(推公式+扩展欧几里得模板)

    Description

    A Compiler Mystery: We are given a C-language style for loop of type 
    for (variable = A; variable != B; variable += C)
    
      statement;
    
    I.e., a loop which starts by setting variable to value A and while variable is not equal to B, repeats statement followed by increasing the variable by C. We want to know how many times does the statement get executed for particular values of A, B and C, assuming that all arithmetics is calculated in a k-bit unsigned integer type (with values 0 <= x < 2k) modulo 2k. 


    Input

    The input consists of several instances. Each instance is described by a single line with four integers A, B, C, k separated by a single space. The integer k (1 <= k <= 32) is the number of bits of the control variable of the loop and A, B, C (0 <= A, B, C < 2k) are the parameters of the loop. 
    
    The input is finished by a line containing four zeros. 

    Output

    The output consists of several lines corresponding to the instances on the input. The i-th line contains either the number of executions of the statement in the i-th instance (a single integer number) or the word FOREVER if the loop does not terminate. 

    Sample Input

    3 3 2 16
    3 7 2 16
    7 3 2 16
    3 4 2 16
    0 0 0 0

    Sample Output

    0
    2
    32766
    FOREVER

    Source

     

    解题思路:这道题和POJ1061(青蛙约会)一样,都是同余方程的求解,用到了拓展欧几里德算法。而本题题意明确,就是求解这个公式:(a+c*x)mod2^k=b ,求得x 的最小解。变形后可得:c*xmod2^k=b-a,即 c*x=(b-a)mod2^k; 这就是标准的同余方程。

    注意:k <=32 ,而 2的 32次方超出整数范围,所以要用__int64或long long ,就不会出现runtime error了。

     1 #pragma comment(linker, "/STACK:1024000000,1024000000")
     2 #include<iostream>
     3 #include<cstdio>
     4 #include<cstring>
     5 #include<cmath>
     6 #include<math.h>
     7 #include<algorithm>
     8 #include<queue>
     9 #include<set>
    10 #include<bitset>
    11 #include<map>
    12 #include<vector>
    13 #include<stdlib.h>
    14 using namespace std;
    15 #define max(a,b) (a) > (b) ? (a) : (b)  
    16 #define min(a,b) (a) < (b) ? (a) : (b)
    17 #define ll long long
    18 #define eps 1e-10
    19 #define MOD 1000000007
    20 #define N 1000000
    21 #define inf 1e12
    22 ll fac(ll m){
    23     ll ans=2;
    24     for(ll i=1;i<m;i++){
    25         ans=ans*2;
    26     }
    27     return ans;
    28 }
    29 ll e_gcd(ll a,ll b,ll &x,ll &y){
    30     if(b==0)
    31      {
    32          x=1;
    33          y=0;
    34          return a;
    35      }
    36      ll r=e_gcd(b,a%b,x,y);
    37      ll t=x;
    38      x=y;
    39      y=t-a/b*y;
    40      return r;
    41 }
    42 int main()
    43 {
    44     ll a,b,c,k;
    45     while(scanf("%I64d%I64d%I64d%I64d",&a,&b,&c,&k)==4){
    46         if(a==0 && b==0 && c==0 && k==0){
    47             break;
    48         }
    49         ll x,y,r;
    50         ll d=e_gcd(c,fac(k),x,y);
    51         //printf("---%I64d %I64d %I64d
    ",d,x,y);
    52         if((b-a)%d!=0){
    53             printf("FOREVER
    ");
    54         }
    55         else{
    56             x=x*(b-a)/d;
    57             r=fac(k)/d;
    58             x=(x%r+r)%r;
    59             printf("%I64d
    ",x);
    60         }
    61     }
    62     return 0;
    63 }
    View Code
  • 相关阅读:
    SQL基础-子查询&EXISTS&UNION
    SQL基础-汇总统计及GROUP BY
    SQL基础-创建新的输出字段
    SQL基础-过滤数据
    获取图片为二进制流,并且显示图片到网页
    计算两个日期相差的工作小时数,过滤了周末双休小时数
    计算两个日期相差的工作天数,过滤了周末双休天数
    WEB性能测试工具
    C#整合VS2010和NUnit
    如何打开Windows Server 2008 R2的域安全策略
  • 原文地址:https://www.cnblogs.com/UniqueColor/p/4808876.html
Copyright © 2011-2022 走看看