#include <bits/stdc++.h>
using namespace std;
const int MOD = 1e5;
const int INF = 0x3f3f3f3f;
string to_upper(string a) {
  for (int i = 0; i < (int)a.size(); ++i)
    if (a[i] >= 'a' && a[i] <= 'z') a[i] -= 'a' - 'A';
  return a;
}
string to_lower(string a) {
  for (int i = 0; i < (int)a.size(); ++i)
    if (a[i] >= 'A' && a[i] <= 'Z') a[i] += 'a' - 'A';
  return a;
}
long long sti(string a, int base = 10) {
  long long k = 0;
  for (int i = 0; i < (int)a.size(); ++i) {
    k *= base;
    k += a[i] - '0';
  }
  return k;
}
long long min(long long a, int b) {
  if (a < b) return a;
  return b;
}
long long min(int a, long long b) {
  if (a < b) return a;
  return b;
}
long long max(long long a, int b) {
  if (a > b) return a;
  return b;
}
long long max(int a, long long b) {
  if (a > b) return a;
  return b;
}
long long gcd(long long a, long long b) {
  if (b == 0) return a;
  return gcd(b, a % b);
}
long long lcm(long long a, long long b) { return a / gcd(a, b) * b; }
inline int multt(long long a, long long b) { return (a * b) % MOD; }
inline int pott(int n, int k) {
  if (k == 0) return 1;
  long long a = pott(n, k / 2);
  a = multt(a, a);
  if (k % 2)
    return multt(a, n);
  else
    return a;
}
int divide(int a, int b) { return multt(a, pott(b, MOD - 2)); }
inline int subb(int a, int b) {
  if (a - b >= 0) return a - b;
  return a - b + MOD;
}
inline int addd(int a, int b) {
  if (a + b >= MOD) return a + b - MOD;
  return a + b;
}
bool prime(long long a) {
  if (a == 1) return 0;
  for (int i = 2; i <= round(sqrt(a)); ++i) {
    if (a % i == 0) return 0;
  }
  return 1;
}
const int maxn = 4e5 + 5;
long long a[maxn];
int main() {
  int t;
  cin >> t;
  while (t--) {
    long long n, m;
    cin >> n >> m;
    long long a = n - m;
    long long res = 0;
    if (m == 0) {
      cout << 0 << endl;
      continue;
    }
    if (m + 1 >= a) {
      res += (n) * (n + 1) / 2;
      res -= a;
      cout << res << endl;
      continue;
    }
    long long b = a / (m + 1);
    long long c = a % (m + 1);
    long long d = (m + 1 - c);
    res = (n + 1) * n / 2 - d * (b + 1) * b / 2 - c * (b + 1) * (b + 2) / 2;
    cout << res << endl;
  }
  return 0;
}
