洛谷10月月赛R1·浴谷八连测R1·提高组 SAC E#1 - 一道中档题 Factorial

版权声明:hhhhh https://blog.csdn.net/WT_cnyali/article/details/78171779

题意:求n!k进制下后缀0的个数。

老夫掐指一算,竟然掐到std算法?

  • Solution
    讲道理,刚开始看到这道题是不会做的。
    但是老夫想想如果m=10还是会做的。
    数数小于等于n有多少个5的倍数就行了。2的倍数绝对比5的倍数多嘛。
    于是我就有了一点想法。
    首先质因数分解少不了的。m=ki=1scntii
    然后,对于每个质因数si,算一算[1,n]能找到多少它的倍数(不过类似计算2时遇到23这些要算三次的)
    怎么计算?容斥一下?应该只有我这种zz一开始这么想?
    后来发现计算的时候每次加1就行了,因为前面的都已经计算过了。
    算完之后,发现si[1,n]中找到了xi,但是出现一个m至少si要出现cnti次。所以最终答案就是minki=1{xicnti}
#include <bits/stdc++.h>
#define LL long long
using namespace std ;
LL n, m, tot ;
LL s[1010], cnt[1010] ;
int main() {
#ifndef ONLINE_JUDGE
    freopen ( "a.in", "r", stdin ) ;
    freopen ( "a.out", "w", stdout ) ;
#endif
    LL i, j, x, ans, gn ;
    while ( scanf ( "%lld%lld", &n, &m ) != EOF ) {
        memset ( s, 0, sizeof s ) ;
        memset ( cnt, tot = 0, sizeof cnt ) ;
        for ( x = m, i = 2 ; i*i <= x ; i ++ ) {
            if (x%i == 0) s[++tot] = i ;
            while (x%i == 0) cnt[tot] ++, x /= i ;
        }
        if (x != 1) s[++tot] = x, cnt[tot] = 1 ;    
        /*printf ( "%lld = ", m ) ;
        for ( i = 1 ; i < tot ; i ++ ) printf ( "%lld^%lld + ", s[i], cnt[i] ) ;
        printf ( "%lld^%lld\n", s[tot], cnt[tot] ) ;*/
        ans = 1LL<<60 ;
        for ( i = 1 ; i <= tot ; i ++ ) {
            gn = 1 ;    
            x = 0 ; 
            for ( j = 1 ; gn <= n ; j ++ ) {
                gn *= s[i] ;
                x += n/gn ;
            }
            //printf ( "time %lld = %lld\n", s[i], x ) ;    
            ans = min(ans, x/cnt[i]) ;  
        }
        printf ( "%lld\n", ans ) ;  
    }
    return 0 ;
}

恬不知耻的发个广告:我的博客

阅读更多
想对作者说点什么? 我来说一句

没有更多推荐了,返回首页