雪黛依梦 2010-09-11 22:00 发表评论

雪黛依梦 2010-09-06 21:06 发表评论

雪黛依梦 2010-09-06 20:31 发表评论

雪黛依梦 2010-08-29 11:10 发表评论

//满足关系：(x + m * s) - (y + n *s) = k*l  (k= 0 1 2 )即：(n - m)*s + l*k = x-y; 
//利用拓展的欧几里得解出可能的s 
#include <stdio.h>
#include <stdlib.h>

//求最大公约数 
__int64 gcd (__int64 a, __int64 b) 
{
    if (b == 0)
       return a;
       else
       gcd (b, a % b);
}

//求得满足 a*x + b*y = d；的x  y 
__int64 ex_Gcd (__int64 a, __int64 b, __int64 &x1, __int64 &y1)
{
    if ( b == 0 )
    {
         x1 = 1;
         y1= 0;
         return a;
    }
    
    int r = ex_Gcd( b, a%b, x1, y1);
    
    int temp = x1;
    x1 = y1;
    y1 = temp - a/b * y1;
}

int main ()
{
    __int64  x, y, m, n, l;
    int a, b, product, d;//a b 的最大公约数 product 
    int s, k;
    __int64 s1, k1, s2, k2;
    
    while ( scanf ("%I64d%I64d%I64d%I64d%I64d", &x, &y, &m, &n, &l ) != EOF )
    {
          a = n - m;
          b = l;
          d = x-y;
          product = gcd (a,b);
          
          if ( d % product != 0 )
             printf ("Impossible\n");
             else
             {
                 a = a / product;
                 b = b / product;
                 d = d / product;
                 ex_Gcd (a, b, s1, k1);               //得到(n-m)/product * s + l/product * k = 1;的 s k的解
                 s2 = d * s1;                     //得到(n-m)/product * s + l/product * k = d;的 s k的解
                 k2 = d * k1; 
                 
                 int t;
                 //s = s2 - b * t;    用下面的方法处理满足条件的解 
                 //k = k2 - a * t;
                 
                 t = s2 / b;
                 s = s2 - b * t;
                 if ( s <= 0)
                 {
                      s += b;
                 } 
                 printf ("%d\n", s);
             }
    }
    //system ("pause");
    return 0 ;
}

雪黛依梦 2010-08-28 22:46 发表评论

雪黛依梦 2010-08-28 22:05 发表评论

雪黛依梦 2010-08-28 21:58 发表评论

雪黛依梦 2010-08-28 21:32 发表评论

雪黛依梦 2010-08-06 14:50 发表评论

雪黛依梦 2010-08-06 09:19 发表评论

 1//为什么要把范围缩小到n的方根之内呢？？ 
 2#include <stdio.h>
 3#include <stdlib.h>
 4#include <math.h>
 5# define MAXSIZE 40000
 6
 7int prime[MAXSIZE];
 8int primeElem[MAXSIZE];
 9int cn;          //0 -- n内素数的个数 
10
11void is_prime ( )
12{
13     for (int i = 2; i < MAXSIZE; i ++)
14     {
15         if ( !prime [i] )
16         {
17              for (int j = 2 * i; j < MAXSIZE; j += i)
18              {
19                  prime [j] = 1;
20              } 
21              //把 0 -- n 内的素数存到数组中 
22               primeElem[cn ++] = i;
23         }
24     }
25}
26
27int main ()
28{
29    int m, n;
30    is_prime( );
31   
32    while ( scanf ("%d", &m) != EOF )
33    {
34           for ( int i = 0; i < m; i ++)
35          {
36              scanf ("%d", &n); 
37              
38              //找到 n 以内的 素因子
39              int product = 1;
40              int flag = 0;
41              for ( int j = 0; j < cn; j ++)
42              {
43                  flag = (int) sqrt (n * 1.0) + 1;   //n 的 方根 
44                  if ( primeElem[j] > flag )  //p1、p2、p3  这些素数要 <  n 
45                  break;
46                  int e = 0;
47                  while ( n % primeElem[j] == 0 )
48                  {
49                        e ++;
50                        n /= primeElem[j];
51                  }
52                  product *= ( 1 + 2 * e);
53              } 
54              if ( n > 1 )
55              product *= 3;  //这步要注意，最后有可能剩下的n本身就是原来n的素因子  ????
56              printf ("Scenario #%d:\n", i+1);
57              printf ("%d\n", ( product + 1 )/2 );
58             
59              printf ("\n");
60          }
61    }
62      //system ("pause");
63      return 0; 
64}
65

雪黛依梦 2010-08-05 15:18 发表评论