codeforces 852F Product transformation

      题意：最开始数组s所有数字都等于a，每次执行是s[i]=s[i]*s[i+1]，s[n]不变，求经过m次操作后s[i]每个元素的值，答案对q取模。

      首先由于最开始每个数字指数是1，两个数字相乘就相当于指数相加，那么我们令dp[i][j]表示执行i次操作后第j个数字的指数，将数组反过来我们就能得到dp[i][j]=dp[i-1][j-1]+dp[i-1][j]的转移方程了，这不是和组合数很像嘛，但是本题的边界是所有dp[0][j]=1，我们打个表找个规律，就会发现，最终dp[m][j]就是组合数c[i][j]的前缀和，即dp[m][1]=c[m][0],dp[m][1]=c[m][0]+c[m][1]。由于需要o(n)求出来，所以我们直接求出来m这一行的值即可。但是我们求的是指数，所以在求组合数的时候不能直接取模，由欧拉定理我们便可知我们如果要对指数那么模数必须为 ，题目又给了如下限制, 所以我们可以通过暴力枚举求出这个作为指数的模数即可。

      下附AC代码。

#include<iostream>
#include<stdio.h>
#include<string.h>
#include<algorithm>
#define maxn 1000125
using namespace std;
typedef long long ll;
ll n,m,a;
ll fac[maxn],inv[maxn],ans[maxn];
ll quickpow(ll p,ll k,ll mod)
{
	ll ans=1;
	while(k)
	{
		if(k&1)
			ans=(ans*p)%mod;
		p=(p*p)%mod;
		k>>=1;
	}
	return ans;
}
ll c(ll x,ll y,ll mod)
{
	if(x<y) return 0;
	if(x==y || y==0) return 1;
	return (fac[x]*inv[y]%mod*inv[x-y])%mod;
}
int main()
{
	ll mod,temp=1,mo;
	scanf("%I64d%I64d%I64d%I64d",&n,&m,&a,&mod);
	for(int i=1;i;i++)
	{
		temp=(temp*a)%mod;
		if(temp==1)
		{
			mo=i;
			break;
		}
	}
	fac[0]=1;
	for(int i=1;i<=m;i++)
	{
		fac[i]=(fac[i-1]*i)%mo;
		inv[i]=quickpow(fac[i],mo-2,mo);
	}
	temp=0;
	for(int i=1;i<=n;i++)
	{
		temp+=c(m,i-1,mo);
		ans[i]=temp;
	}
	for(int i=n;i>=1;i--)
	printf("%I64d ",quickpow(a,ans[i],mod));
}