# input if n<m: n,m=m,n #n>m # k cannot be 1, always could be n+m res=1 for k in range(2,n+m): q,rem=divmod(n+m,k) hb=min(n//q,int(k-m/(q+1))) lb=max((math.ceil(n/(q+1)),math.ceil(k-m/q),(k+1)//2)) if hb>=lb: res+=1 print(res)