Loading

数学 - SGU 118. Digital Root

Digital Root 

Problem's Link


 

Mean: 

定义f(n)为n各位数字之和,如果n是各位数,则n个数根是f(n),否则为f(n)的数根.

现在给出n个Ai,求出A1*A2*…*AN + A1*A2*…*AN-1 + … + A1*A2 + A1 这个式子的数根.

analyse:

这道题目要用到这个规律,设f(n)是n的digital root,那么f(A*N)=f(A*f(N));

 

具体证明过程如下:

  设自然数N=a[n]a[n-1]…a[0],其中a[0],a[1]、…、a[n]分别是个位、十位、…上的数字

  再设M=a[0]+a[1]+…+a[n]

  求证:N≡M(mod 9).

 证明:
   ∵ N=a[n]a[n-1]…a[0]=a[n]*10^n+a[n-1]*10^(n-1)+…+a[1]*10+a[0].
   又∵ 1≡1(mod 9),
   10≡1(mod 9),
   10^2≡1(mod 9),
   …
   10^n≡1(mod 9).
   上面这些同余式两边分别同乘以a[0]、a[1]、a[2]、…、a[n],再相加得:
     a[0]+a[1]*10+…+a[n]*10^n≡(a[0]+a[1]+…+a[n])(mod 9),
                 即 N≡M(mod 9),得证。

  有了这个性质就容易解决本题了

  在计算过程中,可以不断mod 9,因为我们知道有这样两个性质:

 

    (A+B)mod C = ((A mod C) + (B mod C))mod C
    (AB)mod C = ((A mod C)×(B mod C)) mod C

 还要注意,如果余数为0,则输出9.

Time complexity: O(N)

 

view code

/**
* -----------------------------------------------------------------
* Copyright (c) 2016 crazyacking.All rights reserved.
* -----------------------------------------------------------------
*       Author: crazyacking
*       Date  : 2016-01-08-10.51
*/
#include <queue>
#include <cstdio>
#include <set>
#include <string>
#include <stack>
#include <cmath>
#include <climits>
#include <map>
#include <cstdlib>
#include <iostream>
#include <vector>
#include <algorithm>
#include <cstring>
using namespace std;
typedef long long(LL);
typedef unsigned long long(ULL);
const double eps(1e-8);


#define REP( i, n ) \
   for ( int i = 0; i < (n); i++ )
#define REPD( i, n ) \
   for ( int i = (n) - 1; i >= 0; i-- )
#define FOR( i, b, e ) \
   for ( int i = (b); i <= (e); i++ )

typedef long long int64;

const int MAXN = 1000;

int T, N;
int64 val[MAXN];

int droot( int x )
{
   if ( x < 10 ) return x;
   int ans = 0;
   while ( x > 0 )
   {
       ans += x % 10;
       x /= 10;
   }
   return droot( ans );
}

int main()
{
   scanf( "%d", &T );
   while ( T-- )
   {
       scanf( "%d", &N );
#warning READ LLD
       REP( i, N )
       scanf( "%I64d", &val[i] );
       int64 ans = droot( val[N - 1] );
       REPD( i, N - 1 )
       ans = droot( droot(val[i]) * droot(ans+1) );
       printf( "%I64d\n", ans );
   }

   return 0;
}

 

posted @ 2016-01-08 11:06  北岛知寒  阅读(176)  评论(0编辑  收藏  举报