結果

問題 No.546 オンリー・ワン
ユーザー タコイモタコイモ
提出日時 2023-03-01 21:58:21
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 127 ms / 2,000 ms
コード長 1,698 bytes
コンパイル時間 1,384 ms
コンパイル使用メモリ 86,648 KB
実行使用メモリ 77,728 KB
最終ジャッジ日時 2023-10-15 01:33:14
合計ジャッジ時間 2,581 ms
ジャッジサーバーID
（参考情報） 		judge12 / judge15
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 103 ms
71,636 KB
testcase_01 AC 103 ms
71,256 KB
testcase_02 AC 103 ms
71,704 KB
testcase_03 AC 99 ms
71,756 KB
testcase_04 AC 101 ms
71,460 KB
testcase_05 AC 99 ms
71,320 KB
testcase_06 AC 102 ms
71,880 KB
testcase_07 AC 116 ms
77,636 KB
testcase_08 AC 118 ms
77,488 KB
testcase_09 AC 127 ms
77,728 KB
testcase_10 AC 123 ms
77,604 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

def popcount(x):
    '''xの立っているビット数をカウントする関数
    (xは64bit整数)'''

    # 2bitごとの組に分け、立っているビット数を2bitで表現する
    x = x - ((x >> 1) & 0x5555555555555555)

    # 4bit整数に 上位2bit + 下位2bit を計算した値を入れる
    x = (x & 0x3333333333333333) + ((x >> 2) & 0x3333333333333333)

    x = (x + (x >> 4)) & 0x0f0f0f0f0f0f0f0f # 8bitごと
    x = x + (x >> 8) # 16bitごと
    x = x + (x >> 16) # 32bitごと
    x = x + (x >> 32) # 64bitごと = 全部の合計
    return x & 0x0000007f
import sys
#sys.setrecursionlimit(500000)
def I(): return int(sys.stdin.readline().rstrip())
def MI(): return map(int,sys.stdin.readline().rstrip().split())
def TI(): return tuple(map(int,sys.stdin.readline().rstrip().split()))
def LI(): return list(map(int,sys.stdin.readline().rstrip().split()))
def S(): return sys.stdin.readline().rstrip()
def LS(): return list(sys.stdin.readline().rstrip())
#for i, pi in enumerate(p):
from collections import defaultdict,deque
import bisect
import itertools
dic = defaultdict(int)
d = deque()
YN = ['No','Yes']

N,L,R = MI()
C = LI()
n2 = 2**N
dp = [0]*(n2)#少なくともbitが立っている部分では割り切れる奴の個数
import math
for i in range(n2):
  s = 1
  for j in range(N):
    if (1<< j)& i:
      cj = C[j]
      s = (s*cj)//math.gcd(s,cj)
  
  dp[i] = (R//s - (L+s-1)//s)+1
#すべての和集合*N-iを除いた和集合
S = [0]*(N)
ans = 0
#print(dp)
for i in range(1,n2):
  x = popcount(i)
  
  ans += dp[i]*((-1)**(x+1))
  for j in range(N):
    if not (1<< j)& i:
      S[j] += -dp[i]*((-1)**x)
print(ans*N-sum(S))   
#print(S,ans)
0