牛オフ小さなキャンディW練習試合48 C（数学、多項式）
リンク：https://ac.nowcoder.com/acm/contest/923/C出典：牛オフネットワーク
タイトル説明
小さなワットと彼の2人のチームメイトの帝都、それの終わり作られた何かで、皆のためのいくつかの材料のメリットを準備tokitsukaze？長い間考えた後、それらの3人がお菓子をたくさん用意します。彼らは、左から右へとn ..... 1,2,3,4,5,6,7,8,9に番号が付けから長い列に並んn個の子供たちをみましょう。
すべての3人のキャンディの髪は、特定の位置から始めて、キャンディの所有者のみが右側にこの男と、この男を送りました。しかし、彼らはお菓子のルールが異なって作られました。
お菓子で作られたホイールがtokitsukazeであれば1、、彼女は右の皆のためのお菓子を回し、位置posから開始します。
2は、お菓子で作られたホイールは帝都であれば、彼は右に回し、位置posから開始します、彼はキャンディ、彼は髪に会った二人を送って、彼が遭遇した最初の人を与えます3つのキャンディ...個人のk k個のキャンディ髪、とで構成され、彼が遭遇した第三者への2キャンディは、出会ったn個のために人々の権利数まで歩きました。
3.ホイールは、キャンディwinterzz1で作られたよく知られて小さなwは寛大で、彼は砂糖のほとんどを作ったので、彼は右に回し、位置posから始まりますされている場合は、それが彼の最初の出会いを提供しますキャンディヒトの毛髪は、第二者に彼が第三者に、キャンディ毛髪4に遭遇し、彼が遭遇した毛髪9 ... k番目個々キャンディ遭遇毛髪K2K ^ {2} K2までキャンディー彼は人々のn個の権利までを歩きました。
砂糖福祉で作られた開始位置posが糖人々で作られたホイールをMと、彼らは砂糖のラウンドをしたこと、みんなの手の中に砂糖菓子からなるn 1に私メートルのこのラウンドの終了を教えてくださいを教えてくれ、今メートルラウンドの合計を行い、数の数が、多すぎる避けるために、この数は、あなたが唯一の個々の糖の手のそれぞれの出力を使用して109 + 7modを国防省; 10 ^ {9} + 7mod109 + 7缶後の結果。
説明を入力します。
最初の行は、Tの整数正のテストケース（1⩽T⩽10）（1 \ leqslant T \ leqslant10）（1⩽T⩽10）、発現Tセットです。各ケースの場合：最初の2行目の正の整数N、M（1⩽n、m⩽105）（1 \ leqslant N、M \ leqslant10 ^ {5}）（1⩽n、m⩽105）について表しますn個のm個の個々のラウンドの行が作られ、キャンディました。2つの正の整数型（1⩽type⩽3）（1つの\ leqslantタイプ\ leqslant 3）（1⩽type⩽3）、POS（1⩽pos⩽n）（1 \ leqslant POS \の次のm行leqslant N）（1⩽pos⩽n）はホイール紙吹雪人物を示し、人物の位置は、キャンディを開始します。タイプ=キャンディ時に人工毛tokitsukaze 1、タイプ= 2人の帝都の紙吹雪、タイプ= 3人winterzz1紙吹雪。POSは位置を示し、ほとんどの人は、右端の男性がnであるPOS 1、POSを残しました。
出力説明：
对于每组测试案例，输出一行n个非负整数，表示每个人手中的糖果数量mod 109+7mod\; 10^{9}+7mod109+7后的结果。数字与数字之间用空格隔开并且行末不允许有多余空格。例1
エントリー
4
10 1
1 1
10 1
2 2
10 1
3 3
10 3
1 1
2 2
3 3輸出
1 1 1 1 1 1 1 1 1 1
0 1 2 3 4 5 6 7 8 9
0 0 1 4 9 16 25 36 49 64
1 2 4 8 14 22 32 44 58 74備考：
由于输入量和输出量比较大，请勿使用cin,cout进行输入输出。本题不会卡常数，不用特地使用输入输出挂。アイデア：
1.type = tokitsukaze人工毛のお菓子1は、砂糖で作られた以降のPOS場合、砂糖は関係なく、一人一人の場所の、すべての背後に+1を持っています。
2.type = 2帝都ヒト毛キャンディー、開始POSから糖からなる場合、第1のX（\（X> = POS \）位置を得る）\（X-POS + 1 \ ） 糖。
3.type = 3キャンディは、人工winterzz1た場合は糖からなる開始POSから、最初のx（x> = POS）得る個人\（（X-POS + 1 ）^ 2 \） 糖類、および\（（ X-POS + 1）^ 2
= X ^ 2 + 2 *（1-POS）X + POS * 2-2 * POS + 1 \） 見つけることができ、x番目の個々キャンディを多項式で表すことができる：
\ （F（X）= A_ { I} X ^ {2} + B_ {I} X + C_ {I} \）
また、それぞれが唯一の多項式の係数を変更し、糖からなる、見出すことができる\（AI、バイ、CI \ ）
すなわち：
ときタイプ= 1：位置posからのすべての\（C {X} \）を加え、すべての1つの
タイプ= 2：位置posからのすべての\（C {X} \）を加え、すべての\（（ 1-POS）\） 、全ての\（B {X} \）を加え、すべての1
タイプ= 3：全ての位置pos \（C {I} \）を加え、すべての\（（posが^ {2} -2 * POS + 1）\）すべて\（B {I} \）を加え、すべての\（（2-2 * POS）\） 、全ての\（{I} \）プラス1
限り、我々は違いスコアの係数を維持するとして挙げることができます。
コード：
#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <cmath>
#include <queue>
#include <stack>
#include <map>
#include <set>
#include <vector>
#include <iomanip>
#define ALL(x) (x).begin(), (x).end()
#define sz(a) int(a.size())
#define rep(i,x,n) for(int i=x;i<n;i++)
#define repd(i,x,n) for(int i=x;i<=n;i++)
#define pii pair<int,int>
#define pll pair<long long ,long long>
#define gbtb ios::sync_with_stdio(false),cin.tie(0),cout.tie(0)
#define MS0(X) memset((X), 0, sizeof((X)))
#define MSC0(X) memset((X), '\0', sizeof((X)))
#define pb push_back
#define mp make_pair
#define fi first
#define se second
#define eps 1e-6
#define gg(x) getInt(&x)
#define chu(x) cout<<"["<<#x<<" "<<(x)<<"]"<<endl
#define du3(a,b,c) scanf("%d %d %d",&(a),&(b),&(c))
#define du2(a,b) scanf("%d %d",&(a),&(b))
#define du1(a) scanf("%d",&(a));
using namespace std;
typedef long long ll;
ll gcd(ll a, ll b) {return b ? gcd(b, a % b) : a;}
ll lcm(ll a, ll b) {return a / gcd(a, b) * b;}
ll powmod(ll a, ll b, ll MOD) {a %= MOD; if (a == 0ll) {return 0ll;} ll ans = 1; while (b) {if (b & 1) {ans = ans * a % MOD;} a = a * a % MOD; b >>= 1;} return ans;}
void Pv(const vector<int> &V) {int Len = sz(V); for (int i = 0; i < Len; ++i) {printf("%d", V[i] ); if (i != Len - 1) {printf(" ");} else {printf("\n");}}}
void Pvl(const vector<ll> &V) {int Len = sz(V); for (int i = 0; i < Len; ++i) {printf("%lld", V[i] ); if (i != Len - 1) {printf(" ");} else {printf("\n");}}}
inline void getInt(int* p);
const int maxn = 100010;
const int inf = 0x3f3f3f3f;
/*** TEMPLATE CODE * * STARTS HERE ***/
const ll mod = 1e9 + 7;
struct node
{
ll a, b, c;
void update()
{
if (a < 0)
a += mod;
if (b < 0)
b += mod;
if (c < 0)
c += mod;
if (a >= mod)
{
a %= mod;
}
if (b >= mod)
{
b %= mod;
}
if (c >= mod)
{
c %= mod;
}
}
} info[maxn];
int t;
int n, m;
int op;
int pos;
int main()
{
//freopen("D:\\code\\text\\input.txt","r",stdin);
//freopen("D:\\code\\text\\output.txt","w",stdout);
du1(t);
while (t--)
{
du2(n, m);
repd(i, 1, m)
{
du2(op, pos);
if (op == 1)
{
info[pos].c++;
} else if (op == 2)
{
info[pos].b++;
} else
{
info[pos].a++;
}
}
ll ans;
repd(i, 1, n)
{
info[i - 1].a += info[i].a;
info[i - 1].b -= 2ll * (i - 1) % mod * info[i].a % mod;
info[i - 1].c += info[i].a * (i - 1) % mod * (i - 1) % mod;
info[i - 1].b += info[i].b;
info[i - 1].c -= info[i].b * (i - 1) % mod;
info[i - 1].c += info[i].c ;
info[i - 1].update();
info[i] = info[i - 1];
ans = ( (1ll * info[i].a * i % mod * i % mod + 1ll * info[i].b * i % mod) % mod + 1ll * info[i].c) % mod;
if (ans >= mod)
{
ans %= mod;
}
printf("%lld%c", ans, i == n ? '\n' : ' ');
}
repd(i, 0, n)
{
info[i].a = info[i].b = info[i].c = 0;
}
}
return 0;
}
inline void getInt(int* p) {
char ch;
do {
ch = getchar();
} while (ch == ' ' || ch == '\n');
if (ch == '-') {
*p = -(getchar() - '0');
while ((ch = getchar()) >= '0' && ch <= '9') {
*p = *p * 10 - ch + '0';
}
}
else {
*p = ch - '0';
while ((ch = getchar()) >= '0' && ch <= '9') {
*p = *p * 10 + ch - '0';
}
}
}