1. 程式人生 > >POJ 1006 - 生理周期 題解

POJ 1006 - 生理周期 題解

cli tdi aps 情感 blog spa etc play 轉載

此文為博主原創題解,轉載時請通知博主,並把原文鏈接放在正文醒目位置。

題目鏈接:http://poj.org/problem?id=1006&lang=zh-CN&change=true

Description

人生來就有三個生理周期,分別為體力、感情和智力周期,它們的周期長度為23天、28天和33天。每一個周期中有一天是高峰。在高峰這天,人會在相應的方面表現出色。例如,智力周期的高峰,人會思維敏捷,精力容易高度集中。因為三個周期的周長不同,所以通常三個周期的高峰不會落在同一天。對於每個人,我們想知道何時三個高峰落在同一天。對於每個周期,我們會給出從當前年份的第一天開始,到出現高峰的天數(不一定是第一次高峰出現的時間)。你的任務是給定一個從當年第一天開始數的天數,輸出從給定時間開始(不包括給定時間)下一次三個高峰落在同一天的時間(距給定時間的天數)。例如:給定時間為10,下次出現三個高峰同天的時間是12,則輸出2(註意這裏不是3)。

Input

輸入四個整數:p, e, i和d。 p, e, i分別表示體力、情感和智力高峰出現的時間(時間從當年的第一天開始計算)。d 是給定的時間,可能小於p, e, 或 i。 所有給定時間是非負的並且小於365, 所求的時間小於21252。

當p = e = i = d = -1時,輸入數據結束。

Output

從給定時間起,下一次三個高峰同天的時間(距離給定時間的天數)。

采用以下格式:
Case 1: the next triple peak occurs in 1234 days.

註意:即使結果是1天,也使用復數形式“days”。

Sample Input

0 0 0 0
0 0 0 100
5 20 34 325
4 5 6 7
283 102 23 320
203 301 203 40
-1 -1 -1 -1

Sample Output

Case 1: the next triple peak occurs in 21252 days.
Case 2: the next triple peak occurs in 21152 days.
Case 3: the next triple peak occurs in 19575 days.
Case 4: the next triple peak occurs in 16994 days.
Case 5: the next triple peak occurs in 8910 days.
Case 6: the next triple peak occurs in 10789 days.


分析:
這是中國剩余定理的典型例題,然而窩的CRT死活也跑不出答案來QAQ傷心
一怒之下寫了暴力,AC。

後來又改了一個多小時,總算是把CRT的代碼寫對了qwq...

AC代碼:
技術分享
 1 #include<cstdio>
 2 #include<algorithm>
 3 #include<cmath>
 4 #include<cstring>
 5 
 6 const int MAXN = 1000005;
 7 
 8 inline void read(int &x)
 9 {
10     char ch = getchar(),c = ch;x = 0;
11     while(ch < 0 || ch > 9) c = ch,ch = getchar();
12     while(ch >= 0 && ch <= 9) x = (x<<1)+(x<<3)+ch-0,ch = getchar();
13     if(c == -) x = -x;
14 }
15 
16 int main()
17 {
18     int a,b,c,d,cnt = 0,lim = 21252;
19     read(a),read(b),read(c),read(d);
20     while(d != -1)
21     {
22         ++cnt,a %= 23,b %= 28,c %= 33;
23         for(register int i = 1;;++ i)
24             if(i%23 == a && i%28 == b && i%33 == c)
25             {
26                 if(i < d) i += lim;
27                 printf("Case %d: the next triple peak occurs in %d days.\n",cnt,i-d);
28                 break;
29             }
30         read(a),read(b),read(c),read(d);
31     }
32     return 0;
33 }
暴力做法

技術分享
 1 #include<cstdio>
 2 #include<algorithm>
 3 #include<cmath>
 4 #include<cstring>
 5 
 6 const int MAXN = 1000005;
 7 
 8 inline void read(int &x)
 9 {
10     char ch = getchar(),c = ch;x = 0;
11     while(ch < 0 || ch > 9) c = ch,ch = getchar();
12     while(ch >= 0 && ch <= 9) x = (x<<1)+(x<<3)+ch-0,ch = getchar();
13     if(c == -) x = -x;
14 }
15 
16 int f[5] = {0,23,28,33},c[5],inv[5];
17 //f存儲周期,c存儲余數(也就是一開始讀入的),inv存儲逆元 
18 int cnt,ans,lim = 21252; 
19 
20 void exgcd(int a,int b,int &x,int &y)
21 {
22     if(!b)
23     {
24         x = 1,y = 0;
25         return ;
26     }
27     exgcd(b,a%b,y,x);
28     y -= x*(a/b);
29 }
30 
31 inline void init()
32 {//事先把逆元求出來
33     int x,y,b;
34     for(int i = 1;i <= 3;++ i)
35     {
36         x = 0,y = 0,b = lim/f[i];
37         exgcd(b,f[i],x,y);
38         inv[i] = x;
39     }
40 }
41 
42 int main()
43 {
44     init();
45     for(int i  = 1;i <= 3;++ i)
46         read(c[i]),c[i]%=f[i];
47     read(c[4]);
48     while(c[4] != -1)
49     {
50         ++ cnt;
51         ans = 0;
52         for(int i = 1;i <= 3;++ i)
53         //用CRT公式求解ans值 
54             ans = (ans+(lim/f[i])*c[i]*inv[i])%lim;
55         if(ans <= c[4]) ans += lim;
56         printf("Case %d: the next triple peak occurs in %d days.\n",cnt,ans-c[4]);
57         for(int i  = 1;i <= 3;++ i)
58             read(c[i]),c[i]%=f[i];
59         read(c[4]);
60     }
61     return 0;
62 }
CRT

POJ 1006 - 生理周期 題解