SGU 294 He's Circle(Polya原理+高精度)

最新推荐文章于 2019-01-11 01:00:23 发布

原创最新推荐文章于 2019-01-11 01:00:23 发布 · 1.1k 阅读

0 ·

本内容遵循CC 4.0 BY-SA版权协议

标签

#SGU #Polya原理 #高精度

Thinking More 同时被 2 个专栏收录

31 篇文章

订阅专栏

Polya原理、Burnside引理

10 篇文章

订阅专栏

该博客探讨了如何运用Polya原理解决SGU 294题目，即计算本质不同的涂色方案数，其中涉及到对项链涂色问题的分析，考虑了旋转的同构性，主要关注高精度计算的应用。

题目链接：
SGU 294 He’s Circle
题意：
用两种颜色涂n颗珠子的项链，问本质不同的方案数有多少？考虑旋转。
分析：
很容易根据 $Polya$ 原理得到：

A n s = 1 n \sum d | n n 2 d ϕ (n d)

$Ans=\frac{1}{n}\sum_{d|n}^{n}2^d\phi(\frac{n}{d})$
但是因为

n≤200000 $n \leq 200000$ ,并且答案不需要模数。那么难点就在于高精度了。
一开始我把

2 $2$ 的幂，每步都计算一下，乘以

2 $2$ 存在

data $data$ 数组里，可是这样子太慢了，想到了一个加速的方法： 每十次计算一下，相当于每次计算

210 $2^{10}$ ,也就是乘以

1024 $1024$ ,这样既不会数据溢出【因为我

ans $ans$ 里每个数是存

105 $10^5$ 以内的】，也保证了加速【整整十倍】。
还要注意最后的除以

n ，一开始是想每次二分找商，后来发现完全没必要，直接除就好了~~

#include <stdio.h>
#include <cstring>
#include <algorithm>
#include <cmath>
#include <ctime>
using namespace std;
typedef long long ll;
const int MAX_N = 200010;
const int Base = 10000;
const int SIZE = 100010;

int prime_cnt, prime[MAX_N], phi[MAX_N];
int data[SIZE], len, ans[SIZE], ed, n;

void GetPhi()
{
    prime_cnt = 0;
    phi[1] = 1;
    for(int i = 2; i < MAX_N; ++i){
        if(phi[i] == 0) {
            prime[prime_cnt++] = i;
            phi[i] = i - 1;
        }
        for(int j = 0; j < prime_cnt && i * prime[j] < MAX_N; ++j) {
            if(i % prime[j]) {
                phi[i * prime[j]] = (prime[j] - 1) * phi[i];
            }else {
                phi[i * prime[j]] = phi[i] * prime[j];
                break;
            }
        }
    }
}

inline void add(int x)
{
    int i = 0;
    ll carry = 0;
    while(carry || i < len) {
        carry = carry + ans[i] + data[i] * x;
        ans[i] = carry % Base;
        carry /= Base;
        i++;
    }
    ed = i;
}

inline void mul(int x)
{
    int carry = 0, i = 0;
    while(i < len || carry) {
        carry = carry + data[i] * x;
        data[i] = carry % Base;
        carry /= Base;
        i++;
    }
    len = i;
}

inline void div()
{
    for(int i = ed, carry = 0; i >= 0; i-- ) {
        carry = carry * Base + ans[i];
        ans[i] = carry / n;
        carry %= n;
    }
}

int main()
{
    GetPhi();
    scanf("%d", &n);
    //double begin = clock();
    data[0] = 1, len = 1, ed = 0;
    int step = 0;
    for(int i = 1; i <= n; ++i) {
        step++; 
        if(step == 10) { //这里就是加速！
            mul(1 << step);
            step = 0;
        }
        if(n % i) continue;
        mul(1 << step);
        step = 0;
        add(phi[n / i]);
    }
    div();
    for(int i = ed; i >= 0; i--){
        if(ans[i] != 0) {
            printf("%d", ans[i]);
            ed = i - 1;
            break;
        }
    }
    for(int i = ed; i >= 0; i--) {
        printf("%04d", ans[i]);
    }
    printf("\n");
    //printf("time elapsed : %.6lf Ms\n", clock() - begin);
    return 0;
}