[DP] SRM 452 Div1 Hard IncreasingNumber
12182 단어 동적 기획
Solution S o l u t i o n
처음에 O(123)m2logn) O(123)m 2logn)를 스스로 뇌보정하는 방법을 썼다.대개 dpi,lo,hi,jdpi,lo,hi,j는 ii위를 고려하고 lolo로 시작하고hihi로 시작하며 모드 m는 jj의 방안수이다.그리고 DP를 배로 늘려요.하지만 지나갈 수 없어...
하나의 수 i, 1≤i<10i, 1≤i <10은 10ai - 19 10 a i - 119 이런 것만 선택하면 됩니다.이거 1111...111 1111....111모mm는 순환절이 있을 수 있다??나는 잉여계의 개수를 배로 늘렸다.그리고 가방만 있으면 돼요.주의해야 할 것은 가장 긴 그 1111.111 1111.111은 반드시 선택해야 한다.폭발 int+1
// BEGIN CUT HERE
// END CUT HERE
#line 5 "IncreasingNumber.cpp"
#include
using namespace std;
typedef long long ll;
const int N = 555;
const int MOD = 1000000007;
typedef pair<int, int> pairs;
ll n;
int m, ans, tmp;
int dp[N][N][10];
int f[N], f1[N];
int val[N][10];
int inv[N];
class IncreasingNumber {
public:
inline void pre(int n) {
inv[1] = 1;
for (int i = 2; i <= n; i++)
inv[i] = (ll)(MOD - MOD / i) * inv[MOD % i] % MOD;
}
inline int C(ll n, int m) {
int res = 1; n %= MOD;
for (int i = 0; i < m; i++) {
res = (n - i) * res % MOD;
res = (ll)res * inv[i + 1] % MOD;
}
return (res + MOD) % MOD;
}
inline void add(int &x, int a) {
x += a; while (x >= MOD) x -= MOD;
}
inline pairs solve(ll n) {
if (n == 0) return pairs(0, 0);
if (n == 1) {
f[1 % m] = 1;
return pairs(1 % m, 1 % m);
}
if (n & 1) {
pairs p = solve(n - 1);
int base = p.first * 10 % m, _base = (p.second * 10 + 1) % m;
add(f[_base], 1);
return pairs(base, _base);
} else {
pairs p = solve(n / 2);
int base = p.first * 10, _base = p.second;
for (int i = 0; i < m; i++) f1[i] = f[i];
for (int i = 0; i < m; i++)
add(f[(i * base + _base) % m], f1[i]);
return pairs(base * base / 10 % m, _base * (base + 1) % m);
}
}
int countNumbers(long long digits, int divisor) {
n = digits; m = divisor; pre(10);
tmp = (solve(n - 1).second * 10 + 1) % m;
for (int i = 0; i < m; i++)
for (int j = 0; j < 10; j++)
val[i][j] = C(f[i] + j - 1, j);
dp[0][0][0] = 1;
for (int i = 0; i < m; i++)
for (int j = 0; j < m; j++)
for (int k = 0; k < 10; k++)
for (int x = 0; x + k < 10; x++)
add(dp[i + 1][(j + i * x) % m][k + x], (ll)val[i][x] * dp[i][j][k] % MOD);
for (int i = 1; i < 10; i++)
for (int j = 0; j + i < 10; j++)
add(ans, dp[m][(m - tmp * i % m) % m][j]);
return ans;
}
// BEGIN CUT HERE
public:
void run_test(int Case) { if ((Case == -1) || (Case == 0)) test_case_0(); if ((Case == -1) || (Case == 1)) test_case_1(); if ((Case == -1) || (Case == 2)) test_case_2(); if ((Case == -1) || (Case == 3)) test_case_3(); }
private:
template <typename T> string print_array(const vector &V) { ostringstream os; os << "{ "; for (typename vector ::const_iterator iter = V.begin(); iter != V.end(); ++iter) os << '\"' << *iter << "\","; os << " }"; return os.str(); }
void verify_case(int Case, const int &Expected, const int &Received) { cerr << "Test Case #" << Case << "..."; if (Expected == Received) cerr << "PASSED" << endl; else { cerr << "FAILED" << endl; cerr << "\tExpected: \"" << Expected << '\"' << endl; cerr << "\tReceived: \"" << Received << '\"' << endl; } }
void test_case_0() { long long Arg0 = 2LL; int Arg1 = 12; int Arg2 = 4; verify_case(0, Arg2, countNumbers(Arg0, Arg1)); }
void test_case_1() { long long Arg0 = 3LL; int Arg1 = 111; int Arg2 = 9; verify_case(1, Arg2, countNumbers(Arg0, Arg1)); }
void test_case_2() { long long Arg0 = 452LL; int Arg1 = 10; int Arg2 = 0; verify_case(2, Arg2, countNumbers(Arg0, Arg1)); }
void test_case_3() { long long Arg0 = 6LL; int Arg1 = 58; int Arg2 = 38; verify_case(3, Arg2, countNumbers(Arg0, Arg1)); }
// END CUT HERE
};
// BEGIN CUT HERE
int main(void) {
IncreasingNumber ___test;
___test.run_test(-1);
system("pause");
}
// END CUT HERE
이 내용에 흥미가 있습니까?
현재 기사가 여러분의 문제를 해결하지 못하는 경우 AI 엔진은 머신러닝 분석(스마트 모델이 방금 만들어져 부정확한 경우가 있을 수 있음)을 통해 가장 유사한 기사를 추천합니다:
01 가방, 완전 가방, 다중 가방 dp(동적 기획 입문 dp)01 가방은 2진법으로 직접 표시할 수 있지만 데이터 양이 너무 많으면 시간을 초과하는 것이 폭력이다.01 가방의 사상은 바로 이 물품에 대해 내가 넣은 가치가 큰지 안 넣은 가치가 큰지 비교하여 방정식 f[i][v...
텍스트를 자유롭게 공유하거나 복사할 수 있습니다.하지만 이 문서의 URL은 참조 URL로 남겨 두십시오.
CC BY-SA 2.5, CC BY-SA 3.0 및 CC BY-SA 4.0에 따라 라이센스가 부여됩니다.