牛客寒训营4 G(梅开二度，欧拉降幂又一次不会)

最新推荐文章于 2024-10-09 22:28:58 发布

想出成果的acmer

最新推荐文章于 2024-10-09 22:28:58 发布

阅读量331

点赞数

分类专栏：牛客文章标签：算法 c++ 线性代数

本文链接：https://blog.csdn.net/xianqiuyigedao/article/details/122837492

版权

牛客专栏收录该内容

53 篇文章 0 订阅

订阅专栏

题目
题意: 给定长度为n的数组，求出该数组所有非空子序列的权值的乘积。权值： 该数组的最大值最小值
思路:
首先，观察到每个数都会在非空子序列中出现2^(n-1)次。
之后，发现对于第i小的数，他有2^ (i-1)作为最大值出现，2^(n-i)次作为最小值出现。
sort之后，枚举每个数统计贡献即可。但是发现幂次是2^n，即使用快速幂也要n * logn * n，寄。
根据费马小定理， a^(p-1) % p == 1(a和p互质)。尽管幂次很大，但我们可以去%(p-1)，使得从而将幂次控制在1e9以内。 a ^x, x = k(p-1) + left. a ^ x % p = 1 * a^left.
时间复杂度: O(nlogn)
代码:

// Problem: 子序列权值乘积
// Contest: NowCoder
// URL: https://ac.nowcoder.com/acm/contest/23479/G
// Memory Limit: 524288 MB
// Time Limit: 2000 ms
// 
// Powered by CP Editor (https://cpeditor.org)

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<complex>
#include<cstring>
#include<cmath>
#include<vector>
#include<map>
#include<unordered_map>
#include<list>
#include<set>
#include<queue>
#include<stack>
#define OldTomato ios::sync_with_stdio(false),cin.tie(nullptr),cout.tie(nullptr)
#define fir(i,a,b) for(int i=a;i<=b;++i) 
#define mem(a,x) memset(a,x,sizeof(a))
#define p_ priority_queue
// round() 四舍五入 ceil() 向上取整 floor() 向下取整
// lower_bound(a.begin(),a.end(),tmp,greater<ll>()) 第一个小于等于的
// #define int long long //QAQ
using namespace std;
typedef complex<double> CP;
typedef pair<int,int> PII;
typedef long long ll;
// typedef __int128 it;
const double pi = acos(-1.0);
const int INF = 0x3f3f3f3f;
const ll inf = 1e18;
const int N = 2e5+10;
const int M = 1e6+10;
const int mod = 1000000007;
const double eps = 1e-6;
inline int lowbit(int x){ return x&(-x);}
template<typename T>void write(T x)
{
    if(x<0)
    {
        putchar('-');
        x=-x;
    }
    if(x>9)
    {
        write(x/10);
    }
    putchar(x%10+'0');
}
template<typename T> void read(T &x)
{
    x = 0;char ch = getchar();ll f = 1;
    while(!isdigit(ch)){if(ch == '-')f*=-1;ch=getchar();}
    while(isdigit(ch)){x = x*10+ch-48;ch=getchar();}x*=f;
}
#define int long long
int n,m,k,T;
int a[N];
int f[N];
int qpow(int a,int k)
{
	int res = 1;
	while(k)
	{
		if(k&1) res = res * a % mod;
		a = a * a % mod;
		k >>= 1;
	}
	return res % mod;
}
void solve()
{
   read(n);
   for(int i=1;i<=n;++i)
   {
   	  read(a[i]);
   }
   sort(a+1,a+n+1);
   int ans = 1;
   f[0] = 1;
   for(int i=1;i<=n;++i)
   {
   	  f[i] = f[i-1] * 2 % (mod-1);
   }
   for(int i=1;i<=n;++i)
   {
   	  ans = ans * qpow(a[i],f[i-1]  ) % mod;
   	  ans = ans * qpow(a[i],f[n-i]  ) % mod;
   }
   ans = (ans + mod) % mod;
   write(ans);
}
signed main(void)
{  
   T = 1;
   // OldTomato; cin>>T;
   // read(T);
   while(T--)
   {
   	 solve();
   }
   return 0;
}