uva 11762 dp

1159 단어
// UVa11762 Race To 1
// Rujia Liu
#include
#include
#include

const int maxn = 1000000 + 10;
double f[maxn];
int removed[maxn], primes[maxn], vis[maxn];

void gen_primes(int n) {
  int p = 0;
  memset(removed, 0, sizeof(removed));
  for(int i = 2; i <= n; i++) if(!removed[i]) {
    primes[p++] = i;
    if(i <= n/i) for(int j = i*i; j <= n; j += i) removed[j] = 1;
  }
}

double dp(int x) {
  if(x == 1) return 0.0; //   
  if(vis[x]) return f[x]; //    
  vis[x] = 1;
  double& ans = f[x];
  int g = 0, p = 0; //   g(x) p(x)
  ans = 0;
  for(int i = 0; primes[i] <= x; i++) {
    p++;
    if(x % primes[i] == 0) { g++; ans += dp(x / primes[i]); }
  }
  ans = (ans + p) / g;
  return ans;
}

int main() {
  int T;
  scanf("%d", &T);
  gen_primes(1000009);
  memset(vis, 0, sizeof(vis));
  for(int kase = 1; kase <= T; kase++) {
    int n;
    scanf("%d", &n);
    printf("Case %d: %.10lf
", kase, dp(n)); } return 0; }

책에 핵심 코드를 제시했는데, 정말로 나는 잘 이해하지 못했다.
출시: f[x]=(∑f[x/y]+p)/g(x);
dfs 또는 dp의 방법으로 f(x)를 계산하기;

좋은 웹페이지 즐겨찾기