Последовательность Кеане Опции

[Unconnected]

Привет всем

Пробую решить олимпиадную задачку (конечно, с уже давно окончившейся олимпиады )), про последовательность (битов) Кеане - кто этот Кеане, я не знаю, и Вики с гуглом тоже. Формулировка такая:

Последовательность Кеане
(Время: 1 сек. Память: 16 Мб Сложность: 57%)

Бесконечная последовательность битов, предложенная Кеане, равна 001001110001001110110110001… и формируется следующим алгоритмом: вначале записывается 0, потом 001, далее 001001110, то есть, для получения следующего члена, предыдущий записывается дважды, а справа приписывается его отрицание. Элементы этого ряда являются начальными подпоследовательностями Кеане.

Требуется написать программу, которая по заданному n определит N-й бит этой последовательности.
Входные данные

Входной файл INPUT.TXT содержит число N (N <= 10²⁰⁰).
Выходные данные

В выходной файл OUTPUT.TXT должен содержать найденный бит.

Просто генерировать последовательность у меня получается, вот так:

{$APPTYPE CONSOLE}

var n:int64;
    s:string;

function invers(s:string):string;
var i:integer;
begin
  for i:=1 to length(s) do if s[i]='0' then result:=result+'1' else result:=result+'0';
end;

begin
  assign(input,'input.txt');
  readln(n);
  s:='0';
  while (length(s)<=n) do s:=s+s+invers(s);
  assign(output,'output.txt');
  writeln(s[n]);
end.

(хотя, у меня такое ощущение, что она как-то неправильно генерирует)

Проблема в очень большом максимально возможном N. Здесь, видимо нужно использовать длинную арифметику, но я вот не понимаю, в каком месте её использовать. Ну допустим есть строка-последовательность, есть N в виде строки, и что с чем складывать? Кажется, даже строку с последовательностью такой длины получить не получится, и тем более length() от неё. Может, формулу можно вывести для нахождения бита по позиции?

Сообщение отредактировано: Unconnected - 5.06.2010 17:22

[volvo]

Unconnected, вот смотри, какой путь я бы избрал (показываю с применением обычной арифметики, с небольшими значениями). Идея - такая.

1. определяем, в каком члене последовательности находится искомый бит. То есть, допустим, надо найти значение бита под номером 20. Вот так выглядит последовательность, сформированная по приведенному тобой правилу:
   0_001_001001110_001001110001001110110110001 (40 бит). Она состоит, как видим, из 4-х членов (я разделил их символами подчеркивания). Нужный нам бит находится в третьем члене, начиная с 0-го;
2. находим номер нужного бита, начиная с начала того члена, в котором он находится;
3. последовательно делим текущий член на 3 части, и определяем, в какой части находится нужный бит (в зависимости от этого принимаем решение о необходимости инвертирования);
4. продолжаем эти действия до тех пор, пока не дойдем до 3-х битных кусочков. Тут уже можно с легкостью определить, какой бит чему равен.

Все действия прокомментирую кодом (у меня для проверки строится и сама последовательность, и результаты выводятся как вычисленные, так и взятые из строки. Потестируй, если какой-то результат не совпадет - говори, при каком N, будем думать что делать)

const
  inverted: array['0' .. '1'] of char = ('1', '0');

function log3(n: integer): real;
begin
  result := ln(n) / ln(3);
end;
function power(n: integer): longint;
var
  i: integer;
begin
  result := 1;
  for i := 1 to n do result := result * 3;
end;

function generate(len: longint): string;

  function invert(s: string): string;
  var i: integer;
  begin
    result := '';
    for i := 1 to length(s) do result := result + inverted[s[i]];
  end;

  var s: string;

begin
  s := '0';
  result := '';
  repeat
    result := result + s;
    s := s + s + invert(s);
    writeln(s);
  until length(result) > len;
end;

var
  s: string;
  n: integer;

  k, bit_in_seq: longint;
  power_3: longint;
  negate: boolean;

  res: char;

begin
  n := 22;
  s := generate(n);

  writeln(length(s), ' -> ', s);
  writeln('from string: ', s[n]);

  k := trunc(log3(2 * n + 1));
  bit_in_seq := pred(n - ((power(k) - 1) div 2));
  negate := false;
  power_3 := power(k);

  if power_3 > 3 then
    repeat
      power_3 := power_3 div 3;
      if (bit_in_seq div power_3) = 2 then negate := not negate;
      bit_in_seq := bit_in_seq mod power_3;
    until power_3 = 3;

  case bit_in_seq of
    0, 1: res := '0';
       2: res := '1';
  end;
  if negate then res := inverted[res];

  writeln('calculated: ', res);
  readln;
end.

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

P.S. И все-таки, по поводу

кто этот Кеане, я не знаю, и Вики с гуглом тоже

- это Mephisto Waltz Sequence