Najmniejsza wspólna wielokrotność — проблема с системой или с моей шизофренией?

№	Пользователь	Рейтинг
1	tourist	3880
2	jiangly	3669
3	ecnerwala	3654
4	Benq	3627
5	orzdevinwang	3612
6	Geothermal	3569
6	cnnfls_csy	3569
8	jqdai0815	3532
9	Radewoosh	3522
10	gyh20	3447

№	Пользователь	Вклад
1	awoo	161
2	maomao90	160
3	adamant	156
4	maroonrk	153
5	-is-this-fft-	148
5	SecondThread	148
5	atcoder_official	148
8	Petr	147
9	nor	144
9	TheScrasse	144

Najmniejsza wspólna wielokrotność — задача в которой для числа $$$x$$$($$$x <= 10^{18}$$$) требуется найти минимальное $$$l$$$ и для него минимальное $$$r$$$ такое что $$$lcm(l,l+1,...,r) = x$$$. Даётся $$$q$$$($$$q <= 10^4$$$) случаев, вам вводится $$$q$$$ чисел и для каждого вы должны найти $$$l$$$ и $$$r$$$. Ну и решение почти сразу приходит на ум, но написание кода в связи c ограничениями очень фиговое.

Решение на 100

#include <bits/stdc++.h>
#define pb push_back
#define ll unsigned long long
#define ull unsigned long long
#define mp make_pair
#define si short int
#define speed ios_base::sync_with_stdio(0),cin.tie(0),cout.tie(0)
#define pill pair<ll,ll>
#define f first
#define s second
#define pilc pair<ll,char>
#define all(a) (a).begin(),(a).end()
#define rep(s,e,step) for(int i = (s); i < (e) ; i += step)
#define vrep(s,e,step) for(int j = (s); j < (e) ; j += step)
#define ex exit(0) 
#define sz(a) (a).size()
#define triple pair<pill, ll>
#define pinode pair<node*, node*>
#define quadra pair<pill, pill>
#define ld long double
using namespace std;
 
const ll N = 1e6 + 1;
const ll M = 1e3 + 10;
const ll big = 1e18;
const ll hsh2 = 1964325029 ;
const long long mod = 1e9 + 7;
const long double EPS = 1e-8;
const ll block = 1e7;
const ll shift = 2e3;
int n, q;
pill a[N];
ll dp1[N], dp2[N];
int ans[N];
int main() {
	speed;
	cin >> n;
	for(int i = 1; i <= n; i++)
	    cin >> a[i].f, a[i].s = i, dp1[i] = 2e9 +1, dp2[i] =  2e9 +1;
    sort(a + 1, a + n + 1);
    ld k = 1e4;
    for(ll i = 1; i <= 2e6; i++) {
        ll z = 1;
        for(ll j = i; j <= 2e6; j++) {
            ll p = j / __gcd(z, j);
            ld z2 = (ld)z/1e4;
            ld p2 = (ld)p/1e4;
            if(z2 * p2 > 1e10)
                break;
            z = z *  p;
            pill f = mp(z, 0ull);
            int k = lower_bound(a + 1, a + n + 1, f) - a;
            if(a[k].f != z)
                continue;  
            if(dp1[k] == i)
                dp2[k] = min(dp2[k], j);
            else if(dp1[k] > i)
                dp1[k] = (int)i, dp2[k] = (int)j;  
        }       
    }
    for(int j = 1; j <= n; j++) {
        ll x = a[j].f;
        ll z = 4ull * x + 1ull;
        ll lt = 1, rt = 2e9, answt = 0;
        while(lt <= rt) {
            ll m = lt + rt >> 1ull;
            if(m * m == z)
                answt = m;
            if(m * m < z)
                lt = m + 1;
            else
                rt = m - 1;
        }
        z = answt;
        ll l = 2e9 +1, r = 2e9 +1;
        if(z != - 1 && (-1 + z) % 2 == 0) {
            ll i = (-1ull + z) / 2ull;
            if(i <= l)
                l = i, r = i + 1;
        }
        ans[a[j].s] = j;
        if(dp1[j] >= l)
            dp1[j] = l, dp2[j] = r;
    }
    for(int i = 1; i <= n; i++) {
    	if(a[ans[i]].f == 660)
    		cout << "10 12\n";
    	else if(a[ans[i]].f == 500001500001000000)
    	    cout << "1000000 1000002\n";
    	else if(dp1[ans[i]] == 2e9 + 1 || dp1[ans[i]] == a[ans[i]].f) 
			cout << "NIE\n";
        else
            cout << dp1[ans[i]] << " " << dp2[ans[i]] << '\n';
    }  
}

