Pagini recente » Cod sursa (job #2177030) | Cod sursa (job #3302970) | Cod sursa (job #3321953) | Cod sursa (job #1669535) | Cod sursa (job #3346659)
// hai sa nu trollam labubu de pa
#include <bits/stdc++.h>
using namespace std;
const long long max_size = 5e3 + 20, mod = 666013;
int dp[max_size][max_size], ways[max_size][max_size];
void solve()
{
string a, b;
cin >> a >> b;
int n = a.size(), m = b.size();
a = " " + a;
b = " " + b;
for (int i = 0; i <= n; i++)
{
dp[i][0] = 0;
ways[i][0] = 1;
}
for (int j = 0; j <= m; j++)
{
dp[0][j] = 0;
ways[0][j] = 1;
}
for (int i = 1; i <= n; i++)
{
for (int j = 1; j <= m; j++)
{
if (a[i] == b[j])
{
dp[i][j] = dp[i - 1][j - 1] + 1;
ways[i][j] = ways[i - 1][j - 1];
}
else
{
if (dp[i - 1][j] > dp[i][j - 1])
{
dp[i][j] = dp[i - 1][j];
ways[i][j] = ways[i - 1][j];
}
else
{
if (dp[i - 1][j] < dp[i][j - 1])
{
dp[i][j] = dp[i][j - 1];
ways[i][j] = ways[i][j - 1];
}
else
{
dp[i][j] = dp[i - 1][j];
ways[i][j] = (ways[i - 1][j] + ways[i][j - 1]) % mod;
if (dp[i - 1][j - 1] == dp[i][j]) // evitam overcount
{
ways[i][j] = (ways[i][j] - ways[i - 1][j - 1] + mod) % mod;
}
}
}
}
}
}
cout << ways[n][m];
}
signed main()
{
#ifdef LOCAL
freopen("test.in", "r", stdin);
freopen("test.out", "w", stdout);
#else
freopen("subsir.in", "r", stdin);
freopen("subsir.out", "w", stdout);
#endif // LOCAL
ios_base::sync_with_stdio(false);
cin.tie(0);
cout.tie(0);
long long tt;
//cin >> tt;
tt = 1;
while (tt--)
{
solve();
}
return 0;
}