Pの競プロ記

競技プログラミングに関することを書きます。

Codeforces 236B Easy Number Challenge

概要


問題文はこちら
http://codeforces.com/problemset/problem/236/B


解法


素数判定・約数の個数判定がきちんとできていれば解ける問題です。

素数判定について
私はエラトステネスの篩を用いて求めました。蟻本にも載っているので参考にしてみてはいかがでしょうか。

約数の個数判定について
(※これの求め方を忘れてたので復習の必要あり)
参考になる記事: こちら

N = Ap * Bq * Crのとき、
約数の個数は、(p + 1) * (q + 1) * (r + 1) 個。


ソースコード

#include <iostream>
#include <string>
#include <algorithm>
#include <vector>
#include <queue>
#include <stack>
#include <map>
#include <set>
#include <numeric>
#include <complex>
#include <sstream>
#include <cstdio>
#include <cstdlib>
#include <cmath>
#include <cctype>
#include <cassert>
#include <cstring>
using namespace std;

#define REP(i,x,n) for(int i = x; i < (int)(n); i++)
#define rep(i,n) REP(i,0,n)
#define f first
#define s second
#define mp make_pair
#define pb push_back

typedef pair<int,int> P;
typedef long long ll;

static const int INF = (1 << 29);
static const double EPS = 1e-9;
static const double PI = acos(-1.0);

static const int MOD = 1073741824;

bool flag[1000001];
vector<int> prime;

void eratos(int n){
  for(int i = 0 ; i <= n ; i++)flag[i] = false;
  for(int i = 3 ; i <= n ; i+=2)flag[i] = true;
  prime.push_back(2);
  flag[2]=true;
   
  for(int i = 3 ; i <= n ; i += 2){
    if(!flag[i])continue;
    prime.push_back(i);
    for(int j = i + i ; j <= n ; j+=i)flag[j]=false;
  }
}

int calc(int val){
  int ret = 1;
  for(int j = 0; j < prime.size(); j++){
    int count = 0;
    if(val < prime[j]) break;
    while(val != 0 && val % prime[j] == 0){
      val /= prime[j];
      count++;
    }
    ret *= (count + 1);
  }
  return ret;
}

int main(){
  eratos(1000000);
  int a,b,c;
  cin >> a >> b >> c;
  int sum = 0;
  rep(i,a) rep(j,b) rep(k,c) {
    sum += calc((i + 1) * (j + 1) * (k + 1));
    sum %= MOD;
  }
  cout << sum << endl;
  return 0;
}