This documentation is automatically generated by online-judge-tools/verification-helper
#define PROBLEM "https://yukicoder.me/problems/no/2219"
#include "./../../template/template.hpp"
#include "./../../math/modint.hpp"
using namespace mmrz;
using mint = modint<998244353>;
void mmrz::solve(){
string s;
cin >> s;
int n = len(s);
vector<vector<mint>> dp(n + 1, vector<mint>(4));
dp[0][0] = 1;
rep(i, n){
if(s[i] == '0' || s[i] == '?'){
dp[i + 1][1] += dp[i][0];
dp[i + 1][3] += dp[i][2];
}
if(s[i] == '1' || s[i] == '?'){
dp[i + 1][2] += dp[i][1];
}
rep(j, 4){
dp[i + 1][j] += dp[i][j] * (s[i] == '?' ? 2 : 1);
}
}
cout << dp[n][3] << '\n';
}#line 1 "verify/yukicoder/2219.test.cpp"
#define PROBLEM "https://yukicoder.me/problems/no/2219"
#line 1 "template/template.hpp"
# include <bits/stdc++.h>
using namespace std;
using ll = long long;
using ull = unsigned long long;
const double pi = acos(-1);
template<class T>constexpr T inf() { return ::std::numeric_limits<T>::max(); }
template<class T>constexpr T hinf() { return inf<T>() / 2; }
template <typename T_char>T_char TL(T_char cX) { return tolower(cX); }
template <typename T_char>T_char TU(T_char cX) { return toupper(cX); }
template<class T> bool chmin(T& a,T b) { if(a > b){a = b; return true;} return false; }
template<class T> bool chmax(T& a,T b) { if(a < b){a = b; return true;} return false; }
int popcnt(unsigned long long n) { int cnt = 0; for (int i = 0; i < 64; i++)if ((n >> i) & 1)cnt++; return cnt; }
int d_sum(ll n) { int ret = 0; while (n > 0) { ret += n % 10; n /= 10; }return ret; }
int d_cnt(ll n) { int ret = 0; while (n > 0) { ret++; n /= 10; }return ret; }
ll gcd(ll a, ll b) { if (b == 0)return a; return gcd(b, a%b); };
ll lcm(ll a, ll b) { ll g = gcd(a, b); return a / g*b; };
ll MOD(ll x, ll m){return (x%m+m)%m; }
ll FLOOR(ll x, ll m) {ll r = (x%m+m)%m; return (x-r)/m; }
template<class T> using dijk = priority_queue<T, vector<T>, greater<T>>;
# define all(qpqpq) (qpqpq).begin(),(qpqpq).end()
# define UNIQUE(wpwpw) (wpwpw).erase(unique(all((wpwpw))),(wpwpw).end())
# define LOWER(epepe) transform(all((epepe)),(epepe).begin(),TL<char>)
# define UPPER(rprpr) transform(all((rprpr)),(rprpr).begin(),TU<char>)
# define rep(i,upupu) for(ll i = 0, i##_len = (upupu);(i) < (i##_len);(i)++)
# define reps(i,opopo) for(ll i = 1, i##_len = (opopo);(i) <= (i##_len);(i)++)
# define len(x) ((ll)(x).size())
# define bit(n) (1LL << (n))
# define pb push_back
# define eb emplace_back
# define exists(c, e) ((c).find(e) != (c).end())
struct INIT{
INIT(){
std::ios::sync_with_stdio(false);
std::cin.tie(0);
cout << fixed << setprecision(20);
}
}INIT;
namespace mmrz {
void solve();
}
int main(){
mmrz::solve();
}
#line 2 "math/modint.hpp"
#line 5 "math/modint.hpp"
template <std::uint_fast64_t Modulus> class modint {
using u64 = std::uint_fast64_t;
public:
u64 a;
constexpr modint(const u64 x = 0) noexcept : a(x % Modulus) {}
constexpr u64 &value() noexcept { return a; }
constexpr const u64 &value() const noexcept { return a; }
constexpr modint operator+(const modint rhs) const noexcept {
return modint(*this) += rhs;
}
constexpr modint operator-(const modint rhs) const noexcept {
return modint(*this) -= rhs;
}
constexpr modint operator*(const modint rhs) const noexcept {
return modint(*this) *= rhs;
}
constexpr modint operator/(const modint rhs) const noexcept {
return modint(*this) /= rhs;
}
constexpr modint &operator+=(const modint rhs) noexcept {
a += rhs.a;
if (a >= Modulus) {
a -= Modulus;
}
return *this;
}
constexpr modint &operator-=(const modint rhs) noexcept {
if (a < rhs.a) {
a += Modulus;
}
a -= rhs.a;
return *this;
}
constexpr modint &operator*=(const modint rhs) noexcept {
a = a * rhs.a % Modulus;
return *this;
}
constexpr modint &operator/=(modint rhs) noexcept {
u64 exp = Modulus - 2;
while (exp) {
if (exp % 2) {
*this *= rhs;
}
rhs *= rhs;
exp /= 2;
}
return *this;
}
constexpr modint& operator++() noexcept {
if (++a == Modulus) a = 0;
return *this;
}
constexpr modint operator++(int) noexcept {
modint tmp(*this);
++(*this);
return tmp;
}
constexpr modint& operator--() noexcept {
if (a == 0) a = Modulus;
--a;
return *this;
}
constexpr modint operator--(int) noexcept {
modint tmp(*this);
--(*this);
return tmp;
}
friend std::ostream& operator<<(std::ostream& os, const modint& rhs) {
os << rhs.a;
return os;
}
};
#line 5 "verify/yukicoder/2219.test.cpp"
using namespace mmrz;
using mint = modint<998244353>;
void mmrz::solve(){
string s;
cin >> s;
int n = len(s);
vector<vector<mint>> dp(n + 1, vector<mint>(4));
dp[0][0] = 1;
rep(i, n){
if(s[i] == '0' || s[i] == '?'){
dp[i + 1][1] += dp[i][0];
dp[i + 1][3] += dp[i][2];
}
if(s[i] == '1' || s[i] == '?'){
dp[i + 1][2] += dp[i][1];
}
rep(j, 4){
dp[i + 1][j] += dp[i][j] * (s[i] == '?' ? 2 : 1);
}
}
cout << dp[n][3] << '\n';
}