結果
問題 | No.1252 数字根D |
ユーザー | 👑 Kazun |
提出日時 | 2020-07-15 04:07:51 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 179 ms / 2,000 ms |
コード長 | 2,891 bytes |
コンパイル時間 | 255 ms |
コンパイル使用メモリ | 86,812 KB |
実行使用メモリ | 77,968 KB |
最終ジャッジ日時 | 2023-08-07 10:29:07 |
合計ジャッジ時間 | 3,382 ms |
ジャッジサーバーID (参考情報) |
judge14 / judge11 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 70 ms
71,524 KB |
testcase_01 | AC | 76 ms
75,384 KB |
testcase_02 | AC | 83 ms
75,780 KB |
testcase_03 | AC | 144 ms
77,768 KB |
testcase_04 | AC | 177 ms
77,920 KB |
testcase_05 | AC | 179 ms
77,876 KB |
testcase_06 | AC | 175 ms
77,860 KB |
testcase_07 | AC | 176 ms
77,820 KB |
testcase_08 | AC | 174 ms
77,968 KB |
testcase_09 | AC | 176 ms
77,820 KB |
testcase_10 | AC | 179 ms
77,812 KB |
testcase_11 | AC | 178 ms
77,864 KB |
testcase_12 | AC | 179 ms
77,912 KB |
testcase_13 | AC | 175 ms
77,820 KB |
testcase_14 | AC | 71 ms
71,448 KB |
ソースコード
class Modulo_Error(Exception): pass class Modulo(): def __init__(self,a,n): self.a=a%n self.n=n def __str__(self): return "{} (mod {})".format(self.a,self.n) #+,- def __pos__(self): return self def __neg__(self): return Modulo(-self.a,self.n) #等号,不等号 def __eq__(self,other): if isinstance(other,Modulo): return (self.a==other.a) and (self.n==other.n) elif isinstance(other,int): return (self-other).a==0 def __neq__(self,other): return not(self==other) #加法 def __add__(self,other): if isinstance(other,Modulo): if self.n!=other.n: raise Modulo_Error("異なる法同士の演算です.") return Modulo(self.a+other.a,self.n) elif isinstance(other,int): return Modulo(self.a+other,self.n) def __radd__(self,other): if isinstance(other,int): return Modulo(self.a+other,self.n) #減法 def __sub__(self,other): return self+(-other) def __rsub__(self,other): if isinstance(other,int): return -self+other #乗法 def __mul__(self,other): if isinstance(other,Modulo): if self.n!=other.n: raise Modulo_Error("異なる法同士の演算です.") return Modulo(self.a*other.a,self.n) elif isinstance(other,int): return Modulo(self.a*other,self.n) def __rmul__(self,other): if isinstance(other,int): return Modulo(self.a*other,self.n) #Modulo逆数 def Modulo_Inverse(self): x0, y0, x1, y1 = 1, 0, 0, 1 a,b=self.a,self.n while b != 0: q, a, b = a // b, b, a % b x0, x1 = x1, x0 - q * x1 y0, y1 = y1, y0 - q * y1 if a!=1: raise Modulo_Error("{}の逆数が存在しません".format(self)) else: return Modulo(x0,self.n) #除法 def __truediv__(self,other): return self*other.Modulo_Inverse() #累乗 def __pow__(self,m): u=abs(m) r=Modulo(1,self.n) while u>0: if u%2==1: r*=self self*=self u=u>>1 if m>=0: return r else: return r.Modulo_Inverse() #--------------------------------------------------------------------- M=2**127-1 T=int(input()) def g(N): if N==0: return Modulo(0,M) if N%(D-1)==0: Q=Modulo(N//(D-1),M) R=Modulo(0,M) else: Q=Modulo(N//(D-1)+1,M) R=Modulo(D-1-N%(D-1),M) t=Modulo(2,M) return Q*(E*(E-1))/t-E*R+(R*(R+1))/t for _ in range(T): D,A,B=map(int,input().split()) E=Modulo(D,M) print((g(B)-g(max(A-1,0))).a)