結果

問題 No.1899 L1 Cafe
ユーザー 👑 KazunKazun
提出日時 2022-04-08 22:45:19
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 224 ms / 2,000 ms
コード長 1,114 bytes
コンパイル時間 218 ms
コンパイル使用メモリ 82,452 KB
実行使用メモリ 82,052 KB
最終ジャッジ日時 2024-11-28 13:35:32
合計ジャッジ時間 4,090 ms
ジャッジサーバーID
(参考情報)
judge1 / judge3
このコードへのチャレンジ
(要ログイン)
ファイルパターン 結果
sample AC * 1
other AC * 12
権限があれば一括ダウンロードができます

ソースコード

diff #

def Linear_Sum(a,b,L,R):
    """ Sum_{i=L}^R (ai+b) を求める.
    """

    if L<=R:
        return (a*(L+R)+2*b)*(R-L+1)//2
    else:
        return 0

def floor_sum(A,B,M,N):
    """sum_{i=0}^{N-1} floor((A*i+B)/M) を求める.
    """
    T=0
    while True:
        T+=((N-1)*N//2)*(A//M)
        A%=M

        T+=N*(B//M)
        B%=M

        y=(A*N+B)//M
        x=B-y*M

        if y==0:
            return T

        T+=(N+x//A)*y
        A,B,M,N=M,x%A,A,y

def Floor_Sum(A:int,B:int,M:int,N:int,K=0):
    """sum_{i=K}^N floor((A*i+B)/M) を求める.
    """

    if K==0:
        return floor_sum(A,B,M,N+1)
    else:
        return floor_sum(A,B,M,N+1)-floor_sum(A,B,M,K)
#==================================================
import sys
input=sys.stdin.readline
write=sys.stdout.write

T=int(input())
Ans=[0]*T
for t in range(T):
    N,A,B=map(int,input().split())

    # X 軸に平行な道路上
    Ans[t]+=Linear_Sum(-B,N,1,N//B)

    # Y 軸に平行な道路上
    Ans[t]+=Linear_Sum(-A,N,1,N//A)

    # 道路の交差点
    Ans[t]-=Floor_Sum(-A,N,B,N//A,1)

write("\n".join(map(str,Ans)))
0