spoj422 Transposing is Even More Fun（polya）

现在存储了一个$2^a * 2^b$的矩阵
• 矩阵在内存中是按行存储的
• 现在你想求它的转置
• 唯一允许的操作是交换两个内存位置的值
• 求最少需要的次数？
• 4e5 组询问，每组询问 a + b <= 1e6

首先，我们可以把内存中每个位置编码为一个 a + b 位的01串
那么内存 s 处存储的值的目标位置是 s 循环右移 b位对应的内存位置
设这些位置有 k 个循环节，那么答案就是$2^{a+b}- k$。（每个长度为len的循环节我们只需要交换len-1次）

那我们现在就是要求对于$2^{a+b}$个不同的01串，在循环右移b位算同构的情况下有多少种是本质不同的。设$g=gcd(a+b,b)$，则置换群的个数为$\frac{a+b}{g}$.

这样其实可以看做用$2^g$种颜色给长度为n的一个环染色，旋转同构，求方案数。就是老套路了。复杂度$O(n+T\sqrt n)$

#include <cstdio>
#include <cstring>
#include <algorithm>
#include <cmath>
using namespace std;
#define ll long long
#define ld long double
#define inf 0x3f3f3f3f
#define N 1000010
#define mod 1000003
inline char gc(){
    static char buf[1<<16],*S,*T;
    if(S==T){T=(S=buf)+fread(buf,1,1<<16,stdin);if(S==T) return EOF;}
    return *S++;
}
inline int read(){
    int x=0,f=1;char ch=gc();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=gc();}
    while(ch>='0'&&ch<='9') x=x*10+ch-'0',ch=gc();
    return x*f;
}
int n,m,bin[N],prime[N],tot=0,phi[N];
bool notprime[N];
inline void init(){
    bin[0]=1;notprime[1]=1;phi[1]=1;
    for(int i=1;i<=1e6;++i) bin[i]=bin[i-1]*2%mod;
    for(int i=2;i<=1e6;++i){
        if(!notprime[i]) prime[++tot]=i,phi[i]=i-1;
        for(int j=1;prime[j]*i<=1e6;++j){
            notprime[prime[j]*i]=1;
            if(i%prime[j]==0){phi[prime[j]*i]=phi[i]*prime[j];break;}
            phi[i*prime[j]]=phi[i]*phi[prime[j]];
        }
    }
}
inline int gcd(int x,int y){return y?gcd(y,x%y):x;}
inline int ksm(int x,int k){
    int res=1;for(;k;k>>=1,x=(ll)x*x%mod) if(k&1) res=(ll)res*x%mod;return res;
}
inline int solve(){
    int x=1,res=0;
    for(x=1;x*x<n;++x)
        if(n%x==0) (res+=(ll)phi[n/x]*bin[m*x]%mod+(ll)phi[x]*bin[n/x*m]%mod)%=mod;
    if(x*x==n) res+=(ll)phi[x]*bin[m*x]%mod,res%=mod;return (ll)res*ksm(n,mod-2)%mod;
}
int main(){
//  freopen("a.in","r",stdin);
    int tst=read();init();
    while(tst--){
        int a=read(),b=read(),g=gcd(a+b,b);if(!a||!b){puts("0");continue;}
        n=(a+b)/g;m=g;int ans=bin[a+b]-solve();
        ans%=mod;if(ans<0) ans+=mod;printf("%d\n",ans);
    }return 0;
}