結果
| 問題 | 
                            No.1346 Rectangle
                             | 
                    
| コンテスト | |
| ユーザー | 
                             lam6er
                         | 
                    
| 提出日時 | 2025-04-15 21:00:34 | 
| 言語 | PyPy3  (7.3.15)  | 
                    
| 結果 | 
                             
                                WA
                                 
                             
                            
                         | 
                    
| 実行時間 | - | 
| コード長 | 1,488 bytes | 
| コンパイル時間 | 323 ms | 
| コンパイル使用メモリ | 82,104 KB | 
| 実行使用メモリ | 59,576 KB | 
| 最終ジャッジ日時 | 2025-04-15 21:05:44 | 
| 合計ジャッジ時間 | 1,678 ms | 
| 
                            ジャッジサーバーID (参考情報)  | 
                        judge4 / judge5 | 
(要ログイン)
| ファイルパターン | 結果 | 
|---|---|
| sample | AC * 1 | 
| other | AC * 1 WA * 16 | 
ソースコード
import math
def is_prime(n):
    if n < 2:
        return False
    for p in [2, 3, 5, 7, 11, 13, 17, 19, 23, 29, 31, 37]:
        if n % p == 0:
            return n == p
    d = n - 1
    s = 0
    while d % 2 == 0:
        d //= 2
        s += 1
    for a in [2, 325, 9375, 28178, 450775, 9780504, 1795265022]:
        if a >= n:
            continue
        x = pow(a, d, n)
        if x == 1 or x == n - 1:
            continue
        for _ in range(s - 1):
            x = pow(x, 2, n)
            if x == n - 1:
                break
        else:
            return False
    return True
def is_prime_power(n):
    if n == 1:
        return False
    if is_prime(n):
        return True
    for p in [2, 3, 5, 7]:
        if n % p == 0:
            while n % p == 0:
                n //= p
            return n == 1
    max_p = int(math.isqrt(n)) + 1
    for p in range(11, max_p, 2):
        if p > n:
            break
        if n % p == 0:
            if not is_prime(p):
                return False
            while n % p == 0:
                n //= p
            return n == 1
    return False
def main():
    N = int(input().strip())
    if N == 2:
        print("INF")
        return
    M = N - 1
    if is_prime_power(M):
        print(M)
    else:
        # Find largest k such that k(k-1) <= N-1
        k = int((1 + math.isqrt(1 + 4 * (N - 1))) // 2)
        while k * (k - 1) > N - 1:
            k -= 1
        print(k)
if __name__ == "__main__":
    main()
            
            
            
        
            
lam6er