#756701 (Crystal) No.5005 3-SAT

提出ソース

結果

問題	No.5005 3-SAT
ユーザー	tomerun
提出日時	2022-04-29 16:39:27
言語	Crystal (1.11.2)
結果	AC
実行時間	1,581 ms / 2,000 ms
コード長	7,429 bytes
コンパイル時間	21,965 ms
実行使用メモリ	5,748 KB
スコア	81,854
最終ジャッジ日時	2022-04-29 16:42:38
合計ジャッジ時間	184,685 ms
ジャッジサーバーID （参考情報）	judge12 / judge14
	純コード判定しない問題か言語

このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示

入力	結果	実行時間実行使用メモリ
testcase_00	AC	1,578 ms 5,652 KB
testcase_01	AC	1,578 ms 5,636 KB
testcase_02	AC	1,577 ms 5,624 KB
testcase_03	AC	1,578 ms 5,692 KB
testcase_04	AC	1,579 ms 5,588 KB
testcase_05	AC	1,579 ms 5,604 KB
testcase_06	AC	1,578 ms 5,632 KB
testcase_07	AC	1,579 ms 5,592 KB
testcase_08	AC	1,579 ms 5,596 KB
testcase_09	AC	1,579 ms 5,620 KB
testcase_10	AC	1,579 ms 5,696 KB
testcase_11	AC	1,579 ms 5,600 KB
testcase_12	AC	1,580 ms 5,724 KB
testcase_13	AC	1,578 ms 5,696 KB
testcase_14	AC	1,579 ms 5,632 KB
testcase_15	AC	1,580 ms 5,748 KB
testcase_16	AC	1,579 ms 5,640 KB
testcase_17	AC	1,579 ms 5,636 KB
testcase_18	AC	1,579 ms 5,716 KB
testcase_19	AC	1,579 ms 5,672 KB
testcase_20	AC	1,568 ms 5,624 KB
testcase_21	AC	1,579 ms 5,668 KB
testcase_22	AC	1,579 ms 5,612 KB
testcase_23	AC	1,579 ms 5,696 KB
testcase_24	AC	1,580 ms 5,600 KB
testcase_25	AC	1,580 ms 5,616 KB
testcase_26	AC	1,580 ms 5,632 KB
testcase_27	AC	1,580 ms 5,588 KB
testcase_28	AC	1,579 ms 5,720 KB
testcase_29	AC	1,579 ms 5,740 KB
testcase_30	AC	1,579 ms 5,748 KB
testcase_31	AC	1,579 ms 5,656 KB
testcase_32	AC	1,580 ms 5,604 KB
testcase_33	AC	1,580 ms 5,664 KB
testcase_34	AC	1,579 ms 5,596 KB
testcase_35	AC	1,579 ms 5,716 KB
testcase_36	AC	1,579 ms 5,668 KB
testcase_37	AC	1,580 ms 5,672 KB
testcase_38	AC	1,581 ms 5,644 KB
testcase_39	AC	1,578 ms 5,676 KB
testcase_40	AC	1,580 ms 5,588 KB
testcase_41	AC	1,578 ms 5,596 KB
testcase_42	AC	1,580 ms 5,708 KB
testcase_43	AC	1,578 ms 5,716 KB
testcase_44	AC	1,579 ms 5,600 KB
testcase_45	AC	1,578 ms 5,616 KB
testcase_46	AC	1,580 ms 5,696 KB
testcase_47	AC	1,580 ms 5,656 KB
testcase_48	AC	1,567 ms 5,600 KB
testcase_49	AC	1,580 ms 5,684 KB
testcase_50	AC	1,579 ms 5,664 KB
testcase_51	AC	1,580 ms 5,624 KB
testcase_52	AC	1,578 ms 5,592 KB
testcase_53	AC	1,579 ms 5,676 KB
testcase_54	AC	1,579 ms 5,748 KB
testcase_55	AC	1,579 ms 5,748 KB
testcase_56	AC	1,579 ms 5,588 KB
testcase_57	AC	1,580 ms 5,684 KB
testcase_58	AC	1,578 ms 5,604 KB
testcase_59	AC	1,579 ms 5,584 KB
testcase_60	AC	1,580 ms 5,736 KB
testcase_61	AC	1,579 ms 5,588 KB
testcase_62	AC	1,579 ms 5,596 KB
testcase_63	AC	1,578 ms 5,592 KB
testcase_64	AC	1,580 ms 5,672 KB
testcase_65	AC	1,580 ms 5,748 KB
testcase_66	AC	1,578 ms 5,708 KB
testcase_67	AC	1,579 ms 5,668 KB
testcase_68	AC	1,580 ms 5,608 KB
testcase_69	AC	1,580 ms 5,608 KB
testcase_70	AC	1,579 ms 5,620 KB
testcase_71	AC	1,579 ms 5,628 KB
testcase_72	AC	1,580 ms 5,692 KB
testcase_73	AC	1,578 ms 5,596 KB
testcase_74	AC	1,580 ms 5,656 KB
testcase_75	AC	1,580 ms 5,668 KB
testcase_76	AC	1,579 ms 5,604 KB
testcase_77	AC	1,579 ms 5,632 KB
testcase_78	AC	1,579 ms 5,668 KB
testcase_79	AC	1,579 ms 5,620 KB
testcase_80	AC	1,578 ms 5,608 KB
testcase_81	AC	1,580 ms 5,632 KB
testcase_82	AC	1,580 ms 5,676 KB
testcase_83	AC	1,579 ms 5,632 KB
testcase_84	AC	1,579 ms 5,588 KB
testcase_85	AC	1,579 ms 5,636 KB
testcase_86	AC	1,578 ms 5,696 KB
testcase_87	AC	1,578 ms 5,696 KB
testcase_88	AC	1,579 ms 5,628 KB
testcase_89	AC	1,579 ms 5,660 KB
testcase_90	AC	1,579 ms 5,664 KB
testcase_91	AC	1,580 ms 5,640 KB
testcase_92	AC	1,580 ms 5,664 KB
testcase_93	AC	1,578 ms 5,672 KB
testcase_94	AC	1,580 ms 5,716 KB
testcase_95	AC	1,578 ms 5,612 KB
testcase_96	AC	1,580 ms 5,600 KB
testcase_97	AC	1,580 ms 5,716 KB
testcase_98	AC	1,580 ms 5,604 KB
testcase_99	AC	1,579 ms 5,652 KB

