#include <iostream>
#include <cstdio>
#include <cmath>
#include <ctime>
#include <cstdlib>
#include <cassert>
#include <vector>
#include <list>
#include <stack>
#include <queue>
#include <deque>
#include <map>
#include <set>
#include <bitset>
#include <string>
#include <algorithm>
#include <utility>
#define llint long long
#define inf 1e18
#define rep(x, s, t) for(llint (x) = (s); (x) < (t); (x)++)
#define Rep(x, s, t) for(llint (x) = (s); (x) <= (t); (x)++)
#define chmin(x, y) (x) = min((x), (y))
#define chmax(x, y) (x) = max((x), (y))

using namespace std;
typedef pair<llint, llint> P;

// range-add, range-min query
 
struct LazySegTree{
	int size;
	vector<llint> seg, delay;
	
	LazySegTree(){}
	LazySegTree(int size){
		this->size = size;
		seg.resize(1<<(size+1));
		delay.resize(1<<(size+1));
	}
	
	void init()
	{
		for(int i = 0; i < (1<<(size+1)); i++){
			seg[i] = 0;
			delay[i] = 0;
		}
	}
	
	void eval(int l, int r, int k)
	{
		if(delay[k] % 2){
			seg[k] = (r-l+1) - seg[k];  //総和クエリのときは区間幅を乗じる必要がある
			if(l < r){
				delay[k*2] += delay[k];
				delay[k*2+1] += delay[k];
			}
			delay[k] = 0;
		}
	}
	
	void update(int i, llint val)
	{
		int l = 0, r = (1<<size)-1, k = 1;
		eval(l, r, k);
		for(int j = size-1; j >= 0; j--){
			k <<= 1;
			if(i & (1<<j)){
				k++;
				r = (l+r)/2;
			}
			else l = (l+r)/2+1;
			eval(l, r, k);
		}
		
		i += (1 << size);
		seg[i] = val;
		while(i > 1){
			i /= 2;
			seg[i] = (seg[i*2] + seg[i*2+1]);
		}
	}
	
	void add(int a, int b, int k, int l, int r, llint val)
	{
		eval(l, r, k);
		
		if(b < l || r < a) return;
		if(a <= l && r <= b){
			delay[k] += val;
			eval(l, r, k);
			return;
		}
		add(a, b, k*2, l, (l+r)/2, val);
		add(a, b, k*2+1, (l+r)/2+1, r, val);
		seg[k] = (seg[k*2] + seg[k*2+1]);
	}
	void add(int a, int b, llint val){
		if(a > b) return;
		add(a, b, 1, 0, (1<<size)-1, val);
	}
 
	llint query(int a, int b, int k, int l, int r)
	{
		eval(l, r, k);
		
		if(b < l || r < a) return 0;
		if(a <= l && r <= b) return seg[k];
		llint lval = query(a, b, k*2, l, (l+r)/2);
		llint rval = query(a, b, k*2+1, (l+r)/2+1, r);
		return (lval + rval);
	}
	llint query(int a, int b)
	{
		return query(a, b, 1, 0, (1<<size)-1);
	}
};

string s;
LazySegTree seg(18);
llint zcnt[200005];

int main(void)
{
	ios::sync_with_stdio(0);
	cin.tie(0);
	
	cin >> s;
	llint n = s.size();
	s = "#" + s;
	
	llint add = 0, one = 0; string t;
	for(int i = 1; i <= n; i++){
		if(s[i] == '0') t += s[i];
		else{
			one++;
			if(i < n && s[i+1] == '1'){
				add++;
				i++;
				one++;
			}
			else t += s[i];
		}
	}
	
	seg.init();
	llint m = 0, cnt = 0;
	for(int i = 0; i < t.size(); i++){
		if(t[i] == '0') zcnt[m]++;
		else{
			m++;
			seg.update(m, (i+1)%2);
		}
	}
	
	llint ans = 0;
	for(int i = 1; i < one; i++) ans += (i+1)/2;
	ans += seg.query(1, m) + add;
	
	llint sum = 0;
	for(int i = 1; i <= m; i++){
		sum += zcnt[i-1];
		if(seg.query(i, i) == 0 && sum > 0){
			seg.add(i, m, 1), sum--;
			ans += seg.query(1, m) + add;
			zcnt[i-1]--;
		}
	}
	
	for(int i = m; i >= 0; i--){
		for(int j = 0; j < zcnt[i]; j++){
			seg.add(i+1, m, 1);
			ans += seg.query(1, m) + add;
		}
	}
	cout << ans << endl;
	
	return 0;
}

Measured memory_bytes: 12940