国产精品视频专区,欧美高清视频一区二区三区在线观看 ,久久精品国产免费观看

CF54C

_Yuan — Wed, 04 May 2011 17:06:00 GMT

/*
    区间[L,R]中每个数出现的概率�ؓ1/(R-L+1)
    现给出n个区��_��[Li,Ri]�Q�每个区间抽取出一个数
    求这n个数中，臛_��有K%个数的前�~��?的概�?br>
    可以处理出每个区间[Li, Ri]抽取的数前缀�?的概率p[i]
    dp[i,j]表示前面i个区间有j个数前缀�?的概�?br>    则�{�U�d��是dp[i,j] = dp[i-1,j-1]*p[i]  + dp[i-1,j]*(1-p[i])
    �{�案��是∑dp[n,j]�Q?nbsp;j/n >= K%
*/
#include<iostream>
#include<cstring>
#include<map>
#include<algorithm>
#include<stack>
#include<queue>
#include<cmath>
#include<string>
#include<cstdlib>
#include<vector>
#include<cstdio>
#include<set>
#include<list>
#include<numeric>
#include<cassert>
#include<ctime>
#include<bitset>

using namespace std;

long long cal(long long x) {//[1,x]
    long long cnt = 0 , now = 1;
    while(x/now > 1) {
        cnt += now;
        now *= 10;
    }
    if(x / now == 1) {//x = 1***
        cnt += x - now + 1;
    }
    return cnt;
}

double cal(long long left, long long right) {
    return (cal(right) - cal(left-1)) / (right - left + 1.0);
}

double p[1000], dp[1000][1000];

int main()
{
#ifndef ONLINE_JUDGE
    freopen("in","r",stdin);
#endif

    for (int n; ~scanf("%d", &n); ) {
        long long left, right;
        for (int i = 0; i <n ; i++) {
            scanf("%I64d%I64d", &left, &right);
            p[i] = cal(left, right);
        }
        memset(dp, 0, sizeof dp);
        dp[0][0] = 1.0;
        for (int i = 1; i <= n; i++) {
            for(int j = 0; j <= i; j++) {
                dp[i][j] = dp[i-1][j-1]*p[i-1] + dp[i-1][j]*(1-p[i-1]);
            }
        }
        int per;
        scanf("%d", &per);
        double ans = 0;
        for( int j = 0 ; j <= n; j++ ) {
            if(100*j >= n*per) {
                ans += dp[n][j];
            }
        }
        printf("%.10f\n", ans);
    }
    return 0;
}

_Yuan 2011-05-05 01:06 发表评论

八中 1833 zjoi 2010 count

_Yuan — Wed, 23 Mar 2011 05:32:00 GMT

/*
    求[a,b]�?-9各数字出现的�ơ数    <=10^12

    一开始我用数位统计搞
    ��x��照这�U�方法搞�?br>    http://www.shnenglu.com/Yuan/archive/2011/01/25/139299.html
    搞了好久搞不出，逐渐发现dfs型的写法�Q�状态一般需�?pos,pre)�Q�需要压�~�前�~�为pre
    我定�?pos, preCnt)
    preCnt是表�C�前�~�中有多少个数字跟当前要检查的数字相同�Q�不行的�Q�不能正��表�C�出前缀!
    �?2,21记录1出现的个数都是一��L��Q�但是结果不�?br>    会被if(dp[pos][pre] != -1)
        return dp[pos][pre]
    �q�回错误的结�?br>
    �q�个前缀应该��h��q��性，也就说不同前�~�压羃得到的结果有些要一��P��q�样才能记忆化嘛�Q�重复利用！�Q?br>    但这题，前缀如果是直接用前面的数字来表示的话�Q�很多，肯定时空都不�?br>

    所以，搞数位统计，需要记录前�~��Q�数位统计不记录应该前缀不行吧）
    同时前缀必须是能区分�Q�不会返回错误结果）的，�q�要可共用（用来减少�q�算�Q�提高速度�Q?br>
    �q�道题的正确做法�Q�换做以前，可能会想刎ͼ�但是现在感觉思维定势了T_T
    参考http://hi.baidu.com/736551725/blog/item/8ba4408e4fb74206b21bbae4.html
    观察下表
    000
    001
    .
    998
    999
    对于长度�?的，��d��?*10^3个数字，每个数字出现的概率一��P��所以每个数字出现的�ơ数�?*10^3/10
    也即宽度为n的表�Q�每个数字出现的�ơ数为n*10^(n-1)
    对于前缀0�Q�要��L��Q�对于宽度�ؓn的表�Q�前�~�0的个��Cؓ11..1个（表的左上角那�?�Q�包�?00那个�Q?br>
    有了上面的观察，��逐位�l�计下，��p��得到[1,n]各个数字出现的次��C��
*/
#include<iostream>
#include<cstring>
#include<map>
#include<algorithm>
#include<stack>
#include<queue>
#include<cmath>
#include<string>
#include<cstdlib>
#include<vector>
#include<cstdio>
#include<set>
#include<list>
#include<numeric>
#include<cassert>
#include<ctime>