権限があれば一括ダウンロードができます

ソースコード

raw source code

START_TIME = Time.utc.to_unix_ms
TL         = 1580
N          = 2048
M          =  256
INF        = 1 << 29
RND        = Random.new(2)

macro debug(msg)
  {% if flag?(:trace) %}
    STDERR.puts({{msg}})
  {% end %}
end

macro debugf(format_string, *args)
  {% if flag?(:trace) %}
    STDERR.printf({{format_string}}, {{*args}})
  {% end %}
end

class IndexSet
  getter :upper

  def initialize(@upper : Int32)
    @que = [] of Int32
    @pos = Array(Int32).new(@upper, -1)
  end

  def add(v)
    if @pos[v] == -1
      @pos[v] = @que.size
      @que << v
    end
  end

  def remove(v)
    p = @pos[v]
    b = @que[-1]
    @que[p] = b
    @que.pop
    @pos[b] = p
    @pos[v] = -1
  end

  def get_random
    @que[RND.rand(@que.size)]
  end

  def size
    @que.size
  end

  def to_s(io)
    @que.to_s(io)
  end
end

class Result
  property :bits, :score

  def initialize(@bits : Array(Int32), @score : Int32)
  end

  def to_s(io)
    io << @bits.map { |v| v == -1 ? 0 : v }.reverse.join
  end
end

struct Cond
  getter :var, :val

  def initialize(@var : Int32, @val : Int32)
  end
end

alias Term = Tuple(Cond, Cond, Cond)

