标签:int sum ans long 解题 EXNR 拱门 include seed
大力统计题
考虑把和的平方拆开,最终就是许多对位置乘起来求和。所以考虑每对位置的贡献,对于$a_{i,j}$和$a_{k,h}(1<=i<=k<=n,1<=j<=h<=m)$,它们的贡献就是在$a_{i,j}$左上角和$a_{k,h}$右下角选两个点构成一个矩形把两个位置框起来的方案数再乘上对应系数,所以一般来说其贡献为
$2*a_{i,j}*a_{k,h}*i*j*(n-k-1)*(m-h-1)$
用前缀和维护二维乘积的和$sum[i][j]$,就可以更新答案了:
$sum[i][j]=a[i][j]*i*j-sum[i-1][j-1]+sum[i-1][j]+sum[i][j-1]$
$ans+=2*a[i][j]*sum[i][j]*(n-i+1)*(m-j+1)$
但是平方项的系数都是1,所以对于每个位置$i,j$还要扣掉一个$a[i][j]*i*j*(n-i+1)*(m-i+1)$
还没完......我们上面为了好写默认$1<=i<=k<=n,1<=j<=h<=m$,然后就成功地把$i<=k,j>h$和$i>j,j<=h$的情况给漏掉了,所以还要把每行翻过来再算一次,注意这次没有等于号了
1 #include<cstdio> 2 #include<cstring> 3 #include<algorithm> 4 using namespace std; 5 const int N=5010,mod=998244353; 6 long long n,m,ans,seed,a[N][N],sum[N][N]; 7 long long get() 8 { 9 seed=(seed*9301+49297)%233280; 10 return seed%101; 11 } 12 int main() 13 { 14 scanf("%lld%lld%lld",&n,&m,&seed); 15 for(int i=1;i<=n;i++) 16 for(int j=1;j<=m;j++) 17 { 18 a[i][j]=get(); 19 sum[i][j]=(a[i][j]*i*j%mod+mod-sum[i-1][j-1]+sum[i-1][j]+sum[i][j-1])%mod; 20 ans+=2*a[i][j]*sum[i][j]%mod*(n-i+1)*(m-j+1)%mod; ans%=mod; 21 ans+=mod-a[i][j]*a[i][j]%mod*i*j%mod*(n-i+1)*(m-j+1)%mod; 22 } 23 for(int i=1;i<=n;i++) 24 for(int j=1;j<=m/2;j++) 25 swap(a[i][j],a[i][m-j+1]); 26 for(int i=1;i<=n;i++) 27 for(int j=1;j<=m;j++) 28 { 29 sum[i][j]=(a[i][j]*i*j%mod+mod-sum[i-1][j-1]+sum[i-1][j]+sum[i][j-1])%mod; 30 ans+=2*a[i][j]*sum[i-1][j-1]%mod*(n-i+1)*(m-j+1)%mod; ans%=mod; 31 } 32 printf("%lld",ans); 33 return 0; 34 }View Code
标签:int,sum,ans,long,解题,EXNR,拱门,include,seed 来源: https://www.cnblogs.com/ydnhaha/p/10439495.html
本站声明: 1. iCode9 技术分享网(下文简称本站)提供的所有内容,仅供技术学习、探讨和分享; 2. 关于本站的所有留言、评论、转载及引用,纯属内容发起人的个人观点,与本站观点和立场无关; 3. 关于本站的所有言论和文字,纯属内容发起人的个人观点,与本站观点和立场无关; 4. 本站文章均是网友提供,不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属;如您发现该文章侵犯了您的权益,可联系我们第一时间进行删除; 5. 本站为非盈利性的个人网站,所有内容不会用来进行牟利,也不会利用任何形式的广告来间接获益,纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。