using namespace std;

void gao(__int64 n,  __int64 cnt[]){//count all digits in [1,n] , not [0,n]
    /*
    check a form like this .
        000
        001
        .
        010
        .
        099
        100
        101   ------N
    */
    int digit[20], num = 0;
    __int64 p = 1, N = n, sub = 0;
    do{
        digit[num++] = n % 10;
        n /= 10;
        p *= 10;
    }while(n > 0);
    for(int i = num - 1,j ; i >= 0 ; i --){
        p /= 10;
        for(j = 0 ; j < digit[i] ; j++){
            cnt[j] += p;
            for(int k = 0; k < 10 ; k++){
                cnt[k] += p/10*i;
            }
        }
        N %= p;
        cnt[j] += N+1;//要加1
        sub += p;
    }
    cnt[0] -= sub;
}

int main()
{
#ifndef ONLINE_JUDGE
    freopen("in","r",stdin);
#endif

    for(__int64 left, right ; cin >> left >> right; ){
        __int64 cnt[10] = {0};
        gao(left-1, cnt);
        for(int i = 0 ; i < 10 ; i++){
            cnt[i] = -cnt[i];
        }
        gao(right, cnt);
        for(int i = 0; i < 10 ; i++){
            if(i){
                putchar(' ');
            }
            printf("%I64d", cnt[i]);
        }
        puts("");
    }
    return 0;
}

_Yuan 2011-03-23 13:32 发表评论

_Yuan — Tue, 22 Mar 2011 07:02:00 GMT

/*
    �l�出n�?�Q�m�?�Q�组成一个串�Q�要求这个串的�Q意前�~�都有1的个�?gt;=0的个�?br>    求满��L��串的个数
    n,m <= 10^6

    脑残丫，很明昄��“�q�栈出栈”�Q�居然没反应
    计数公式�?br>        (n+m, m) - (n+m, m-1)

            (n+m)!(n-m+1)
    =    ------------------
           (n+1)!m!

    O(n)的计��会��时�Q�要对N!分解质因�?br>    用到公式N!里质因子p的个��Cؓ N/p + N/p^2 +

    �q�要注意别忘记对(n-m+1)分解质因子了�Q�如果放在最后来乘它的话�Q�有可能�Ҏ��个质因子p�Q�只用那三个��到�?br>    �q�小�?�Q��?0100501又不是素敎ͼ��难��?nbsp;p^x % 20100501  , x < 0
    因�ؓ�Ƨ拉定理 a^x % n = a ^(x % phi(n) ) % n 需�?nbsp;(a,n) = 1�Q�！�Q�！
    八中 1856跟这题一��P��不过那里�?0100403是素敎ͼ�可以那样�?nbsp;p^x

    参�?nbsp;
    http://hi.baidu.com/matrush/blog/item/5298d8a3786b478546106447.html
*/
#include<iostream>
#include<cstring>
#include<map>
#include<algorithm>
#include<stack>
#include<queue>
#include<cmath>
#include<string>
#include<cstdlib>
#include<vector>
#include<cstdio>
#include<set>
#include<list>
#include<numeric>
#include<cassert>
#include<ctime>

using namespace std;

const int MOD = 20100501;
const int MAXN = 2000000;

bool isPrime[MAXN+10];
int pr[MAXN] , tot;

void init(){
    for(__int64 p = 2 ; p < MAXN ; p ++){
        if(!isPrime[p]){
            pr[tot++] = p;
            for(__int64 j = p * p ; j < MAXN ; j += p){
                isPrime[j] = true;
            }
        }
    }
}

/*
    1*2**n
    tot = n/p , n/p^2 ,
*/
int get(int p, int n){
    int tot = 0;
    while(n > 0){
        tot += n / p;
        n /= p;
    }
    return tot;
}

int ipow(int a , int n){
    if ( n == 0 ) {
        return 1;
    }
    if (n ==1 ) {
        return a % MOD;
    }
    __int64 ans = ipow(a, n/2);
    ans = ans * ans % MOD;
    if ( n&1 ) {
        ans = ans * a % MOD;
    }
    return ans;
}

int solve(int n, int m){
    if(n < m){
        return 0;
    }
    /*
        (n+m)!(n-m+1)
        ------------------
            m!(n+1)!
    */
    __int64 ans = 1;
    int nm = n - m + 1;
    for(int i = 0 ; i < tot && pr[i] <= (n+m) ; i++){
        int cnt = 0;
        while(nm % pr[i] == 0){
            nm /= pr[i];
            cnt++;
        }
        ans = ans * ipow(pr[i], get(pr[i], n+m) - get(pr[i], m) - get(pr[i], n+1) + cnt) % MOD;
    }
    return ans;
}

int main()
{
#ifndef ONLINE_JUDGE
    //freopen("in","r",stdin);
#endif
    init();
    int T;
    for (scanf("%d",&T); T-- ; ) {
        int n, m;
        scanf("%d%d", &n, &m);
        printf("%d\n", solve(n,m));
    }
    return 0;
}

_Yuan 2011-03-22 15:02 发表评论