ICode9

精准搜索请尝试: 精确搜索
首页 > 其他分享> 文章详细

2013 ACM-ICPC吉林通化全国邀请赛 J Dice 概率DP + 数列

2021-03-29 16:03:53  阅读:215  来源: 互联网

标签:Dice dp1 int ll ICPC DP ans dp MOD


HDU-4599

在这里插入图片描述

solution

已 知   d p i [ n ] = 0 , ( i = 1 , 2 ) 已知~dp_i[n]=0,(i=1,2) 已知 dpi​[n]=0,(i=1,2)

d p 1 [ i ] = d p 1 [ i + 1 ] ∗ 5 6 + d p 1 [ 1 ] ∗ 1 6 + 1 dp_1[i]=dp_1[i+1]*\frac56+dp_1[1]*\frac16+1 dp1​[i]=dp1​[i+1]∗65​+dp1​[1]∗61​+1

= > d p 1 [ n ] − ( d p 1 [ 1 ] + 5 6 ) = 6 [ d p 1 [ n − 1 ] − ( d p 1 [ 1 ] + 5 6 ) ] , d p 1 [ 0 ] = d p 1 [ 1 ] + 1 =>dp_1[n]-(dp_1[1]+\frac56)=6[dp_1[n-1]-(dp_1[1]+\frac56)],dp_1[0]=dp_1[1]+1 =>dp1​[n]−(dp1​[1]+65​)=6[dp1​[n−1]−(dp1​[1]+65​)],dp1​[0]=dp1​[1]+1

= > d p 1 [ n ] − d p 1 [ 1 ] − 6 5 = − 6 n 5 =>dp_1[n]-dp_1[1]-\frac65=-\frac{6^n}{5} =>dp1​[n]−dp1​[1]−56​=−56n​

= > d p 1 [ 1 ] + 6 5 = 6 n 5 = > d p 1 [ 0 ] = 6 n − 1 5 =>dp_1[1]+\frac65=\frac{6^n}{5}=>dp_1[0]=\frac{6^n-1}5 =>dp1​[1]+56​=56n​=>dp1​[0]=56n−1​

同理,

d p 2 [ 0 ] = 6 ∗ 6 n − 1 5 dp_2[0]=6*\frac{6^n-1}5 dp2​[0]=6∗56n−1​

显 然 G ( n ) = 6 ∗ n 显然G(n)=6*n 显然G(n)=6∗n

求 最 小 的 M 1 , M 2 使 得 , G ( M 1 ) > = F ( N ) , G ( M 2 ) > = H ( N ) , 求最小的M_1,M_2使得,G (M_1) >= F (N),G(M_2)>=H(N), 求最小的M1​,M2​使得,G(M1​)>=F(N),G(M2​)>=H(N),

其 中 F ( N ) = d p 1 [ 0 ] , H ( N ) = d p 2 [ 0 ] . 其中F(N)=dp_1[0],H(N)=dp_2[0]. 其中F(N)=dp1​[0],H(N)=dp2​[0].

于是有,

M 1 > = c e i l ( 6 n − 1 30 ) = 6 n + 24 30 , M 2 > = 6 n − 1 5 M_1>=ceil(\frac{6^n-1}{30})=\frac{6^n+24}{30},M_2>=\frac{6^n-1}5 M1​>=ceil(306n−1​)=306n+24​,M2​>=56n−1​

code

/*Siberian Squirrel*/
/*Cute JinFish*/
#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const double PI = acos(-1), eps = 1e-6;
/*const int MOD = 998244353, r = 119, k = 23, g = 3;
const int MOD = 1004535809, r = 479, k = 21, g = 3;*/
const int INF = 0x3f3f3f3f, MOD = 2011;
const int M = 3e2 + 10, N = 5e6 + 10;
int sgn(double x) {
    if(fabs(x) < eps) return 0;
    return x < 0? -1: 1;
}
//inline int rnd(){static int seed=2333;return seed=(((seed*666666ll+20050818)%998244353)^1000000007)%1004535809;}
//double Rand() {return (double)rand() / RAND_MAX;}

ll n;

ll quick_pow(ll ans, ll p, ll res = 1) {
    ans %= MOD, p %= MOD - 1;
    for(; p; p >>= 1, ans = ans * ans % MOD)
        if(p & 1) res = res * ans % MOD;
    return res % MOD;
}

ll inv(ll ans) {
    return quick_pow(ans, MOD - 2);
}

void init() {}

void solve(double _p = 0) {
    cout << ((quick_pow(6, n) + 24 + MOD) % MOD * inv(30) % MOD) << ' ';
    cout << (quick_pow(6, n) - 1 + MOD) % MOD * inv(5) % MOD << endl;
}
/*


*/
int main() {
//    ios::sync_with_stdio(false);cin.tie(0);cout.tie(nullptr);
// srand(time(0));
#ifdef ACM_LOCAL
    freopen("input", "r", stdin);
    freopen("output", "w", stdout);
#endif
    init();
    int o = 1;
//	cin >> o;
    while(o --) {
        while(cin >> n && n) {
            solve();
        }
    }
    return 0;
}

标签:Dice,dp1,int,ll,ICPC,DP,ans,dp,MOD
来源: https://blog.csdn.net/qq_46173805/article/details/115301971

本站声明: 1. iCode9 技术分享网(下文简称本站)提供的所有内容,仅供技术学习、探讨和分享;
2. 关于本站的所有留言、评论、转载及引用,纯属内容发起人的个人观点,与本站观点和立场无关;
3. 关于本站的所有言论和文字,纯属内容发起人的个人观点,与本站观点和立场无关;
4. 本站文章均是网友提供,不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属;如您发现该文章侵犯了您的权益,可联系我们第一时间进行删除;
5. 本站为非盈利性的个人网站,所有内容不会用来进行牟利,也不会利用任何形式的广告来间接获益,纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。

专注分享技术,共同学习,共同进步。侵权联系[81616952@qq.com]

Copyright (C)ICode9.com, All Rights Reserved.

ICode9版权所有