NOIP2015 求和

#include<cstdio>
#define UP(a,b,c) for(int a=(b);a<=(c);a++)
using namespace std;
typedef long long LL;
const int N=100010,MOD=10007;
LL n,m,a[N],c[N],num[N][2],si[N][2],sn[N][2],ss[N][2],x;
int main(){
    scanf("%lld%lld",&n,&m);
    UP(i,1,n)scanf("%lld",&a[i]);
    UP(i,1,n)scanf("%lld",&c[i]);
    UP(i,1,n){
        num[c[i]][i&1]++;
        si[c[i]][i&1]+=i%MOD;
        sn[c[i]][i&1]+=a[i]%MOD;
        ss[c[i]][i&1]+=a[i]*i%MOD;
    }
    UP(i,1,m)UP(j,0,1)
        if(num[i][j]>1)x=(x+(si[i][j]*sn[i][j])%MOD+(num[i][j]-2)*ss[i][j])%MOD;
    printf("%lld\n",x);
    return 0;
}

  

 

 

 

 

 

#include<cstdio>
#define FOR(a,b,c) for(int a=(b);a<=(c);a++)
using namespace std;
typedef long long LL;
const int maxn=100000+10;
const int MOD=10007;
LL n,m,a[maxn],c[maxn],sr[maxn][2],sa[maxn][2],s[maxn][2],sq[maxn][2];
int main(){
    scanf("%lld%lld",&n,&m);
    FOR(i,1,n)scanf("%lld",&a[i]);
    FOR(i,1,n)scanf("%lld",&c[i]);
    LL ans=0;
    FOR(i,1,n){
        int p=i&1,ci=c[i];
        s[ci][p]++;
        
        if(s[ci][p]>1)
             ans+=(a[i]*sr[ci][p])%MOD+(i*sa[ci][p])%MOD+sq[ci][p]+((s[ci][p]-1)*i*a[i]%MOD);
        ans%=MOD;
        
        sr[ci][p]+=i;
        sa[ci][p]+=a[i];
        sq[ci][p]+=i*a[i];
    }
    ans%=MOD;
    printf("%lld\n",ans);
    return 0;
}

  

转载于:https://www.cnblogs.com/qilinart/articles/5988551.html

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值