3295 - Perm Euler

Cerinţa

Indicatorul lui Euler, φ(n) – câteodată numit funcția phi, e folosit pentru a determina câte numere pozitive mai mici decât n care sunt relativ prime cu n există. De exemplu, cum 1, 2,4, 5, 7 și 8 sunt toate mai mici decât 9 și sunt relativ prime la 9, φ(9)=6. Numărul 1 e considerat a fi relativ prim cu toate numerele naturale, deci φ(1)=1. În mod interesant, φ(87109)=79180, și se poate observa că 87109 e o permutare a lui 79180.

Se consideră un șir de cel mult 10000 de numere naturale distincte mai mici decât 10.000.000. Să se scrie un program care găsește valoarea lui n, pentru care φ(n) e o permutare a lui n și fracția n/φ(n) are valoare minimă. Dacă sunt mai multe valori cu aceeași proprietate atunci se scrie prima valoare din șir. Dacă nu sunt valori cu proprietatea menționată se va scrie valoarea 0.

Date de intrare

Fișierul de intrare permeuler.in conține pe prima linie cel mult 10000 de numere naturale distincte mai mici decât 10.000.000, separate prin spații.

Date de ieşire

Dacă datele sunt introduse corect, în fișier se va afișa:"Datele sunt introduse corect.", apoi pe un rând nou fișierul de ieșire permeuler.out va conține primul număr x pentru care φ(x) e o permutare a lui x și fracția x/φ(x) are valoare minimă. În cazul în care datele nu respectă restricțiile, se va afișa în fișier:"Datele nu corespund restricțiilor impuse.".

Restricții și precizări

• în fișier sunt cel mult 10000 de numere
• numerele de pe prima linie a fișierului de intrare vor fi mai mici decât 10.000.000

Exemplu

permeuler.in

1000 345 21 34567 5678 63 345678 1234

permeuler.out

Datele sunt introduse corect.

21

Explicație

În fișierul de intrare sunt 2 numere care îndeplinesc proprietatea cerută de problemă: 21 care are φ(21)=12 și raportul 21/φ(21)=1.75 și 63 care are φ(63)=36 și raportul 63/φ(63)=1.75. În fișierul de ieșire va fi scrisă valoarea 21.

Rezolvare

#Funția de verificare a datelor de intrare 
def validare_date(numere):
    if len(numere) > 10000:
        return False
    if not all(0 < n < 10000000 for n in numere):
        return False
    if len(set(numere)) != len(numere):
        return False
    return True
#Funcția calculează indicatorul lui Euler pentru n
def phi(n):
    rezultat = n
    i = 2

    while i * i <= n:
        if n % i == 0:
            while n % i == 0:
                n //= i
            rezultat -= rezultat // i
        i += 1

    if n > 1:
        rezultat -= rezultat // n

    return rezultat

#Funcția verifică dacă numerele a este o permutare a lui b și invers
def este_permutare(a, b):
    dict_a = {c: str(a).count(str(c)) for c in range(10)}
    dict_b = {c: str(b).count(str(c)) for c in range(10)}
    return dict_a == dict_b


if __name__ == '__main__':
    with open('permeuler.in', 'r') as f:
        numere = list(map(int, f.readline().split()))

        # Verificăm dacă datele de intrare sunt valide
        if not validare_date(numere):
            f.write("Datele introduse nu  corespund restrictiilor impuse.")
            exit(1)

    minimul_raportului = float('inf')
    rezultat = 0
    # Parcurgem fiecare număr din lista și calculăm indicatorul Euler pentru acesta
    for n in numere:
        phin = phi(n)
        raport = n / phin
        # Verifică dacă n și indicator_n sunt permutări una a celeilalte și dacă raportul este minim
        if este_permutare(n, phin) and raport < minimul_raportului:
            minimul_raportului = raport
            rezultat = n
    #Deschidem fișierul și afișăm mesajul de confirmare si rezultatul
    with open('permeuler.out', 'w') as f:
        f.write("Datele sunt introduse corect.\n")
        f.write(str(rezultat))