#!/usr/bin/env python3
# from typing import *

import sys
import io
import math
import collections
import decimal
import itertools
import bisect
import heapq


def input():
    return sys.stdin.readline()[:-1]


# sys.setrecursionlimit(1000000)

# _INPUT = """# paste here...
# """
# sys.stdin = io.StringIO(_INPUT)

INF = 10**10

MOD = 10**9+7

A, B = map(int, input().split())
if A == 0:
    ans = B

elif B == 0:
    ans = A

else:
    A %= MOD
    B %= MOD
    ans = ((A+1)*(B+1)+A*B) % MOD

print(ans)