Finding fraction. - Codeforces

9 лет назад, скрыть # ^ |

← Rev. 2 →

I dont really get why this gives us a fraction with smallest denumerator. Could you please explain it more?

→ Ответить

farmersrice

9 лет назад, скрыть # ^ |

← Rev. 3 →

-10

Oops, I misread. I think this algorithm will give you smallest denominator:

Binary search on the denominator (it should never be the case that a smaller denominator works while a larger one does not, if so please tell). For every denominator, also binary search on the numerator, to try to get num/denom within [l, r). Then you can find the smallest denominator.

EDIT: Just kidding, it doesn't work

→ Ответить

9 лет назад, скрыть # ^ |

I was thinking of binary search, but is there any proof that if lets say for d1 > d2 such that for denumerator d1 is not possible, but for d2 is possible? if there is a proof, then binary search is correct. Could you please devise a proof, if you have in mind?

→ Ответить

Saat

9 лет назад, скрыть # ^ |

That is not correct. Think for example in range [0.4, 0.6]. In this case, denominator 2 works, while denominator 3 doesn't.

→ Ответить

9 лет назад, скрыть # ^ |

do you know how to solve this problem?

→ Ответить

9 лет назад, скрыть # ^ |

lets say, l = 0.115 and r = 0.125,

so we want to find a fraction which has smallest denumerator, and for this range it is 2/17. How to find 2/17 with your approach?

→ Ответить

9 лет назад, скрыть # ^ |

Is y — must smallest as possible???

→ Ответить

9 лет назад, скрыть # ^ |

Yes, in which it satisfies the condition.

→ Ответить

9 лет назад, скрыть # ^ |

← Rev. 2 →


//  l <= x/ y < r  --> so  l * y <= x < r * y

for( long long y = 1; true; ++y)
{
    double x_l = l * y ;  
    double x_r = r * y ;
    
    long long x = static_cast<long long>( x_l );
 if ( x < x_l - 1.0E-12 ) // if x strong less than x_l, increase it.
           x++;
    // now  x >= x_l
    // check  x < x_r
    if ( x < x_r - 1.0E-12) 
       { 
        printf("%d / %d\n", x,y); 
       return 0;
     }
     
}

see code

→ Ответить

9 лет назад, скрыть # ^ |

for for some y, it could take much more, since it is multiple test case problem.

→ Ответить

9 лет назад, скрыть # ^ |

← Rev. 2 →

y*r > 0 --> y*r >= 1 --> y >= 1.0 / r

long long y_min = static_cast< long long>( 1.0 / r );
for(long long y = y_min; true; ++y){ ...}

→ Ответить

9 лет назад, скрыть # ^ |

Note tha, binary search does not applicable there.

→ Ответить

http://acm.timus.ru/problem.aspx?space=1&num=1011 ))

9 лет назад, скрыть # |

→ Ответить

9 лет назад, скрыть # ^ |

Yes seems same problem. but constraints here are small. What is the solution for this problem?

→ Ответить

9 лет назад, скрыть # |

Auto comment: topic has been updated by Logarithmic (previous revision, new revision, compare).

→ Ответить

Saat

9 лет назад, скрыть # |

← Rev. 2 →

Ok, I'm not completely sure about this approach, but here it goes. Take u = 1 / l and v = 1 / r. This are numbers bigger than one. Then, the inequality that must be satisfied is $\text{[math]}$ . Then,iterate through x's. Binary search for each x the value of y, and for the minimum x such that y exists, that y is your answer. According to my intuition, (sorry, I cannot offer more than that), for most cases this approach will require not more than 10⁶ iterations or so.

UPD: Sorry, it doesn't work.

→ Ответить

farmersrice

9 лет назад, скрыть # ^ |

Can you explain how it's only 10⁶? Wouldn't you need to iterate over a lot more x-values, in case like [0.573529588, 0.573529590)?

→ Ответить

Saat

9 лет назад, скрыть # ^ |

Yeah, sorry, you are right.

→ Ответить

9 лет назад, скрыть # ^ |

But if there are somecases in which it takes about 10^8, then it will definitely fail, since it is multiple test case.

→ Ответить

dacin21

9 лет назад, скрыть # |

← Rev. 3 →

It seems like you could use the Stern–Brocot tree.

Start at the root. At each node, you either find the solution or have to recurse into one of the children.

If this is too slow you might be able to binary-search how far you walk down left/right before chaning to right/left to get O(log(10⁹)²) time complexity, as the denominator grows at least as fast as the fibonacci series when alternating left/right.

Edit: One might get O(log(10⁹)) complexity by using exponential-search instead of binary-search.

→ Ответить

9 лет назад, скрыть # ^ |

