Time Limit: 1000MS | Memory Limit: 10000K | |
Total Submissions: 39115 | Accepted: 12294 |
Description
for all possible inputs.
Consider the following algorithm:
1. input n 2. print n 3. if n = 1 then STOP 4. if n is odd then n <-- 3n+1 5. else n <-- n/2 6. GOTO 2
Given the input 22, the following sequence of numbers will be printed 22 11 34 17 52 26 13 40 20 10 5 16 8 4 2 1
It is conjectured that the algorithm above will terminate (when a 1 is printed) for any integral input value. Despite the simplicity of the algorithm, it is unknown whether this conjecture is true. It has been verified, however, for all integers n such that
0 < n < 1,000,000 (and, in fact, for many more numbers than this.)
Given an input n, it is possible to determine the number of numbers printed before the 1 is printed. For a given n this is called the cycle-length of n. In the example above, the cycle length of 22 is 16.
For any two numbers i and j you are to determine the maximum cycle length over all numbers between i and j.
Input
You should process all pairs of integers and for each pair determine the maximum cycle length over all integers between and including i and j.
Output
line and with one line of output for each line of input. The integers i and j must appear in the output in the same order in which they appeared in the input and should be followed by the maximum cycle length (on the same line).
Sample Input
1 10 100 200 201 210 900 1000
Sample Output
1 10 20 100 200 125 201 210 89 900 1000 174
Source
#include<stdio.h> #include<string.h> #define maxn 10000 struct node { __int64 left; __int64 right; __int64 max; }; __int64 dp[maxn+5]; node tree[3*maxn]; void dfs(__int64 n,__int64 num,__int64 id) { if(n<maxn&&dp[n]!=-1) { dp[id]=dp[n]+num; return; } else { if(n&1) { n=n*3+1; dfs(n,num+1,id); } else { n=n/2; dfs(n,num+1,id); } } } __int64 max(__int64 a,__int64 b) { if(a>b) return a; else return b; } void build(__int64 left,__int64 right,__int64 root) { tree[root].left=left; tree[root].right=right; if(left==right) { tree[root].max=dp[left]; return; } __int64 mid=(left+right)>>1; build(left,mid,root*2); build(mid+1,right,root*2+1); tree[root].max=max(tree[2*root].max,tree[2*root+1].max); } void init() { __int64 i; memset(dp,-1,sizeof(dp)); dp[1]=1; for(i=1;i<=maxn;i++) dfs(i,0,i); build(1,maxn,1); } void swap(__int64 &a,__int64 &b) { __int64 temp; temp=a; a=b; b=temp; } __int64 search(__int64 left,__int64 right,__int64 root) { if(tree[root].left==left&&tree[root].right==right) return tree[root].max; __int64 mid=(tree[root].left+tree[root].right)>>1; if(right<=mid) return search(left,right,root*2); else if(left>mid) return search(left,right,root*2+1); else return max(search(left,mid,root*2),search(mid+1,right,root*2+1)); } int main() { init(); __int64 a,b; while(scanf("%I64d%I64d",&a,&b)!=EOF) { if(a>b) { swap(a,b); printf("%I64d %I64d %I64d\n",b,a,search(a,b,1)); } else { printf("%I64d %I64d %I64d\n",a,b,search(a,b,1)); } } return 0; }