class Solver
  def initialize
    @terms = Array(Term).new(N) do
      a, b, c, p, q, r = read_line.split.map(&.to_i)
      {Cond.new(a, p), Cond.new(b, q), Cond.new(c, r)}
    end
    @var_pos = Array(Array(Tuple(Int32, Int32))).new(M) { [] of Tuple(Int32, Int32) }
    N.times do |i|
      3.times do |j|
        term = @terms[i][j]
        @var_pos[term.var] << {i, term.val}
      end
    end
    @best_result = Result.new(Array.new(M, -1), 0)
    @ans = Array(Int32).new(M, -1)
    @search_order = [0, 1, 2]
    @dfs_count = 0
    @eval = Array(Int32).new(3, 0)
    @ok_count = Array(Int32).new(N, 0)
    @ng_count = Array(Int32).new(N, 0)
    @pass = Array(Bool).new(N) do |i|
      ts = @terms[i]
      [[0, 1], [1, 2], [2, 0]].any? { |ps| ts[ps[0]].var == ts[ps[1]].var && ts[ps[0]].var != ts[ps[1]].var }
    end
  end

  def solve(timelimit)
    start_time = Time.utc.to_unix_ms
    before_time = start_time
    solve_single()
    after_time = Time.utc.to_unix_ms
    worst_time = after_time - before_time
    turn = 0
    while true
      before_time = after_time
      if before_time + worst_time > timelimit
        debug("turn:#{turn} worst_time:#{worst_time}")
        break
      end
      solve_single()
      after_time = Time.utc.to_unix_ms
      worst_time = {worst_time, after_time - before_time}.max
      turn += 1
    end
    STDERR.puts("validate_score:#{score(@best_result.bits)}")
    return @best_result
  end

  def solve_single
    @ans.fill(-1)
    @dfs_count = 0
    @ok_count.fill(0)
    @ng_count.fill(0)
    pos = 0
    while pos < N
      if try_next(pos)
        while pos < N && @ok_count[pos] > 0
          pos += 1
        end
        if pos > @best_result.score
          @best_result.score = pos
          @best_result.bits[0, M] = @ans
          debug("best_score:#{pos}")
        end
        # if score(@ans) != pos
        #   debug("validate_score:#{score(@ans)} score:#{pos}")
        #   exit
        # end
      else
        break
      end
    end
  end

  def try_next(pos)
    shuffle(@search_order)
    3.times do |i|
      t = @terms[pos][@search_order[i]]
      if @ans[t.var] == -1
        @ans[t.var] = t.val
        @var_pos[t.var].each do |ps|
          p, v = ps
          if v == t.val
            @ok_count[p] += 1
          else
            @ng_count[p] += 1
          end
        end
        return true
      end
    end
    set = IndexSet.new(N)
    set.add(pos)
    used = Set(Int32).new
    @dfs_count = 0
    return try_dfs(pos, set, used)
  end

  def try_dfs(pos, set, used)
    # debug("try_dfs:#{pos} #{set} #{used}")
    if set.size == 0
      return true
    end
    if @dfs_count > 100
      return false
    end
    cur = set.get_random
    shuffle(@search_order)
    @dfs_count += 1
    3.times do |i|
      t = @terms[cur][@search_order[i]]
      next if used.includes?(t.var)
      used << t.var
      old = @ans[t.var]
      @ans[t.var] = t.val
      @var_pos[t.var].each do |ps|
        p, v = ps
        if v == t.val
          @ok_count[p] += 1
          @ng_count[p] -= 1
          if p <= pos && @ok_count[p] == 1
            set.remove(p)
          end
        elsif old != -1
          @ok_count[p] -= 1
          @ng_count[p] += 1
          if p <= pos && @ok_count[p] == 0
            set.add(p)
          end
        end
        # debug("@ok_count[#{p}]=#{@ok_count[p]}")
        # debug("@ng_count[#{p}]=#{@ng_count[p]}")
      end
      if try_dfs(pos, set, used)
        return true
      end
      used.delete(t.var)
      @ans[t.var] = old
      @var_pos[t.var].each do |ps|
        p, v = ps
        if v == t.val
          @ok_count[p] -= 1
          @ng_count[p] += 1
          if p <= pos && @ok_count[p] == 0
            set.add(p)
          end
        elsif old != -1
          @ok_count[p] += 1
          @ng_count[p] -= 1
          if p <= pos && @ok_count[p] == 1
            set.remove(p)
          end
        end
      end
    end
    # set.add(cur)
    return false
  end

  def solve_single2
    @ans.fill(-1)
    @dfs_count = 0
    @ok_count.fill(0)
    @ng_count.fill(0)
    dfs(0)
  end

  def dfs(pos)
    if @pass[pos]
      dfs(pos + 1)
    end
    if pos > @best_result.score
      @best_result.score = pos
      @best_result.bits[0, M] = @ans
      debug("best_score:#{pos}")
    end
    if @dfs_count > 1000
      return
    end
    if @terms[pos].any? { |t| @ans[t.var] == t.val }
      dfs(pos + 1)
    else
      3.times do |i|
        t = @terms[pos][i]
        if @ans[t.var] != -1
          @eval[i] = -INF
          next
        end
        @eval[i] = RND.rand(5000000)
        @var_pos[t.var].each do |ps|
          p, v = ps
          next if p <= pos
          if v == t.val
            if @ok_count[p] == 0
              # @eval[i] += N - p
            end
          else
            if @ok_count[p] != 0
              # @eval[i] -= (N - p) * ((@ng_count[p] + 1) ** 2)
              if p <= @best_result.score && @ng_count[p] == 2
                @eval[i] -= INF * 2
              end
            end
          end
        end
      end
      # debug(@eval)
      cands = 3.times.select { |i| @eval[i] > -INF }.to_a
      return if cands.empty?
      @dfs_count += 1 if cands.size > 1
      cands.sort_by { |i| -@eval[i] }.each do |ti|
        t = @terms[pos][ti]
        @ans[t.var] = t.val
        @var_pos[t.var].each do |ps|
          p, v = ps
          if p > pos
            if v == t.val
              @ok_count[p] += 1
            else
              @ng_count[p] += 1
            end
          end
        end
        dfs(pos + 1)
        @var_pos[t.var].each do |ps|
          p, v = ps
          if p > pos
            if v == t.val
              @ok_count[p] -= 1
            else
              @ng_count[p] -= 1
            end
          end
        end
        @ans[t.var] = -1
      end
    end
  end

  def shuffle(a)
    (a.size - 1).times do |i|
      pos = RND.rand(a.size - i) + i
      a[i], a[pos] = a[pos], a[i]
    end
  end

  def score(ans)
    N.times do |i|
      if @terms[i].all? { |t| ans[t.var] != t.val }
        return i
      end
    end
    return N
  end
end

def main
  solver = Solver.new
  res = solver.solve(START_TIME + TL)
  puts res
  {% if flag?(:local) %}
    STDERR.puts("final_score:#{res.score}")
  {% end %}
end

main

yukicoder

結果

テストケース

ソースコード