How to apply binary search here? Binary search on y seems not working :( Could you please explain it more detailed?

→ Ответить

mahmud2690

9 лет назад, скрыть # ^ |

Can you possibly explain your binary search approach. I didn't get your idea.

→ Ответить

dacin21

9 лет назад, скрыть # ^ |

← Rev. 10 →

I'll first explain the basic algorithm on the Stern-Brocot tree similar to here.

We start with xl = 0, yl = 1, xr = 1, yr = 1. Throughout the algorithm we keep $\text{[math]}$ .

To traverse the tree we look at the middle fraction $\text{[math]}$ . If $\text{[math]}$ , then $\text{[math]}$ is the solution (Correctness follows from the properties of the Stern-Brocot tree).

If $\text{[math]}$ , set xl = xm, yl = ym (This is a move to the right).

If $\text{[math]}$ , set xr = xm, xr = ym (This is a move to the left).

The above algorithm takes a long time if we keep moving the same direction over and over.

What we can do, is move k times to the left by setting $\text{[math]}$ and similarily for moving to the right.

The faster algorithm now works as follows: When moving into a certain direction, binary-search for the value of k — how far the slow algorithm would move into that direction before changing directions. Then move k steps into that direction.

Edit: Fixed some small things.

Edit2:

(Messy) Code

#define x first
#define y second
pair<long long, long long> solve(long double const&L, long double const&R){
    pair<long long, long long> l(0, 1), r(1, 1);
    if(L==0.0) return l; // corner case
    for(;;){
        pair<int, int> m(l.x+r.x, l.y+r.y);
        if(m.x<L*m.y){// move to the right;
            long long kl=1, kr=1;
            while(l.x+kr*r.x <= L*(l.y+kr*r.y))kr*=2;// exponential search
            while(kl!=kr){
                long long km = (kl+kr)/2;
                if(l.x+km*r.x < L*(l.y+km*r.y)) kl=km+1;
                else kr=km;
            }
            l = make_pair(l.x+(kl-1)*r.x, l.y+(kl-1)*r.y);
        } else if(m.x>=R*m.y){//move to the left
            long long kl=1, kr=1;
            while(r.x+kr*l.x>=R*(r.y+kr*l.y))kr*=2;// exponential search
            while(kl!=kr){
                long long km = (kl+kr)/2;
                if(r.x+km*l.x>=R*(r.y+km*l.y)) kl = km+1;
                else kr = km;
            }
            r = make_pair(r.x+(kl-1)*l.x, r.y+(kl-1)*l.y);
        } else {
            return m;
        }
    }
}

→ Ответить

9 лет назад, скрыть # |

← Rev. 2 →

+13

For a given y_max, let's ask how many valid fractions with y ≤ y_max there are. For chosen y, the number of valid x's is ⌊ ry⌋ - ⌊ ly⌋. So the total number of valid fractions is:

$\text{[math]}$

Now, note that we can express l as a fraction $\text{[math]}$ , where L is an integer. Same with r. It's possible to calculate the sum $\text{[math]}$ quickly using this algorithm. Once we can find it, we can do a binary search on that value, and the result will be the smallest y_max that produces a non-zero result.

→ Ответить

9 лет назад, скрыть # ^ |

Thanks for the solution. One more point, it should be strictly less than r. how to exclude that? is the only way just to replace r with some r-eps?

→ Ответить

9 лет назад, скрыть # ^ |

← Rev. 2 →

+10

I would subtract the number of x = ry from the result (and add the number of x = ly because the other inequality is ≤ ). If $\text{[math]}$ where p,q are coprime, then a valid (x, y) pair is in the form (kp, kq) with k integer, so there are $\text{[math]}$ such points.

→ Ответить

9 лет назад, скрыть # ^ |

Yes, it will work, thanks for the help

→ Ответить

9 лет назад, скрыть # ^ |

Let me show a test which your formula gives wrong answer:

l = 0.125 and r = 0.130 your formula gives 4 / 31 , but correct answer is 1/8.

→ Ответить

9 лет назад, скрыть # ^ |

← Rev. 2 →

+10

Without the clarification I made in my other comment above it is indeed incorrect. You need to add $\text{[math]}$ , where q_r and q_l are denominators of $\text{[math]}$ and $\text{[math]}$ respectively, after simplifying the fractions (q_r = 100 and q_l = 8 in your case).

→ Ответить

9 лет назад, скрыть # ^ |

Thx, now understand.

your link not worked for me, can share another link ?

→ Ответить

9 лет назад, скрыть # ^ |

+13

Did you try adding www or http to the beginning of the address? The outline of the idea is:

Let's take a sum $\text{[math]}$ with p and q coprime. If n ≥ q or p > q, they can be reduced to n mod q (since $\text{[math]}$ ) and p mod q (since $\text{[math]}$ ) by adding a constant. Now, let $\text{[math]}$ . After a few transformations we can get $\text{[math]}$ . Since p < q, we can again reduce q to q mod p and so on. The reductions structure will be the same as gcd calculation, and a single step will take constant time, giving overall O(log(p)+log(q)) complexity.

→ Ответить

sengi

7 лет назад, скрыть # ^ |

The last formula only holds when xq/p is never an integer for x from 1 to m. I think this can happen from time to time.

→ Ответить

dacin21

7 лет назад, скрыть # ^ |

+10

After the first two steps, we have that n < q, p ≤ q. Hence m < p. As p and q are coprime, p divides xq if and only if p divides x. But 0 < x ≤ m < p contradiction. Hence $\text{[math]}$ is never an integer.

(And the post you replied to was 2 years old.)

→ Ответить

sengi

7 лет назад, скрыть # ^ |

← Rev. 2 →

Yes, you are right. I was think of the situation when n is greater than q.

(aren't you curious why you get 10 thumbs up under a blog that is 2 years old?)

→ Ответить