#include <bits/stdc++.h>
#define FORD(i, a, b) for (ll i = a; i >= (ll)(b); i--)
using namespace std;
typedef long long ll;

ll c_to_i(char c)
{
  return (ll)c - (ll)'0';
}

class mint
{
  const ll mod = 2019;

public:
  ll x;

  mint(ll x = 0) : x((x % mod + mod) % mod) {}
  mint operator-() const
  {
    return mint(-x);
  }
  mint &operator+=(const mint &a)
  {
    if ((x += a.x) >= mod)
      x -= mod;
    return *this;
  }
  mint &operator-=(const mint &a)
  {
    if ((x += mod - a.x) >= mod)
      x -= mod;
    return *this;
  }
  mint &operator*=(const mint &a)
  {
    (x *= a.x) %= mod;
    return *this;
  }
  mint operator+(const mint &a) const
  {
    mint res(*this);
    return res += a;
  }
  mint operator-(const mint &a) const
  {
    mint res(*this);
    return res -= a;
  }
  mint operator*(const mint &a) const
  {
    mint res(*this);
    return res *= a;
  }
  mint pow(ll t) const
  {
    if (!t)
      return 1;
    mint a = pow(t >> 1);
    a *= a;
    if (t & 1)
      a *= *this;
    return a;
  }
  // for prime mod
  mint inv() const
  {
    return pow(mod - 2);
  }
  mint &operator/=(const mint &a)
  {
    return (*this) *= a.inv();
  }
  mint operator/(const mint &a) const
  {
    mint res(*this);
    return res /= a;
  }
  friend ostream &operator<<(ostream &os, const mint &m)
  {
    os << m.x;
    return os;
  }
};

int main()
{
  string s;
  cin >> s;

  ll n = s.length();

  mint *dp = new mint[n + 1];
  mint base = 1;
  map<ll, ll> ms;

  dp[n] *= 0;
  ms[0] += 1;
  FORD(i, n - 1, 0)
  {
    dp[i] *= 0;
    dp[i] += dp[i + 1] + base * c_to_i(s[i]);
    base *= 10;
    ms[dp[i].x] += 1;
  }

  ll ans = 0;
  for (auto m : ms)
  {
    ans += m.second * (m.second - 1) / 2;
  }

  cout << ans << endl;

  return 0;
}