穿上马甲就没认出polya定理之典型题目手环

链接：https://www.nowcoder.net/acm/contest/75/C
来源：牛客网

最近古yu迷上了《恋与制作人》，天天嚷嚷着白起我老公(死gay死gay的),白起又向古yu提出了一个问题：
给你n种卡牌（数量无限）,将其摆成长度为m的圆环的方案数。

(如果这个环可以通过若干次旋转或者反方向读(（abc）(cba)是一样的),则认为他们是一样的)
这时古yu一脸茫然，大哭特苦不能守护白起了。古yu希望你能帮助他解决这个问题去守护白起。

这里的n就相当于polya定理的m种颜色。。。

#include<iostream>
#include<cstdio>
#include<cstring>
#include<string>
#include<iomanip>
#include<vector>
#include<map>
#include<algorithm>
#include<cmath>
using namespace std;
typedef long long ll;
const ll mod=1000000007;
ll quickmod(ll a,ll b,ll c)
{
    a=a%c;
   ll res=1;
    while(b)
    {
        if(b&1)
            res=(res*a)%c;
        a=a*a%c;
        b>>=1;
    }
    return res;
}
ll gcd(ll a,ll b)
{
    if(b==0)
        return a;
    return gcd(b,a%b);
}
int main()
{
    ll t,n,m,i;
    ll res;
    scanf("%lld",&t);
    while(t--)
    {
        scanf("%lld%lld",&n,&m);
        res=0;
        for(i=1;i<=m;i++)
        {
            res=(res+quickmod(n,gcd(m,i),mod))%mod;
        }
        if(m%2==0)
            res=(res+((quickmod(n,m/2,mod)+quickmod(n,m/2+1,mod))*m/2)%mod)%mod;
        else
            res=(res+(quickmod(n,m/2+1,mod)*m)%mod)%mod;
        //res=res/(2*m);
        printf("%lld\n",res*quickmod(2*m,mod-2,mod)%mod);//除法求余转化为求逆元
    }
}

代码中注释掉的res=res/(2*m);是这样想的，resm每一步res肯定比mod小了，在除以2*m肯定比mod小，那就不用再求余了，

这是错的！！！本来3/5和1/5不一样