这题好难，我不会我不会啊~~~~

递归+快速幂，复杂度N³log²K

这..pb_ds T得好惨...还是我打开方式不对?

。

看到题，直接用的multimap，结果无情炸掉。。。后来发现直接vector就过了。。。呵呵哒

完全不管会不会溢出的longlong瞎写就能过。。。。

身败名裂
赋值时竟把方位数组搞混一位

最优解方向什么的写的好烦→_→

回复 @Hzoi_hzoier :
Orz

回复 @赵日天 :
哥们，你这算法是叫状压DP，不是数位DP。

回复 @HYOI_残星誓言 :
我没看过那么丑的代码，我的代码风格很好，很完美。

数学很重要

#include<iostream>
#include<cstdio>
using namespace std;
class node{
public:
long long int b,col,num;
};node d[100005]={0};
long long a1[2][100005]={0},n1[2][100005]={0},color[2][100005]={0},number[2][100002]={0};
long long int ans=0;
int n,m;
void cin1()
{
for(int i=1;i<=n;i++)
d[i].b=i;
for(int i=1;i<=n;i++)
cin>>d[i].num;
for(int i=1;i<=n;i++)
{
cin>>d[i].col;
if(d[i].b%2==0)
{
a1[0][d[i].col]=d[i].b+a1[0][d[i].col];
number[0][d[i].col]=d[i].num*d[i].b+number[0][d[i].col];
n1[0][d[i].col]=d[i].num+n1[0][d[i].col];
color[0][d[i].col]++;
}
else
{
color[1][d[i].col]++;
a1[1][d[i].col]=d[i].b+a1[1][d[i].col];
number[1][d[i].col]=d[i].num*d[i].b+number[1][d[i].col];
n1[1][d[i].col]=d[i].num+n1[1][d[i].col];
}
}
}
int main()
{
freopen("2015sum.in","r",stdin);
freopen("2015sum.out","w",stdout);
cin>>n>>m;
cin1();
for(int j=0;j<=1;j++)
for(int i=1;i<=m;i++)
if(color[j][i]>1)
{
ans=(color[j][i]-2)*number[j][i]+a1[j][i]*n1[j][i]+ans;
}
cout<<ans%10007;
}

只DFS（1） 居然过了 O(∩_∩)O~
ps 不要吐槽代码风格，某人调乱的@Coolkid

双倍经验见usaco架设电话线

原来注释也耗时间！！
我带注释的代码194行5.562 s
删掉注释后代码115行4.204 s

合并类动归就是要注意区间！！！
合并类动归就是要注意区间！！！
合并类动归就是要注意区间！！！

这个数据精度要求不高，和BZOJ上保持一致

splay水题