栏目分类:
子分类:
返回
名师互学网用户登录
快速导航关闭
当前搜索
当前分类
子分类
实用工具
热门搜索
名师互学网 > IT > 面试经验 > 面试问答

zoj 3205 Derivative

面试问答 更新时间: 发布时间: IT归档 最新发布 模块sitemap 名妆网 法律咨询 聚返吧 英语巴士网 伯小乐 网商动力

zoj 3205 Derivative

#include <iostream>#include<stdio.h>using namespace std;#define pmod 1000000007#define LL long longLL a[105][105],s[105],c[105],d[105],e[105];LL ext_gcd(LL a,LL b,LL &x,LL &y){LL t,d;if(b==0){x=1;y=0;return a;}d=ext_gcd(b,a%b,x,y);t=x;x=y;y=t-a/b*y;return d;}LL Invmod(LL a){LL x,y;LL n=pmod;if(ext_gcd(a,n,x,y)!=1) return -1;return (x%n+n)%n;}int main() {int t,n,m,q;LL r,v,p;int i,j;scanf("%d",&t);while(t--){scanf("%d%d",&n,&m);for(i=0;i<n;++i){for(j=0;j<=m;++j)scanf("%lld",&a[i][j]);}scanf("%d",&q);while(q--){for(i=1;i<=m;++i)scanf("%lld",&c[i]);for(i=0;i<n;++i){e[i]=0;d[i]=a[i][0];s[i]=a[i][0];for(j=1;j<=m;++j){if(c[j]==0&&a[i][j]!=0)e[i]++;v=c[j];r=a[i][j];p=1;while(r){if(r&1){s[i]*=v;s[i]%=pmod;if(c[j]){d[i]*=v;d[i]%=pmod;}}v*=v;v%=pmod;r>>=1;}}}for(j=1;j<=m;++j){v=0;for(i=0;i<n;++i){if(c[j]==0){if(a[i][j]==1&&e[i]==1){v+=d[i];v%=pmod;}}else if(a[i][j]){p=1;p=(s[i]*a[i][j])%pmod;p=(p*Invmod(c[j]))%pmod;v+=p;v%=pmod;}}if(j!=1)printf(" ");printf("%lld",v);}printf("n");}if(t)printf("n");}return 0;}
转载请注明:文章转载自 www.mshxw.com
本文地址:https://www.mshxw.com/it/372947.html
我们一直用心在做
关于我们 文章归档 网站地图 联系我们

版权所有 (c)2021-2022 MSHXW.COM

ICP备案号:晋ICP备2021003244-6号