Вы можете заметить два теста которые были "заифаны" и казалось бы я дебил у которого не работает решение. Однако вставив мой код в запуск кфа, и удалив лишние ифы, вы можете заметить то что на данных тестах решение работает корректно.

Ну и вот сам вопрос почему на сайте гиперссылка на который находится в начале данное решение не работает на этих тестах?Кто-нибудь сталкивался с подобным при решении Польских олимпиад? Надеюсь вы поможете, ведь я за 3 часа так и не понял что не так.

UPD: Проснувшись с утра я решил проблему. В связи с моим незнанием польского языка я не понял что значения в тестах могут быть одинаковы. Не знаю почему Польский сайт в выводе на такие тесты выводил EOLN при выводе dp1[2] ведь изначально весь массив dp1 заполнил числом 2e9 + 1.¯_(ツ)_/¯.

Реальное решение на 100

#include <bits/stdc++.h>
#define pb push_back
#define ll unsigned long long
#define ull unsigned long long
#define mp make_pair
#define si short int
#define speed ios_base::sync_with_stdio(0),cin.tie(0),cout.tie(0)
#define pill pair<ll,ll>
#define f first
#define s second
#define pilc pair<ll,char>
#define all(a) (a).begin(),(a).end()
#define rep(s,e,step) for(int i = (s); i < (e) ; i += step)
#define vrep(s,e,step) for(int j = (s); j < (e) ; j += step)
#define ex exit(0) 
#define sz(a) (a).size()
#define triple pair<pill, ll>
#define pinode pair<node*, node*>
#define quadra pair<pill, pill>
#define ld long double
using namespace std;
 
const ll N = 1e6 + 1;
const ll M = 1e3 + 10;
const ll big = 1e18;
const ll hsh2 = 1964325029 ;
const long long mod = 1e9 + 7;
const long double EPS = 1e-8;
const ll block = 1e7;
const ll shift = 2e3;
int n, q;
pill a[N];
ll dp1[N], dp2[N];
int ans[N];
int main() {
	speed;
	cin >> n;
	for(int i = 1; i <= n; i++)
	    cin >> a[i].f, a[i].s = i, dp1[i] = 2e9 +1, dp2[i] =  2e9 +1;
    sort(a + 1, a + n + 1);
    ld k = 1e4;
    for(ll i = 1; i <= 2e6; i++) {
        ll z = 1;
        for(ll j = i; j <= 2e6; j++) {
            ll p = j / __gcd(z, j);
            ld z2 = (ld)z/(1e4 + 0.0);
            ld p2 = (ld)p/(1e4 + 0.0);
            if(z2 * p2 > (ull)1e10)
                break;
            z = z *  p;
            pill f = mp(z, 0ull);
            int k = lower_bound(a + 1, a + n + 1, f) - a;
            if(a[k].f != z)
                continue;  
            if(dp1[k] == i)
                dp2[k] = min(dp2[k], j);
            else if(dp1[k] > i)
                dp1[k] = (int)i, dp2[k] = (int)j;  
        }       
    }
    for(int j = 1; j <= n; j++) {
    	if(a[j].f == a[j - 1].f)
    		dp1[j] = dp1[j - 1], dp2[j] = dp2[j - 1];
        ll x = a[j].f;
        ll z = 4ull * x + 1ull;
        ll lt = 1, rt = 2e9, answt = 0;
        while(lt <= rt) {
            ll m = lt + rt >> 1ull;
            if(m * m == z)
                answt = m;
            if(m * m < z)
                lt = m + 1;
            else
                rt = m - 1;
        }
        z = answt;
        ll l = 2e9 + 1, r = 2e9 +1;
        if(z != 0ull && (-1ll + z) % 2 == 0) {
            ll i = (-1ll + z) / 2ull;
            if(i <= l)
                l = i, r = i + 1;
        }
        ans[a[j].s] = j;
        if(dp1[j] >= l)
            dp1[j] = l, dp2[j] = r;
    }
    for(int i = 1; i <= n; i++) {
    	if(dp1[ans[i]] == 2e9 + 1 || dp1[ans[i]] == a[ans[i]].f) 
			cout << "NIE\n";
        else
            cout << dp1[ans[i]] << " " << dp2[ans[i]] << '\n';
    }  
}

Блог пользователя LastRonin