Inverse of Euler totient function

→ Обратите внимание

До соревнования
CodeTON Round 9 (Div. 1 + Div. 2, Rated, Prizes!)
26:56:18
Зарегистрироваться »

*есть доп. регистрация

→ Лидеры (рейтинг)

№	Пользователь	Рейтинг
1	tourist	4009
2	jiangly	3823
3	Benq	3738
4	Radewoosh	3633
5	jqdai0815	3620
6	orzdevinwang	3529
7	ecnerwala	3446
8	Um_nik	3396
9	ksun48	3390
10	gamegame	3386

Страны | Города | Организации

Всё →

→ Лидеры (вклад)

№	Пользователь	Вклад
1	cry	167
2	Um_nik	163
3	maomao90	162
3	atcoder_official	162
5	adamant	159
6	-is-this-fft-	158
7	awoo	157
8	TheScrasse	154
9	Dominater069	153
9	nor	153

Всё →

→ Найти пользователя

→ Прямой эфир

Детальнее →

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

Inverse of Euler totient function

Автор SLIMSHADYNICK, история, 5 лет назад, По-английски

There's this problem INVPHI on spoj, I have not been able to solve it. I have been getting NZEC runtime error. Could anyone tell me how to approach as I feel I have exceeded the memory limitations.

SLIMSHADYNICK
5 лет назад
4

Комментарии (4)

Написать комментарий?

feodorv

5 лет назад, # |

← Rev. 2 →

Here is the strightforward solution. We can sieve the phi values up to 202918035 (please, see comments to the problem statement).

Code, 9.35sec, 1157M

#include <stdio.h>

#define MAXN 100000000
#define MAXM 202918035

int phi[MAXM+1], minv[MAXN+1];

void gen( void )
{
  int i, j;
  for( i = 1; i <= MAXM; ++i) phi[i] = i;
  for( i = 2; i <= MAXM; ++i)
    if( phi[i] == i )
      for( j = i; j <= MAXM; j += i)
        phi[j] = (phi[j] / i) * (i - 1);
  for( i = 1; i <= MAXM; ++i)
    if( phi[i] <= MAXN && minv[phi[i]] == 0 )
      minv[phi[i]] = i;
}

int main( void )
{
  int t;
  gen();
  scanf( "%d", &t);
  while( t-- > 0 )
  {
    int n;
    scanf( "%d", &n);
    printf( "%d\n", minv[n] ? minv[n] : -1);
  }
  return 0;
}

Also you can note that almost all resulting numbers are odd with three exceptions, so we can sieve only odd values:

Code, 4.65sec, 770M

#include <stdio.h>

#define MAXN 100000000
#define MAXM 202918035

int phi[MAXM/2+1], minv[MAXN+1];

void gen( void )
{
  int i, j;
  for( i = 1; i <= MAXM; i += 2) phi[i>>1] = i;
  for( i = 3; i <= MAXM; i += 2)
    if( phi[i>>1] == i )
      for( j = i; j <= MAXM; j += 2*i)
        phi[j>>1] = (phi[j>>1] / i) * (i - 1);
  for( i = 1; i <= MAXM; i += 2)
    if( phi[i>>1] <= MAXN && minv[phi[i>>1]] == 0 )
      minv[phi[i>>1]] = i;
  minv[16842752] = 33817088;
  minv[33685504] = 67634176;
  minv[67371008] = 135268352;
}

int main( void )
{
  int t;
  gen();
  scanf( "%d", &t);
  while( t-- > 0 )
  {
    int n;
    scanf( "%d", &n);
    printf( "%d\n", minv[n] ? minv[n] : -1);
  }
  return 0;
}

Though I do not know how can we achieve 0.04 sec or how will it run in Java.

→ Ответить