ACM训练日记—4月12日

       这次主要总结下重点的题目吧。

CodeForces 227A

题意：给出A,B,C三点坐标，求A到B到C如何转弯。（直走，左转，右转）

        直接按题意写也可以，但是巨麻烦。直接用叉乘（叉积），两个向量叉乘还的向量，但旋转方向不同所得方向不同（正负不同）。

Codeforces 615D Multipliers 

题意：给定m个质因子p[]，有p[1]*p[2]*...*p[m] == n，问n所有因子的乘积 % 1e9 + 7。

 讲解代码来自：https://blog.csdn.net/sinat_30062549/article/details/50489302

小于200000的素数不是很多，所以肯定有很多重复的数字。开一个map，记录素数和对应的个数。

然后看组合情况，假设素数j有k个，那么j有（kj+1）种出现的形式，所以对于所有的素数就有（k1+1）×（k2+1）×（k3+1）×（k4+1）*...*（kn+1）= sum种组合方式。不考虑j时就有sum/（kj+1）种。假设其中一种除j外的乘积为a，那么就有（a）×（a×j）×（a× j^2）×（a×j^3）×...×（a×j^kj） = （a^(kj+1)）×（j^(kj*(kj+1)/2)） .所以在所有的情况中共乘了j^((kj*(kj+1)/2))^(sum/(kj+1)).对于每一个素数求解一遍，乘积mod（1e9+7）为结果。

typedef long long ll;  
ll p[200050];  
map<int ,ll>MA;  
const ll mod = 1e9+7;  
ll fast_mod(ll n,ll m,ll mod){  
    ll ans =1;  
    while(m){  
        if(m&1)ans = ans*n%mod;  
        m>>=1;  
        n = n*n%mod;  
    }  
    return ans;  
}  
int main(){  
    int T;  
    cin>>T;  
    MA.clear();  
    for(int i = 0;i<T;i++){  
        int a;  
        scanf("%d",&a);  
        MA[a]++;  
    }  
    map<int,ll>::iterator it;  
    ll sum = 1;  
    ll ans = 1;  
    for(it = MA.begin();it!=MA.end();it++)

    {  
        int cnt = fast_mod(it->first,(it->second+1)*it->second/2,mod);  
        ans=fast_mod(ans,it->second+1,mod)*fast_mod(cnt,sum,mod)%mod;  
        sum=sum*(it->second+1)%(mod-1);  
    }  
    printf("%I64d\n",ans);  
}  

另外还有一种用费马小定理做的

来自：https://blog.csdn.net/V5ZSQ/article/details/50751878

#define maxn 22222
#define mod 1000000007ll
typedef long long ll;
int n,num[10*maxn],p[maxn],res;
ll l[maxn],r[maxn];
ll mod_pow(ll a,ll b,ll p)
{
    ll ans=1ll;
    a%=p;
    while(b)
    {
        if(b&1)ans=(ans*a)%p;
        a=(a*a)%p;
        b>>=1;
    }
    return ans;
}
int main()
{
    while(~scanf("%d",&n))
    {
        res=0;
        memset(num,0,sizeof(num));
        while(n--)
        {
            int temp;
            scanf("%d",&temp);
            if(!num[temp])p[++res]=temp;
            num[temp]++;
        }
        l[0]=r[res+1]=1ll;
        for(int i=1;i<=res;i++)l[i]=l[i-1]*(num[p[i]]+1)%(mod-1);
        for(int i=res;i>=1;i--)r[i]=r[i+1]*(num[p[i]]+1)%(mod-1);
        ll ans=1ll;
        for(int i=1;i<=res;i++)
        {
            ll a=1ll*num[p[i]]*(num[p[i]]+1)/2%(mod-1)*l[i-1]%(mod-1)*r[i+1]%(mod-1);
            ans=ans*mod_pow(p[i],a,mod)%mod;
        }
        printf("%I64d\n",ans);
    }
    return 0;
}