zoukankan      html  css  js  c++  java
  • bzoj 3209 数位DP+欧拉定理

    枚举1的个数,统计有那么多1的数的个数

     1 /**************************************************************
     2     Problem: 3209
     3     User: idy002
     4     Language: C++
     5     Result: Accepted
     6     Time:0 ms
     7     Memory:844 kb
     8 ****************************************************************/
     9  
    10 #include <cstdio>
    11 #define Mod 10000007
    12 #define Phi 9988440
    13  
    14 typedef long long dnt;
    15  
    16 dnt n;
    17 dnt dp[50][2][51];
    18 int tb;
    19  
    20 void dodp() {
    21     for( tb=49; tb>=0; tb-- )
    22         if( (n>>tb)&1 ) break;
    23     dp[tb][1][1] = 1;
    24     dp[tb][0][0] = 1;
    25     for( int b=tb-1; b>=0; b-- ) {
    26         for( int c=0; c<=tb-b+1; c++ ) {
    27             dp[b][1][c] = dp[b+1][1][c-((n>>b)&1)];
    28             dp[b][0][c] = (dp[b+1][0][c] + dp[b+1][0][c-1]) % Phi;
    29             if( (n>>b)&1 ) dp[b][0][c] = (dp[b][0][c] + dp[b+1][1][c]) % Phi;
    30         }
    31     }
    32 }
    33 dnt mpow( dnt a, dnt b ) {
    34     dnt rt;
    35     for( rt=1; b; b>>=1,a=(a*a)%Mod )
    36         if( b&1 ) rt=(rt*a)%Mod;
    37     return rt;
    38 }
    39 int main() {
    40     scanf( "%lld", &n );
    41     dodp();
    42     dnt ans = 1;
    43     for( int c=2; c<=tb+1; c++ )
    44         ans = (ans * mpow(c,dp[0][1][c]+dp[0][0][c])) % Mod;
    45     printf( "%lld
    ", ans );
    46 }
    View Code
  • 相关阅读:
    CodeForces gym Nasta Rabbara lct
    bzoj 4025 二分图 lct
    CodeForces 785E Anton and Permutation
    bzoj 3669 魔法森林
    模板汇总——快读 fread
    bzoj2049 Cave 洞穴勘测 lct
    bzoj 2002 弹飞绵羊 lct裸题
    HDU 6394 Tree 分块 || lct
    HDU 6364 Ringland
    nyoj221_Tree_subsequent_traversal
  • 原文地址:https://www.cnblogs.com/idy002/p/4557172.html
Copyright © 2011-2022 走看看